1 


Université Ibn Tofail 
Faculté des Sciences - Kenitra 


Cours de 

PHYSIQUE STATSTIQUE 

Filière : Science de la Matière Physique (SMP) 

- Ss- 


2014-2015 


Pr . Ahmed QACHAOU 
Laboratoire de Physique de la Matière condensée (LPMC) 


exosup.com 


page facebook 


2 


Table des Matières 


Ch 1 : Ensembles et distributions à l’équilibre 

Ch 2 : Description des mélanges statistique 

Ch 3 : Distributions statistiques à l’équilibre thermique 

APPLICATIONS 

Ch 4 : Réversibilité microscopique et irréversibilité macroscopique 
Ch 5 : Distribution des vitesses de Maxwell 

Ch 6 : Systèmes de particules identiques et indépendantes.- Distributions 
de Maxwell-Boltzmann, de Fermi-Dirac, ET de Bose-Einstein. 

Ch 7 : Eléments de théorie de bandes d’énergie 

Ch 8 : Gaz de bosons en nombre non conservé- Distribution de Planck 


p : 1-11 

p : 12-25 
p : 26-40 


p : 41-63 
p : 64-75 

p : 76-94 
p : 95-116 
p : 117-128 


Note : 

Le cours dispensé dans ce polycopié correspond à la totalité du programme de la 
Licence de Physique ancienne accréditation (4ans). Il couvre et déhorde par conséquent 
le programme de la Licence Science de la Matière Physique nouvelle accréditation S y 

(3ans) -SMP -. 

Les parties en plus du programme de SMP sont des compléments ollrant aux 
étudiants volontaires une vue sur certaines applications du langage et formalisme de la 
physique statistique et permettant aussi une possibilité de lecture plus large. 


Kenitra 24/1/2014 


Pr. A.Qachao u 


exosup.com 


page facebook 


3 


Chapitre 1 : 

ENSEMBLES ET DISTRIBUTIONS STATISTIQUES A L’EQUILIBRE 


Introduction : Il existe trois ensembles statistiques correspondant à trois distributions 
statistiques : 

- Ensemble microcanonique ; 

- Ensemble canonique ; 

- Ensemble grand canonique. 

Dans la suite de ce chapitre, nous allons définir et étudier ces différents types de distributions 
et aussi introduire l’outil de raisonnement de physique statistique à l’état d’équilibre sur un 
problème d’ordre général tel que la marche au hasard qui représente un exemple simple de 
phénomènes naturels dépendant d’un grand nombre de paramètres. 

I- Description de la marche au hasard 

Soit un fou qui se déplace linéairement sur un axe x’ox (fig), on admet que : 

• Ses pas sont de même longueur unité. 

• Il effectue un mouvement linéaire et discret : Chaque pas est tantôt à droite (d) 
tantôt à gauche (g); son déplacement est donc irrégulier. 

• Soit n pas, dont k pas à droit et (n-k) pas à gauche =>son déplacement est 
erratique. 



1°) Définitions : 

a- Le déplacement m au bout de n pas est : m = k-(n — k)=2k-n 

b- Une complexion est une suite de n pas complètement déterminée (n, k et l’ordre des 

pas sont connus). 

Exemple : dggdcldgddggggg 

c- Une configuration est une suite de n pas où les nombres k et (n-k) des pas (d) et (g) 
sont connus, mais l’ordre des (d) et (g) n ’ est pas connu. 

2°) Remarque : Deux suites de n pas ne sont pratiquement jamais identiques =>le 
déplacement m n’est pas le même. Cependant m reste toujours petit devant le déplacement 
maximum (m max ) possible pour n très grand : C’est cette stabilité relative de m qui peut être 
décrite par une méthode statistique. Il suffit d’attribuer une mesure de probabilité à m. La 
stabilité de m résulte alors de la loi des grands nombres. 

II- Probabilisation de la marche au hasard 

Probabilisation possible sous conditions : 

• Caractère aléatoire de la marche 

• Etablissement d’une loi de probabilité. 
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1°) Caractère aléatoire de la marche 

Un pas n’est pas toujours dans la même direction du pas précédent : deux suites de n pas ne 
sont presque jamais identiques. 

La même expérience ne conduit pas au même résultat ! Le résultat est dû au hasard. 

Ensemble statistique : 

D’après ce qui précède, il est ainsi plus facile de raisonner sur le cas de n fous supposés tous 
‘'semblables” et faisant chacun un pas (au lieu d’un seul fou avec n pas). Leur marche sera 
bien entendu différente (il y aura toujours des petits détails différents d’un fou à l’autre). On 
retrouve cette situation dans toute expérience réalisée sur un système macroscopique 
dépendant d’un grand nombre de paramètres dont seul un petit nombre d’entre eux sont 
reproduits dans la préparation d’une expérience. On considère donc n expériences 
‘'semblables” (c’est à dire, préparées de la ‘ même manière’). Elles constituent un ensemble 
statistique donnant n résultats pas tous identiques. 

2°) Loi de probabilité 

On utilise la méthode des probabilité à priori, c.a.d., à partir de l’information dont on 
dispose sur la marche du fou on forme un ensemble complet de possibilités (toute 
complexion appartient à une et une seule possibilité) pour lesquelles il n’y a aucune raison 
qu’elles aient des probabilités différentes. Cette équiprobabilité des possibilités détermine la 
loi de probabilité. 

Exemple : Un dé non pipé et lancé différemment à chaque jet, n’a aucune raison a priori que 
les 6 résultats possibles ne soient pas équiprobables (p = 1/6 ). 

Propriété : La loi de probabilité dépend de l’information dont on dispose. 

Cette information est en général très incomplète (jusqu’ à quel point le fou a-t- il perdu le 
sens?,.. ). Le calcul des valeurs moyennes des grandeurs observables est possible une fois la 
loi de probabilité est choisie. Nous étudions quatre lois de probabilité correspondant à 
différents types d’informations possibles sur la marche du fou . 


a- Absence totale d’information : (loi binomiale) 

a ) Loi de probabilité 

L’ignorance sur la marche est totale. Il n’y a pas de raison de supposer que le fou a 
tendance à se déplacer plus vers la droite que vers la gauche. On suppose donc 
l’ équiprobabilité des pas (d) et (g) : p d = p =1/2 (p+q =1 avec p = p d et q = p g ). Ceci est 

équivalent au fait que toutes les complexions possibles C n sont équiprobables => La 
probabilité d’une complexion est : P(C n ) =l/W(n) . Dans une suite de n pas il y a donc : 
W{n) = 2" complexions possibles =>P(C (1 ) = 2“". De même la probabilité d’une configuration 

(k, n-k) est : P n (. k ) = C k n p k q"' k ou : 

C k n est le nombre de fois d’arrangements de k pas à droite parmi n total de pas. p 1 et q" k 
sont des probabilités des k pas à droite et ( n - k ) pas à gauche respectivement. On généralise 
à p et q quelconques. 


n+m n+m n—m 

Q„{m)=C n 2 p 2 q 2 
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0 ) Valeurs moyennes et fluctuations 

Définitions : i) La valeur moyenne d’une variable aléatoire k est : k = Y, kP ^ C n) 

C n 

ii)Le s fluctuations de la variable aléatoire k , notées (j(k) , est : 

<y{k)= yj(k-k') = ^ k 2 - k j 


Applications : i) Valeur moyenne du nombre k de pas à droite :k=^ j kP(c n )=^ j kP n {k) = 

c n k = 0 


t<P k r k =P—Y J C nP k <l"- k =P—{p + ‘iy=np(p + ‘l)"- 1 « k=np 
k= 0 dp k=0 dp 


Soit k = n/2 pour p- 1/2 ; p et q étant considérés comme variables indépendantes. 

ii) Valeur moyenne de m: m = 2k- n= n(p-q) => m = 0 pour p-q - 1/2 

iii) Les fluctuations cr(k) de la variable aléatoire k: 


Considérons 

p 2 #t£c„W* 
dp à'-o 

alors on a 

r) 2 00 

: Vz c .v<r‘ 

< p t7i 1 

00 

X ' 7 2 s~ik k n—k 

2/ C nP d 

k = 0 

00 

X ’ / k k n—k / 2 i 

~L kC nP Q =k -k 

k = 0 

. Or, 

r) 2 00 

2 u 'SP r^k k r 

P -^rL C nP <7 
dp' k = o 


-p 2 -/p(p^r- 

k 2 -k = p 2 n(n- 1 ) et a(k)= (k 2 - k ) - J{p 2 n{n -l) +np - n 2 p 2 ) = yjnpq ; p + q = 1. 


Par conséquent : cr(m) = 2cr(&) = 2 pipq . Pour p-q-l/2=> cr{k) - et a(rn)-\fn . 


Déf: Dispersions relatives. Les fluctuations ou dispersions relatives d’une variable 

<j(m) <j(m) 

aléatoire m sont définies par : (ou quelques fois — = - .) 


m„ 


m 


Convergence de la loi de probabilité : cr(m)de m est très petite devant m max dès que n est 
cx(m) 1 

très grand : — - — - = — — n >r -> 0 . On dit alors que la loi de probabilité est convergente. 

"Vax Vf! 


y ) L’ entropie S 

L’absence totale d’information sur la marche du fou décrit l’existence du plus grand 
désordre possible dans la suite des pas. Cependant une prévision précise est obtenue sur le 
déplacement m : celui-ci est presque sûrement négligeable devant m max possible, quand n est 
grand. 

Définition : Le manque d’information ou encore le degré de désordre, est mesuré par 


la quantité : 


S = LogW («) 


appelée entropie. 


Une définition équivalente : 


S= - LogP(C n ) = - ^ P (C„)LogP(C„) . 

C n 


Commentaire : 

Plus la marche est mal connue, plus W(n)est grand et plus le désordre est grand. 
L’entropie est maxima quand toutes les complexions sont possibles et également probables; 
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alors: S-knLog2. L’entropie est proportionnelle au nombre n de pas: C’est une 

caractéristique générale de l’entropie. 

S 



Cas contraire : Le fou est guéri et prend le chemin le plus cours! Tous ses pas sont 
dans la même direction. C’est l’état le plus ordonné possible. L’information sur cet état est 
complète. Il n’y a plus qu’une complexion possible (W (n) = l) et l’entropie S est nulle. 

Cependant, on a jamais une information complète ; ce qui fait que la loi binomiale n’a pas un 
grand intérêt. T/information dont on dispose est partielle ; elle correspond par conséquent aux 
trois distributions statistiques micro canonique ; canonique et gand-canonique ci-dessous. 


b- Distribution microcanonique - Ensemble microcanonique 


a j Loi de probabilité 

Le déplacement m au bout de n pas est supposé connu => Seules les complexions C n 


correspondant à m sont possibles. Le nombre de ces C n possibles est: W(n,m ) =C „ 2 . Si 
absence de toute autre information, alors ces C n sont équiprobables et alors : 

loi de distribution micro canonique. 


P{C n )=W-'{n,m) 


P ) Entropie microcanonique 


n+m 

S n (. m)= LogW (: n,m)= Log C n 2 


= Log n ! - Log 



! -Log 



i 


f n + m ^ 

n + m n - m r 

^ n - m N 

n — m 

{ 2 J 

Log 

2 2 

l 2 J 

2 


Or n est grand et m n , on peut utiliser la formule de Stirling: Log n ! - n Log n-n. Alors : 
S„ ( m ) — n Log n-n - — + — Log 


La forme asymptotique de S: 

Log 


S n (m ) —n Log n - " ^ log 


f n + m’ 

n — m T 

n — m\ 

l 2 J 

2 a 

l 2 J 


n ± m 
2 y 


2 2 n 


1 TTl 

— Log n +Log I — ± — I . Donc, la forme asymptotique de l’entropie micro- 


canonique est: S n (m)- -n( p d Logp d + Pg Log Pg ) avec Pd = - + et Pg = - - . 

Remarque : On retrouve encore la proportionnalité de S et de n, soit: S n (ni) = nS x avec 

.S', = - ( Prl Log Pd + Pg Log Pg j entropie d’un seul des n pas conforme à la définition générale: S — — LogP où 
P est la probabilité de chaque possibilité (il y en a deux : p = p d ou p = p g ). 
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La définition : S — — LogP généralise la définition S — LogW valable lorsque les W possibilités sont 
équiprobables. 


c- Distribution canonique - Ensemble canonique. 


a } Loi de probabilité 

Le déplacement M au bout de N pas est supposé connu => Comme précédemment, les 
complexions C N correspondantes sont équiprobables. 

Considérons maintenant seulement les n premiers pas de chaque suite de N pas («« iV) . Soit 
C n une complexion correspondant au déplacement m au bout de n pas. Le nombre de 
complexions C N contenant la complexion C n est égal au nombre de manières de choisir les 
(N-n) autres pas; soit: W (N -n,M -m) . 

W(^N — n,M — m ) W (n' , m') 


La probabilité d’une complexion C„ est alors: P (C„ ) = - 


avec n’ = (N - n) et m’ = (M - m). Soit encore : 


W(N,M) W(N,M) 

n' + m ' / / 

P(C„)= C^/c, 


N+M 

2 


On peut écrire: P(C„ ) oc exp S n , ( m ') où S„, ( m ') = LogC n , 2 est l’entropie micro canonique de 
la suite complémentaire de (N - n) autres pas. 

Dans la limite: h<kN , S n ,(m ') est donnée par la forme asymptotique précédente. Elle est de 


la forme : S n ,(m')=f 


^ M — m—m'^ 
y N — n — n' j 


«N => m « N et (n + m) « N, soit m « n’. Alors: 

a S 


et elle peut être développée au 1 er ordre en m/n ' car n 

8S„ 


s A™') - S n ,{M)~ 


m - 


dm' 


En posant: (3 -- 


dm ' 


on a alors: P(C„)cc exp ( ~[3m ) 


Finalement, la loi de probabilité de la distribution canonique d’une complexion C n est: 


1 


P(C n )_ex P (-^) 


Définition : La constante de normalisation Z (ft) est appelée fonction de partition 
canonique. Elle est donnée par : 


z(/ 9 = 2 > 


- (3 m 


Explicitement, on a : Z(/?) = ^<? Pm =^W ( n,m)e ,3m = ^ C„ 2 e pm 

C n m m 

en passant vers la somme sur les valeurs de m. Soit encore : 

n+m n+m n—m n+m n+m n—m n+m n+m 

Z ( fi )- = I Cj- (e-> )- (e-> )— = I C.“ (e~ f )“ (‘ f )~ = £ ' C,“ (e" )“ C )" 


Par conséquent : 


Z{fi)= [e- p +e p )" « Z (p)= (2ch p) n 


0 ) Valeurs moyennes et fluctuations 

i) Cas où n«N : 

• Valeur moyenne du déplacement m: 
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m=Y. mP {C n )= 


Explicitement : 


me 


—J3m 


le 


-P m — 


Z(P)t Z{P)dPT n 

dLog ( 2chp ) 


1 dZ(p) 
z (p) dp 


<=> 


dLogZ(p) 


m — — 


dp 


m — —n - 


dp 


<=> 


m — — nt h p 


m 


• Fluctuations de m : 


= Z ,?ï2p ( c »)= vZ ,î?2 ^ 


(m) = 


m — m 


pm _ 1 5 n _p m _ 1 5 Z _ 5 

( 1 5Z) 

+ 

f 1 ÔZ 1 

z dp 1 z dp 2 dp 

[zdp J 

Izôp) 


Donc : 
Explicitement : 


a 2 (m)— m 2 -m = 

K ’ dp 


( 1 5Z A 


Z dp 


<=> 


2 ( \ dm 

a (m)- 

V ’ dp 


( m ) = 


ch P 


i • <j(m) 1 

• Fluctuations relatives: — d=-Z = 


m 


yfn shp 


->• 0 (loi des grand nombres). 


ii) Cas : 1« n «N : Etude de la probabilité P n (m) du déplacement m au bout de n pas: 
P n ( m )=W(n,m)P(C n ) oc cf" +m)/2 •C*?' +m ' )/2 ^ P n (rn)oze s " H+sAM ~ m) 

Pour n et n’ grands, cette loi présente un maximum très aigu dans le voisinage de m=m . La 
valeur moyenne m est donc obtenue en cherchant ce maximum m + ; solution de l’équation : 

d dS n ( m ) 55„, [M -m ) 

— [ S n (m) + S n , (M - m)l = 0 ^ — = . Car : m'=M-m = 

dm L J dm dm ' 


55,., 55, , 


dm ' dm 


Défi 


■ i 


4 H 


m 


•î'^- ; 

\ 

I.pH 


0 " M 


p=^ 

dm ' 


m 


caractéristique de la suite complémentaire de n’ = (N -n) autre 


dS 

pas p'=—±- 
dm 


caractéristique de la suite de n pas. 


ôS n f m ) dS , (m - m ) 

Rq: L’équation — - = — — exprime l’égalité: p=p' pouvant être fonnellement 

dm dm ' 

interprétée comme l’équilibre de la suite des n pas avec la suite complémentaire des (N - n) 
autres pas ( analogue à 1’ équilibre thermique de 2 systèmes matériels macroscopiques). 
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Etude de P n (m) au voisinage de m + - m : On développe S n (m ) et S n , (m')=S n , (M -m) 

^m — rnj 




jusqu’ au 2 eme ordre. Ainsi: S n (m)=S n ( m 'j + 


dm 


l -\ 1 ô 2 5„ 

I m -ml + - 


2! dm 2 




20 

De même : 5„, (m') = S„, (m') + — - 

v 7 dm 

P) O 

) + 5 1 ,, . ( m')=S n ('«)+ S n , (m'Y 
= S lt (m) + S n , [m ') 


1 > 2! dm' 2 


[ni'— m'j 


+ .... et 


dm 


dS n , 

dm' 


[tn—rn^ + 


m =m 

d 2 S„ 


dm" 


d 2 S„ 


dm ' 


[m—m^j +... 


1 

+ — 

d 2 S n 

d 2 S„. 


2 ! 

dm 2 

_ m - dm ' 2 

"m '-m ' 


[m-nij + ... 


Par conséquent : 


p n(m) = P n(™) 


-\ 1 
exp — 


~d 2 S n 

dm 2 

+ 3 2 5„, 

- dm ' 2 

m=m 

"m-m' _ 

| m-m j 


Conclusion : La loi de probabilité P n ( m ) est une loi normale ( appelée encore loi de 


Rq : D’après la forme asymptotique de S n (m) on a : 


d 2 S, , 


d 2 5„ 

d 2 S„, 


dm 2 

dm ' 2 

m=m 

"m '=m ' _ 


-n 1 



dm 


2 2 

— n —m n 
m=m 


pour n »1 


et 


d \' 
dm ’ 2 


-n' 1 1 

2 2 

— : n' —m' n N —n 


pour N » n . Donc : o- 


1 1 

1 

n N -n 


Il en résulte encore que les fluctuations relatives : 


CF 

m, 


(m) 1 

\fn 


H— >00 


->0 . 


y j Equivalence entre les distributions microcanonique et canonique 

• Entropie en distributions canonique S c : 

s c =- LogP ( C K ) = - ^ P ( C K ) L°8 p ( c „ ) = ~dÆj e ~ Pm [LogZ + fîm) • Soit: S c = LogZ + j3m 

c n ^ c n 

Explicitement : S c = nLog (2 ch P) - n fith fi 

• Entropie en distributions microcanonique S n (m) 

S n (m)^-n[p d Logp d + p g Logp g ] avec p dg =\ ± Y’ 

Prenons : m = — nth fi , alors on a : S n ( m — — nth fi) — nLog ( 2 ch fi ) —n/3thf3—S c . 

Conclusion : Pour n grand ; les fluctuations de m sont négligeables en distribution 
canonique. Il revient donc de même de considérer m fixe et égal à sa valeur moyenne 
m = m . On est donc ramené au cas de la distribution micro canonique qui est ainsi 
équivalente à la distribution canonique pour m-m . D’où l’égalité des deux entropies. 

.. a (m) 

Les fluctuations relatives sont aussi équivalentes : — - — > 0 

m 

max 

Rq: Il n’en est pas de même lorsque n est petit. Même pour m-m les fluctuations différent 
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considérablement : elles sont nulles en distribution micro canonique (<r m - fli n _> () > oj tandis 
qu’elles sont grandes en distribution canonique (°v = Jü) ■ 


d- Distribution grand- canonique - Ensemble grand canonique 

a ) Loi de probabilité 

Le déplacement M au bout de N pas est toujours supposé connu. Considérons le 
déplacement m au bout de n premiers mais quelconques de pas. Cette fois n n’est pas fixé : n 
et m varient donc d’une complexion à une autre. Le nombre de complexions C N contenant la 

complexion C n correspondant à m et n est toujours: W (N -n,M -m) . Sa probabilité est 


n ' + m ' / N + M 


encore : P(C n )= C n , 2 / C N 1 ce e s " {m} avec n'=(N-n ) et m'={M -m). 

Pour n<KN ,il faut développer maintenant S n , (m 1 ) = f{n\ m') S n , {m ') aux voisinages de M et 


N: S„. (m') = 


as, 


dn' 


dS„ 

— m 

-n'=N dm 

m' = M 


. Posons: a = - 

n'=N dn' 


et fl=^ 

- if=N dm' 

m'=M 


, on a : 


" n'=N 
m'=M 


r(c,)= 


(a, fi) 


exp («/? - fl m) 


: loi de probabilité de la distribution grand-canonique. 


Définition : La grande fonction de partition est définie par : Z G (a,/?) = ^V 


an—flm 


soit : 


Z G [a,ff) - ^ e an ^ e />m =^e an Z(fl) où Z (fi) = ^e fSm est fonction de partition canonique. 

c„ 

1 


Explicitement : Z G = Ve" (2chfl)‘ - - 

n 1-2 e a chfl 

fl ) Valeurs moyennes et fluctuations 

- Valeurs moyennes: 


(développement de — ). 

1-x 


1 


1 d 


• m = 'J'mP(C )= — ^ m e an = ^ 


y an-flm 


‘‘G C 


z G dfl- 


1 dZ G (a,fl ) 
Zg àfl 


<=> 


- dLogZ G (a,fl) 


m= — 


dp 


i =Z np { c )~H nQa 


n— fl m 


'‘G C 


1 dZ G (a,fl) 
Z n da 


<=> 


n— + 


dLogZ G (a,fl ) 


da 


Explicitement : m = - nth fl 


Fluctuations de m et n : a 2 (m)= - 


dm 

~dfl 


et 


2 ( \ dn 

a (n)= + — 
V ’ da 


->1 


<j(m ) cr ( 77 ) 

Exercice : Venner que les fluctuations relatives: — =-Z et — — 

m n 

Conclusion : L’équivalence avec la distribution canonique est évidente pour n-n. Cette 
équivalence est limitée aux valeurs moyennes. Il n 'en est pas de même pour les fluctuations. 
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Chapitre 2 : 

DESCRI PTI ON DES MELANGES STATI STI QUES 


I- Introduction 

On ne peut jamais avoir une information complète sur un système physique X qui 
dépend d’un grand nombre de paramètres. On peut juste connaître quelques grandeurs de ce 
système; telles que l’énergie total, le volume,... 

Définitions : i) Un microétat de Hou état microscopique est une complexion. L’information 
pour un microétat est complète. 

ii) Un macroétat de X ou état macroscopique est déterminé par les grandeurs 
macroscopiques ou observables de ce système. L’information pour un 
macroétat est par contre toujours incomplète. ! 

Exemple : Les valeurs de B et de M caractérisent un macroétat d’un sel paramagnétique. Un 
microétat de ce sel est par contre connu si l’on connaît l’état de spin de chaque ion. 

Propriété : Un macroétat correspond généralement à plusieurs (souvent à un grand nombre) 

de microétats => Ensemble statistique. 

Définition d’un ensemble statistique '. N systèmes semblables ne sont pas tous identiques ; ils 
sont distribués, à t donné, suivant les différents microétats possibles. Ces N systèmes 
constituent un ensemble statistique. Chaque système de cet ensemble est décrit par une 
méthode statistique. A chaque microétat possible on associe une probabilité. La loi de 
probabilité tiendra compte de l’information dont on dispose. 

Définition d’un mélange statistique : L’ensemble des microétats possibles est appelé : 
mélange statistique. On dit aussi que l’état macroscopique du système est une superposition 
incohérente d’états microscopiques. 


II- Opérateur densité - Statistique quantique 
A- Cas d’un état (ou complexion) pur : Etat complètement déterminé 
1°) Opérateur densité : 

Soit {|n„)j une base orthonormée de l’espace de Hilbert des états. L’état dynamique 
du système est décrit par un ket | y/) : 

• k) = 2X(0k> telle que: C„ (/)=(«„ \yr) 

n 

• Normalisation : imom 

n 

• Valeur moyenne d’une observable A(t) : (A(r)}= (ii/\a\i//) = 'Y j C* (t)C p (t) A np 

n,p 

avec A np =(u n \A\u p ). 
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L’évolution de l’état dynamique est régie par l’équation de Schrôdinger : 


ôt 


Définition : L’état dynamique du système peut être décrit d’une manière équivalente, à 
l’aide de l’opérateur densité p(t ) défini par: p{t) - i//(tfj(i//(t') de matrice densité: 


p np (t ) = (u„\p(t)\u p } = C n (t)C* (t) . p nn : population de l’état \u n ) et p np : cohérence 

entre les états | u n ) et | u^j . 

Propriétés : 

• Normalisation: Trp = Y d P m = Y\ C n (Of = 1 

n n 

• Valeur moyenne d’une observable A(t) : Tr(pA) 


(A(t))= ty\A\y/) = Y,C* ( t)Ç p ( t)A np =J^(ys\u n )(u p \y/)(u n \A\u p ) = J](u p \y)(y/\u n )(u n \A\u p 

n,p n,p n,p 

= H{ U p H \ A \ U p) = Tj{ U P |ÆI “»)(“» \ A \ U n) = Yj{ U P \f>^ A \ U n) = Yj{ U P \p A \ U n) = Tr(pA) 

n,p p n p p 

avec A np =(u n |a|« p ) . 

• L’évolution de l’état dynamique est régi par l’équation équivalente à 

.•* 8 P 


l’équation de Schrôdinger : 




En effet : 


ih ^^- {y/{t)\ + ^ ~ =X\yr){yr\- \y/)(y\X= Hp- pjf = [jf (t) , p]. 

• Herméticité: p + - p . 

En effet: p + = (| |) + =k)(^| = / 7 - 

• Propriétés de p(t) spécifiques du cas d’un état pur: p 1 = p et Trp 2 = 1 . 

2°) Application : Lumière complètement polarisée. 

Considérons un faisceau de lumière parallèle monochromatique (flpm) se propageant 
dans la direction k (trièdre orthononné : in , u 2 , k : voir figure). 



Le champ électrique E est transverse : f?(r) = u l a l cos cot + ii 2 a 2 (coscot + (p) . 

En notation complexe: £(f)= 'JÏRes où e^c^m +c 2 m je' 0 " où c 1 =—j= et c 2 =-^e ,<p 

v2 v 2 

L’intensité I de E décrit généralement une ellipse dans le plan ( u x ,u 2 ). Elle est définie par : 

1 pT | — ► . 2 1/2 2\||2 ||2 

I— ~J 0 kv) dt — — \a x +a 2 J — |cj| + |c 2 | 

La description de la polarisation, donnée ici, est purement classique. Cependant, elle 
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est équivalente à une description quantique et elle illustre bien le formalisme de l’opérateur 
densité. L’état de polarisation est décrit par le ket \y/) de composantes c, et c 2 , sur la base 

des deux états de polarisation rectilignes, suivant u x et u 2 : | i//) = c l u l +c 2 u 2 . 

Intensité : I = (i//\i//) (la normalisation de | y/) est ici I et non pas 1). L’opérateur densité : 


p{t) = |K 0 )M 0 |- Sa matriee est p = ( Pnp ) = ]- 




^a x a 2 e ’ 


; p np {t) = c„{t)c*{t). 


La matrice p est décomposée sur la base canonique des matrices carrées : 


1 


O Oô (0 L f 0 -A (1 O'' 


0 1 


et 


1 0 


0 -1 


telle que: 


p = |(l +C 7-r) 


où : cr r= cr, • r, + a 2 ■ r 2 + a 3 ■ r 3 . Par identification, on obtient : 
1 1 • 1 l 2 2 \ 

r x = —a x a 2 cos^j ; r 2 = —a x a 2 smcp\ r 3 = — (aj -a 2 J 


Définition : Les sont les composantes d’un vecteur r d’un espace abstrait à 3 

dimensions, r est appelé vecteur de Stokes. Il caractérise complètement l’état de polarisation 
du (flpm). 



PôleD : 
Pôle G : 


r 

a | a 2 , 

V 


K 



f 

a | ci , , 
v 


n 



Equation : cp=kn : 

Hémisphère supérieure : 
Hémisphère inférieure : 


Lumière polarisée circulairement à droite. 

-> Lumière polarisée circulairement à gauche. 

-> Polarisations rectilignes. 

Polarisations elliptiques droites. 
Polarisations elliptiques gauches. 


B- Superposition incohérente d’états (ou macroétat) quantiques purs : 

Etat incomplètement déterminé 


1°) Opérateur densité : 

Considérons un système physique pouvant se trouver dans n états possibles 
(indépendants entre eux) | y/ i )} affectés de probabilités {/?,}. 

Définition : L’ opérateur densité p(t) est défini par : 

Soit : p(t) = YjPiP: où : Pi = Wi)(Vi | • 

i 

La matrice densité : p np (t ) = ^ (u„ \ (//, (0) Pi ('/L (0 \ U P ) = X C m (0 Pi C P * (0 

i i 


n 

P(f) =Zk/(0)^(^i(0 
1=1 
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Propriétés : 








Normalisation : Trp = Y J P tm = X| c « (Of = 1 

n n 

Valeur moyenne : (A(t)}= Tr(pA) . 

L’évolution de l’état dynamique: ih-^-=\jK [t),p] ■ 

Herméticité: p + = p . 

p est défini positif : V|i/}: Ui\p\u\>0. 


Remarques : i) Il faut faire une distinction entre la superposition incohérente ci-dessus (décrite par 

J/: . 

P np (t ) =ÏC ni (t)p i C Jjj (f) ) et la superposition cohérente suivante : \y/')='L^p i \yi) qui est décrite par 

i i 

l’opérateur densité : p = 'Lp l p l + X JpÏpj \ ) Wj • 

En général, p * IL p, p s . Cependant, si on superpose un grand nombre d’états (/a) qui différent entre 


eux seulement par une phase aléatoire Aq. r . alors X i// 

i*j 


V i 


. , . . A <p- 

= X |y/)((/qe " =0 . La superposition 

i*j 


devient incohérente. 

ii) Dans un mélange statistique quantique, le caractère aléatoire d’une mesure apparaît à deux niveaux : 
il y a d’abord l’ indéterminisme de la mécanique quantique (le résultat d’une mesure est aléatoire, même 
pour un état initial parfaitement connu). Ensuite, l’information incomplète sur l’état du système conduit 
à un traitement statistique de tous les états possibles. 


2°) Applicationi : Lumière partiellement polarisée. 

On superpose plusieurs (flpm), de même longueur d’onde, et complètement polarisés 
tel que : et |i//)=^|i// i ) . En général, l’intensité lumineuse I n’est pas la 

i i 


somme des intensités partielles : I = (i// \ y/) = X ( V- 1 Vj ) = X A + X ( V- 1 V ) ■ L ’ interférence 

i,j i i*j 

(y/i | ¥ , ) n’est généralement pas nulle et / ^ ^ / . Si (i// ( | y/j) =0 , les flpm (i) et (j) ont des 


polarisations opposées. Si on superpose un grand nombre de (flpm), de même intensité /„ 
avec des phases <p t distribuées au hasard, alors : («, | '//,}=('„ e" p< et 

c n c p -0 . Dans ce cas aussi il n’y a pas d’interférence : 

I = n I 0 . On dit qu’il y a superposition incohérente. La polarisation résultante est décrite par 
l’opérateur densité : p = ^p , , . Chaque /?. est proportionnelle à l’intensité 7 0 . Celle-ci 


\ 

r \ 

Xki)W v w= 


i*j / 

v ‘*J / 


constitue le poids de chaque flpm dans la superposition incohérente. Elle joue donc le rôle de 
la probabilité p t . Le vecteur de Stokes résultant est : r = — V r t , barycentre des p affectés 

n i 

des masses 7 0 . La région intérieure de la sphère de Poincaré correspond à cette polarisation 
partielle, super-position incohérente de lumières complètement polarisées. La polarisation est 
nulle si r = 0 . 


Définition : le degrés de polarisation est défini par: P = |r| < 1 . 

Inversement, toute lumière partiellement polarisée peut être considérée comme la 
superposition de plusieurs (flpm) complètement polarisés, et ceci d’une infinité de manière. 
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3°) Application2: Faisceau de protons polarisés 

L’état de spin \y/) des protons, particules de spin 1/2 , est décrit sur la base ||+); |-)j de 
vecteurs propres de S. . 

a- Cas pur : Polarisation totale dans la direction u (0 , cp) : L’état de cette polarisation est 

■ . Q if. 

représenté par: = cos — e 2 1+) + sin— e 2 1-/ et l’opérateur densité 

p[6,(p) de matrice: p(9,(p) = 

\ A A AJ 

Le vecteur de polarisation ou vecteur de stokes est défini par: r{0 , <p) =Tr(pa) =u (0, <p) 

Valeur moyenne de spin : ( S^=^Tr[pcr ) =^r(0 , <p) . 


2 & 
cos — 


• O a _■ 
sm — cos — e v 


. e 0 

sm — cos — e 


+l(p 


sm — 


b- Superposition incohérent d’états purs représentée par : p = -^~ jdQP(0, cp)p(0, cp) 
P{0 , (p) est la probabilité de polarisation dans la direction u{0 , <p). Q est l’angle solide. 
Le vecteur de Stokes est maintenant : r =-^-J dQP(0, <p)r(0, cp ) . 


C- Partie d’un système. Trace partielle 


Soient les parties (1) et (2) d’un système global : (1) + (2). £, et £ 2 espaces des états 
correspondant : £ = £, ® £ 2 pour (1) + (2). Si{|w„)] est la base dans ce, et j|v p ^j la base dans 

E 2 , alors {|“„v p }| est base dans £ telle que : |m k v ;; )= \ U „)\ V P ) ■ 

A , observable agissant dans £, => A, son prolongement dans £ tel que : A, = A, <S> 1 2 
Valeur moyenne : ^A, ^j-Tr^pA l j . On peut écrire aussi: ^A, S j-Tr(p 1 A,) . 


p x = Tr 2 p : opérateur densité dans £, . p 2 = T\\ p : opérateur densité dans £ 2 . 

Résultat : En général p* p x ® p 2 . Ceci est dû aux corrélations entre les états des deux parties. 


D- Assemblée de spin 


2 


Soit un sel paramagnétique de n ions Cu 2+ . Soit | C n ) = | S zl S z2 S z3 S m } un état 


quantique ou complexion, possible de n ions Cu 2+ ; où 



. Les kets | C n ) , vecteurs 


propres de l’opérateur énergie magnétique 'H , constituent une base de l’espace des états à 2" 
dimensions. Cette assemblée de n spins sera alors dans une complexion | C„), panni les 2" 

avec une probabilité P(|C„)) , qui dépend de l’information dont on dispose. 

L’opérateur densité correspondant est : p=^|C K ) P(|C„))(C n | . 
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p est diagonale sur la base 1 1 C„ )} : 


p=—e 

Z 


-ptt 


avec : 


Z=Tr(e~ fix ) 


Exemple: Prenons le cas simple où il y a un seul ion en équilibre thennique avec un 


thermostat. La matrice densité correspondante est alors: p- — 

z 


( 0 -f>* 


\ 


0 

+/? s 


ou 


0 é- ^ 

e +Pe (il y a deux complexions possibles de spin pour cet ion : spin 
‘haut’ T et spin 'bas' i). 


z=e Pe 


Conclusion : Tous les états d’équilibre étudiés en physique statistique sont représentés par 
une matrice densité diagonale sur la base des vecteurs propres de l’hamiltonien ttf . 


III- Mélange statistique classique 
A- Etat complètement déterminé ou microétat classique 


Considérons le cas général d’un système ayant / degrés de liberté. L’état du système est 

décrit par 21 coordonnées généralisées q(q l ,q 2 ,q i , , q s ) et moments généralisés 

p{p l ,p 2 ,p i , , p, ) . Leur évolution est donnée par les 2 1 équations de 

* dO ’f * dJf 

Hamilton : q, = ; p j = ; î =1,2, ,/. jf{q, p,t) est le hamiltonien. L’espace des 

dPi £q, 

phases Q. est un hyperespace cartésien de dimension 2 / . L’état du système à un instant donné 
est représenté par un point M (q, p,t) e O. . L’unicité de la solution des équations de 

ÔJ~C 

Hamilton => par un point M (q, p,t) s Q , il passe une et une seule trajectoire. Si = 0 

dÔK 

=> tK (q, p,t ) est conservé (i.e. — = 0) => le point représentatif M se déplace sur une 

hypersurface : Jf=cste . Les trajectoires possibles sont sur cette hypersurface. 


B- Etat classique incomplètement déterminé (ou macroétat classique) 

l°) Densité en phase 

La loi de probabilité sur les différents microétats possibles devient en statistique classique une 
densité continue de probabilité p(q, p,t) dans l’espace des phases Q . p(q, p,t) remplace 
l’opérateur densité utilisé en statistique quantique. 

Définitions : La densité en phase p ( q, p,t) est définie par : 

p{q, p,t)dT = probabilité que M ( q ; p,t) e dT à l’instant t. dT -dqdp . 

2°) Propriétés : 

• a- Normalisation : J p (q, p,t ) dT = 1 

• b-Valeur moyenne d’une grandeur A(q, p,t ) : 
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( A (q, p,t)) = \ p (q, p,t)A(q, p,t)dT . 

• c- L’évolution est donnée par le théorème de Liouville: 

”Le volume Y d’un domaine de l’espace des phases Q reste constant 
lorsqu’on le suit au cours du temps 


Un point M (q, p,t) e Q. représente un système, donc un ensemble de points M, (q, p,t) e Q 
représente un ensemble de systèmes. Ainsi : 

‘ 'La densité p{q, p,t) des systèmes au voisinage d’un système donné dans l’espace 


des phases Q reste constante dans le temps ’ ’. Autrement dit : = 0 « ^ (f ) , /?] . 

Démonstration : 

En plus d’une démonstration basée sur les propriétés des transformations canoniques, 
on peut obtenir une démonstration équivalente en considérant Q comme analogue à un fluide 

incompressible de densité volumique p et de vitesse q, p \ . Alors l’équation de 

continuité s’écrit: + V |/n>j =0 (pas de création ni d’annihilation de particules : fluide 

dp — * — — * — 

incompressible) <^> — + V p • v + pV • v = 0 . Or, l’incompressibilité du fluide s’exprime par 

dt 

le fait que la vitesse v de ses particules ne diverge pas ; c'est-à-dire que : 

r . . a 


dp 




Ôq L+ Ôjq 

8q, dp, 


=1 


( d 2 tt 


d z K 


dp i dq i dqfip i 


= 0 . Donc : 


V^.Ï=Z <fq, 

i \ d q< 
Conséquences : 


f dp • dp • ^ 


dp, 


Pi 


=Z 


y 


r dp dtK dp d!K ' 

K dq t d Pj dp i dq i , 


dp 


[P>K] et: — - ] c.q.f.d. 


dp 


• p est unifonne: Vp- 0 => — = 0 <=> La densité des points représentatifs dans O est 

dt 

constante <=> Le nombre des systèmes dans un état donné est constant dans le 
temps L’ensemble des systèmes est en équilibre statistique. Alors : 

• Equilibre statistique <=>[/?,#"] = 0 <=> p= fonction des constantes du mouvement. 

• Equilibre statistique <=> Tous les états accessibles sont équiprobables. 


3°) Limite classique d’un mélange quantique : 

a- A la limite fi -» 0 ; la mécanique quantique est identifiée à la mécanique 
classique. La description d’un mélange statistique classique doit pouvoir être dérivé d’un 

mélange quantique. Si .4 (/;,/?,, r 2 , /y , , r v , p N ) est une fonction définie sur O , variant 

lentement et A est l’opérateur correspondant sur l’espace des états, obtenu en remplaçant 
p, par —V, ; la limite classique s’exprime par : TrA= Z< u n | A K) 

l n 

b- Normalisation de dY : Soit un gaz formé de N molécules identiques dont les 
degrés de liberté internes sont négligés : seuls interviennent les états de mouvement de 
translation des molécules. 
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Lemme : Le nombre d’états quantiques, ou complexions, dans un volume unité de 

l’espace des phases(aux relations d’incertitudes près) est asymptotiquement égal à : ^ 3jV 

quand ce nombre est très grand. Alors : V rflcO, d\ doit être normalisé tel que : 

dr = N \h w d " >]d ' P ' di, ' Nd3 Pn ■ 

• Le facteur h~ 3N s’explique par les relations d’incertitude de Heisenberg. Pour chaque d° 
de liberté on a : Ap i Aq i » h . Par conséquent, l’état de chaque molécule peut être représenté 
dans Q par une cellule de volume « h 3 . Si les distances intermolécules sont grandes en 
moyenne (A q t ») , on peut négliger le recouvrement de ces N cellules. La densité des 

complexions est donc » (l/h 3 ) W = h~ 3N . 

• résulte de l’indiscernabilité des molécules. Suivant le principe de Pauli, l’état 

quantique des N molécules doit être symétrique ou antisymétrique par rapport au groupe des 
permutations des N molécules entre elles. Il est symétrique pour des bosons et antisymétrique 
pour des fermions. Par contre en mécanique classique les molécules sont discernables, au 
moins par la position qu’elle occupent. A un ensemble de positions {r x ,p x ,r 2 ,p 2 , ,r N ,p N ^j il 

correspond N ! états distincts obtenus par permutations des molécules entre elles. En fait, ils 
correspondent tous à un même état quantique. Le volume d rcQ est donc toujours trop riche 
en complexions d’un facteur AM. Plus généralement, pour / degrés de liberté, l’élément de 

volume: 

comporte N { molécules de type(l) et N 2 molécules de type(2) ; S = N l ! N 2 ! . 


i 

n 

i = 1 


dq i dp i 

h 


doit être divisé par un facteur de symétrie S. Par exemple: si le système 


c- Condition de validité de la mécanique classique(MC) : 

"La longueur d’onde À. associée à chaque molécule doit être beaucoup plus petite que la 
distance moyenne d entre molécules : A <sc d ”. 


Soit pour un gaz avec d - 


V_ 


et A - 


h 


yj3mk B T 


(Â : longueur d’onde thermique des 


particules), la condition de validité de la M.C s’exprime par : 


N 

n - — n n 

V e 


3mk B T 


\ 3/2 


n est la concentration des molécules, n Q est leur "concentration quantique” et m leur masse. 

Il faut donc que la concentration n soit faible et la température T grande pour que la M.C 
soit valable. 


4°) Densité réduite : Comme pour la densité partielle en distribution quantique, on ne fait 

intervenir généralement que des observables à un corps et parfois à deux corps. On n’aura 
donc pas besoin de toute la densité p dans O mais seulement des densités en phase réduites à 

un corps f[r,p jet à deux corps f 2 (r ,p ,r' définies par: 
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f(r,p) = jjjp(r,p,r 2 ,p 2 , ,r N ,p N }dT N _, 

< 

f 2 ( r ,p,r',p''j =—jp(r,p,r',p',r 3 ,p 3 , ,r N ,p N jdT N _ 2 

Ces densités seront donc directement introduites dans les calculs des valeurs moyennes. 
Exemple : l’énergie interne U du système à N molécules identiques est : 


U = $f(r,p) 


(p-eA) 


2 m 


Vf) 


d 3 rd 3 p + — J / 2 (r ,p ,r' ,p'^w[^ r - r' j d 3 rd 3 pd 3 r'd 3 


P ■ 


IV- Complément : Entropie statistique 
A-Information statistique 


1°) Introduction : 

La marche au hasard fournit un exemple simple introductif à la théorie de 


l’information. La bijection : 


pas à droite : 
pas à gauche : 


<-> 1 
<-> 0 


constitue un code. Il s’ agit ici d’un code binaire, c’est à dire formé seulement de deux 
symboles (1 ,0). Ce doublet est appelé un "bit” (binary digit). 

Une suite de n pas est représentée dans ce code par un mot de n bits : 01 10101001 .. ..01. 
C’est un message codé pouvant être transmis par les différents moyens de communication 
(Radio, Téléphone,...). Exemple : Un mot du langage courant est codé par l’alphabet et peut 
être transmis par l’impression. 


2°) Définition de l’information statistique : 

a- Un mot de n bits permet de transmettre une quantité d'information I telle que : 
I = kLogW (n) où W (77) -2". 

b- Unités de I : En théorie de l’information, on pose k- — - — et 1 - Log 2 W (n) . 

Log 2 


Pour un mot de n bits : I = Log 2 2" = n => La quantité d’information est égale au nombre de 
bits du mot => ’Le bit est donc l’unité de l’information’. 

c- Nature de l’information : Avant d’observer la marche au hasard, il y a Vf (n)- T 

possibilités équiprobables pour une suite de n pas. Après observation, la suite des n pas est 
connue (complexion). Elle est représentée par un mot unique de n bits. Ceci conduit à deux 
points de vue équivalents : 

(i) Point de vue aposteriori : L’information apportée par l’observation réduit le nombre 
W (n)= 2" de possibilités à une seule (celle qui a été observée). Cette information est la plus 


grande possible. L’information serait moins grande si on n’observait que le nombre de pas à 
droite. D’une manière générale si W’ est le nombre de possibilité qui restent après 


W 

transmission du message codé, la quantité d’infonnation transmise est : I = Log 2 — 


(ii)Point de vue apriori : Avant l’observation des n pas, il y a W[n)= 2" possibilités 


équiprobables. On ignore complètement quelle sera la complexion observée. La quantité I 
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mesure ce manque d’information. Ce point de vue rejoint celui de l’entropie. C’est lui qui 
prévaut en physique statistique. 

d- Désordre statistique : 

Le désordre statistique traduit l’existence de plusieurs possibilités pour un événement 
donné. Le désordre croit avec le nombre de possibilités. L’information statistique est donc une 
mesure quantitative du désordre. 


e- Généralisation: Théorème de Shannon : 

W 

Si les W possibilités ont des probabilités p x , p x ,..., p w , alors: l(p x , p x ,..., p w )=-k'^ j p i Logp i . 

i = 1 


Si W possibilités sont équiprobables => 


> Pw) = 


\W’ W 



f- Propriétés de I 

i) Lemme : Si {<?, } et {/?, }, i e{l ,...., n} , sont deux suites de nombres réels (>0) telles 

que: Y J q i =TjPi alors : Za Lo £4, • 

i i i i 


Démonstration : il suffit d’utiliser Vx ; Logx< x-l. 

ii) I est maximum quand toutes les possibilités sont équiprobables. 

iii) I >0 ; sauf s’il y a une seule possibilité seulement, alors : I = 0 . 

iv) I croit avec W croissant ( I croit avec la longueur du message). 

v) Les évènements impossibles (p; = 0) ne contribuent pas à l’information I : 

l(pt,p 2 ,....,p w ,0,0 I(p it p 2 ,...., p w ) 

vi) I est symétrique . 

vii) Additivité séquentielle : Soient A et B deux séquences telles que : 

V W 

A = {p i ,p 2 ,....,p v }^>q A =YjPi et B = {Pv*’Pv + 2 ,--’P w }^>< 1 b= Z Pi alors: 

1=1 i=V+ 1 


l(pi,p 2 ,...,p v ,Pv ,-,Pw) = l{q A ,q B ) + Va 1 


r 

Px 

Pv 1 

r 

(q A ’ 

’ q a ) 

+ q B I 

V 


Pv + 1 
q,< 


Pw ' 
C 1b y 


Les informations successives s’ajoutent. 

7 = ~Y.Pi L °8Pi = ~dA Tj — l °S—- YjP> Lo sq A 

k i i q A q A i 


V 1 Pi T Fi 

q B Zj — L o s — 
q a 


Pi 


V+l 1 


Z 


v+i 


P,Logq B . 


viii) Corrélations (sous additivité) : (x i )'" et (y ; )" m et n possibilités. Chaque 
expérience réalise un couple ( y , y , ) avec une probabilité p ij . 

• x ; seuls sont réalisés avec une probabilité : p, PiLogp , . 

j i 

• Les y j seuls sont réalisés avec : P j = Z Pu 1 y = ~ k Z PjLogPj . 

i j 

• La somme des deux infonnations s’écrit : I x + I y - - k Z Pi j Logp i p j . 

ÿ 

• Aux mn possibilités (y, y, ) -> l’information : I xy = - ^ Z Pu kogp,, 

ij 

Or d’après le lemme => 7 n , <I X + I y . 

• Conclusions : 

1- Les corrélations entre x, et y. diminuent l’information. 

2- / -I x +1 si x,. et y j sont indépendants , alors p tJ = p l p ; . 
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3- Si la corrélation est complète, c.a.d si les jc. sont toujours liés à y. 

avec m = n, p tj - p ^ => / rv =I X , l’information sur une seule suite (x, )"' 
de possibilités suffit. 

ix) Mélange ( concavité) : Soit une urne contenant n boules rouges, l’information 
concernant la couleur des boules est nulle. De même pour une autre ume contenant n boules 
bleues, l’information concernant la couleur des boules est nulle. Mélangeons les deux urnes 
=> /= kLog 2 , pour une boule tirée, puisque les deux couleurs sont équiprobables. Le mélange 
des deux urnes a donc augmenté l’information manquante. 

Généralisation : Soit une urne contenant n x boules dont p l rouges et ( «, - p, ) bleues, les 
probabilités de tirage : p R - — et p B = — — — . Pour une deuxième urne de ru boules dont p 2 

n i 77j 

rouges et ( n 2 - p 2 ) bleues: q R =—\ q B =— — — . Pour le mélange des 2 urnes on a: 

n 2 n 2 

n 2 

r R =\p R + À 2 q R et r B =À l p B + Â 2 q B avec : 2, = — - — ; i = 1, 2 ; V 2 =1 .L’information statistique 

>h + n 2 i=1 

résultant du mélange des 2 urnes vérifie l’inégalité : I (r fi , r B ) > \I (p R , p B ) + X,/ [q R , q B ) . 


B- Entropie statistique d’un mélange statistique quantique 


1 °) Définition : i) Pour tout système physique, on définit une entropie, identique à 
l’information statistique sur ce système. 


Remarque : Du fait de la perturbation introduite lors des mesures précisant un état quantique (microétat) du 
système, le point de vue a posteriori de l’information apportée par ces mesures est donc de peu d’intérêt. On lui 
préfère le point de vue apriori, selon lequel l’entropie, ou information statistique, mesure l’information 
manquante sur le système. 

ii) Soit un macroétat représenté par p diagonal sur ||t')j : p ='^ j \f/p i (j\. L’entropie 

i 

statistique S(p ) ou quantité d’information manquante est: 

Von-Neumann (1927). 


s {p)^- k YjPi L °SP, = ~kTr(pLogp ) 


Dém. : (n\p\ n)= p nn = ^ (n | i)p t ( i \ n ) = p n => S (p) = - p ii Logp ii = -kTr ( pLogp ) . 

i i 

Remarques : 

• s(p) => l’entropie mesure le d° de mélange statistique ou encore le désordre de la 

préparation incomplète du système physique. Mais s(p)ne mesure pas le désordre dû à l’ indéterminisme 
quantique (une grandeur peut prendre plusieurs valeurs dans un état pur) puisque s(p)=0 pour un état pur 
( LogPj =0 pour p i =1). 

• 5 (p) n’est pas une grandeur mesurable expérimentalement (observable). 11 n’est pas 
représentée par un opérateur hermétique. 

• La définition de S (p) s’étend aussi pour n -> œ , ou pour Z -» J . 

2 °) Propriétés de S ( p) : 

• S (p)<- kTr (pLogp') pet p 1 deux opérateurs densités. 
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S est maxima si tous les états propres /') sont équiprobables. 

S (p) > 0 . (S = 0 pour un état pur) 

S ( p) est invariante pour toute transfonnation unitaire sur p : p' = U pU~ l 

telle que : UU~' =1 . L’entropie ne dépend pas de la base 

Sous-additivité: (1) et (2) parties d’un système telles que =Tr.p , alors: 


S(p)< S (p x ) + S (p 2 ) ; la corrélation entre (1) et (2) diminue le désordre. Si 


(1) et (2) sont indépendantes => 5 (/?) = S (p 1 ) + 5 (p 2 ) , on a alors additivité. 

• Concavité : Soient p l et p 2 deux opérateurs densités du même système, alors : 


S(A l p i + A 2 p 2 )> À i S(p i ) + À 2 S(p 2 ) avec 2, 2 >0; A, + Â 2 =1. Elle exprime 


que la superposition incohérente de 2 mélanges statistiques se traduit en 
général par une perte supplémentaire d’information. 

Exemple : 2 états purs => 5 (p x ) = 5 (p 2 ) = 0 . Alors que la superposition 
=> p=\p x + A 2 p 2 => S(/?)>0 . Il y a égalité si p x = p 2 et \ 2 >0 
• L’additivité séquentielle reste construite dans le temps. 

C- Entropie statistique d’un mélange statistique classique 

1°) Définition : S cl (p)- -k ^pLogpdY, si p représente bien un opérateur densité et si 

1 


dV comporte 


NI h 3 


complexions. 


2 °) Propriétés : S cl a les mêmes propriétés que S en mécanique quantique, sauf une seule: 
S cl n ’est pas bornée inférieurement. 

Considérons une particule se déplaçant sur Ox. Supposons que sa position xe [x,x + Ax] et 


son impulsion p e [p , p + Ap] . La densité: p = 


ApAx 


(car : pAT = 1= p ) e t 


S cl = -kpLogpAT = 


— k - 


h 


h 


Log 

ApAx ApAx 


ApAx 


<=> S,=-kLog 


h 


h ApAx 

Si ApAx < h o S cl < 0 . Mais p est alors en contradiction avec la relation d’incertitude de 
Heisenberg. Cet exemple montre que toutes les densités de probabilité ne sont pas observées 
dans 

la nature ; car certaines violent le principes de la mécanique quantique. Ces dernières 
correspondent à S pouvant être négative. 


Remarque: S cl n’est définie qu’à une constante additive ( —k log h ) près. 
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Chapitre 3 : 

DI STRI BUTI ON STATI STI QUE A L'EQUI Ll BRE THERMI QUE 

-I- Choix de l’opérateur densité 


1°) Rappel : 

• L’information sur un système macroscopique Z est incomplète: il y a un grand 
nombre de complexions (microétats) possibles. Ce nombre de complexions est limité 
par l’information partielle dont on dispose (qui peut être, par exemple l’énergie du 
système, le volume V, le nombre de particules, . . .). 

• Un ensemble statistique est formé de N exemplaires semblables de ce système. Ces N 
exemplaires sont répartis aléatoirement sur les complexions possibles avec une 
certaine loi de probabilité. 

• Le macroétat du système est alors un mélange statistique de complexions possibles 
affectées d’un poids statistique égal à la probabilité. Ce mélange statistique est décrit 
par l’opérateur densité p permettant le calcul des valeurs moyennes (statistique et 
quantique) de toute grandeur du système. Il suffit pour cela de déterminer la loi de 
probabilité. 

• Pour un système isolé à l’équilibre thennique, cette loi est fournie par le principe 
d’équiprobabilité . 

2°) Principe d’équiprobabilité : 

” Toutes les complexions possibles d’un système isolé à l’équilibre thermique sont 
équiprobables”. 

Ainsi si W est le nombre de complexions possibles, la probabilité d’une complexion est: W" 1 

3°) Définitions : 

a- L’ équilibre thermique d’un système isolé est un macroétat tel que toutes les 
grandeurs macroscopiques caractérisant ce système restent constantes dans le temps. 

b- Tout système isolé évolue vers un état d’équilibre thermique avec une constante de 
temps appelée : temps de relaxation, et ceci quelques soient les conditions initiales 
(sauf très rares exceptions). 

c- Un système isolé à l’équilibre thermique n’évolue plus macroscopiquement. Cet 
équilibre thermique est un macroétat stationnaire. 

Remarque : Cet équilibre macroscopique n’exclue pas l’évolution du microétat où se trouve le système 
à un instant t. Bien plus, à l’échelle microscopique, le microétat va évoluer de façon aléatoire à cause des 
interactions entre les différents constituants du système et aussi à cause des petites perturbations dues à l’action 
du milieu extérieur (ME), si le système n’est pas très exactement isolé. 

4°) Hypothèse érgodique : 

Le comportement stationnaire à l’échelle macroscopique et fluctuant à l’échelle 
microscopique est précisé par l’hypothèse érgodique constituée de trois conditions : 

i) le caractère aléatoire de l’évolution du microétat ; 

ii) la loi de probabilité apriori et, 

iii) la propriété de la moyenne temporelle de toute grandeur A(t) du système : 
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A = 


lim — 

A '->“ A t 


t + At 

| A(t^dt 


A(t) est un signal aléatoire dû aux fluctuations microscopiques ; sa valeur moyenne 

1 t + At 

temporelle (valeur moyenne dans le temps) est définie par: lim — f A(t)dt . 

Aï->co ' 

La résolution des appareils de mesure étant faible (finie), toute mesure de la grandeur A(t) est 

une moyenne temporelle : lim — J A(t)dt sur un intervalle A t grand devant la durée des 

fluctuations de A(t) . Celles-ci se compensent mutuellement et la moyenne temporelle ci- 
dessus est constante et est égale à la moyenne statistique A de A(t). 


Remarques : 

i) La résolution finie des appareils de mesure => même si le système contient un petit 

nombre de constituants, la mesure de A(t) doit être égale à A ; bien que les fluctuations deviennent 
importantes à cette limite ! 

ii) Dans le cas d’un système physique interagissant avec un ME, le principe d’équiprobabilité est 
appliqué au système global V u (me\ supposé en équilibre thermique. 

La méthode statistique => trois cas selon l’information dont on dispose => 3 distributions: 
microcanonique, canoniques et grand canonique. 


-II- Distribution microcanonique 

1°) Définition de l’ensemble microcanonique 

En pratique un système physique ne peut pas être exactement isolé. 
a-Définition : Un système microcanonique est un système presque isolé, c’est à dire 
que les variations de son énergie E restent faibles et ne produisent aucun effet cumulatif 
d’échange d’énergie. On suppose alors que E appartienne à un petit intervalle: 
U <E <U + AU , avec AU « U . 


b- Opérateur densité microcanonique : 

L’espace des kets possibles (espace de Hilbert) est réduit au sous espace engendré par les 
vecteurs propres | c) de :H' de valeurs propres qui doivent être comprises entre U et U + AU . 


Si la dimension de ce sous espace est finie et égale à W , le principe d’équiprobabilité donne 


l’opérateur densité : 


P 


w i \ 

Vie) — (c | = — -este, si U<E<U+AU 


{ 0 . 


ailleurs. 


2°) Entropie microcanonique : 

a- Définition : E’ entropie microcanonique S m est définie par : S m -kLogW 


b- Propriétés : 

i) L’entropie S m est maxima par rapport à toutes les autres distributions (voir 
propriétés de S ( p) précédentes). 

ii) Pour un système presque isolé, l’entropie est maxima à l’équilibre thermique. 
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Remarque : Cette propriété peut être prise pour un principe : Le principe de désordre maximum. La 
deuxième propriété de 5 ( /?) précédente permis de déduire que toutes les complexions possibles sont 
équiprobables. Les principes de désordre maximum et d’équiprobabilité sont donc équivalents . 


III- Distribution canonique 


1°) Ensemble canonique : 

Supposons que le système Z, d’énergie E, interagisse en échangeant de l’énergie 
avec un autre système d? , d’énergie E ( , et ayant un nombre de degrés de liberté beaucoup 
plus grand. C 5 ' est appelé thermostat. I et ^ peuvent en particulier être de même nature. 
Z sera une partie d’un grand échantillon. E r » E : l’énergie est grosso-modo 

proportionnelle au nombre de degrés de liberté. ^ joue le rôle de réservoir d’énergie 
inépuisable pour Z => l’énergie E peut donc prendre n’importe quelle valeur. 

Soit le système global T.A>df comme étant isolé (si ce n’est pas le cas, on peut 
toujours élargir le thennostat pour qu’il en soit ainsi, en incluant s’il le faut tout l’univers). 
Les complexions possibles |c) de Z et |c c ^ de d? sont tous les états propres de leurs 

hamiltoniens LC et LC c , dans lesquels leur interaction mutuelle est négligeable, quelques 
soient les valeurs propres E et E c soumises à la condition de conservation de l’énergie : 

E totale = E + E c . 

2°) Loi de probabilité : 

Les complexions possibles du système global sont : |c)|c c ^ . Elles sont équiprobables 
(selon le principe d’équiprobabilité). Cherchons P(c) d’une complexion de Z, état propre 
de LC de valeur propre E. P{c) est proportionnelle au nombre de complexions | «y ) du 
thermostat qu’on peut associer à |c) pour former une complexion |c)|c c ^>du système global 

d’énergie totale : E totale = E+ E c . Si W c ( E c ) est le nombre de complexions de d? pour E c , 
alors : P(c ) oc W c (E c = E lotale - E ) . Considérons que le thermostat est presque isolé, on a : 

S c {E c ) = kLogW c (E c )=kLogW c {E c =E totale -E) et P(c) oc E) . L’hypothèse: 


E c » E , c’est à dire: E «c E lolale conduit à: P(c) oc e pE 
probabilité canonique : P(c) = 


„-P E 


AP) 


OU 


APh 2 > 


-PE 


1 dS c , . , 

avec /?= — . La loi de 

k dE c 

est la fonction de partition. 


3°) Opérateur densité : 

Sur la base ||c}} de vecteurs propres de ‘K , l’opérateur densité est diagonal; il s’écrit : 


P = 


AP) 


-pH 


avec 


Z(/?)= Tr(e- px ) 


Remarque : Les valeurs propres de LC peuvent être dégénérées. Dans ce cas, il y a plusieurs complex- 
ions pour une même énergie. Celles-ci sont équiprobables. 
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4°) Distribution en énergie 

Lorsque Z a lui même un très grand nombre de degrés de liberté (i.e de constituants) ; 


la densité de ses niveaux d’énergie: 




dW 


dE 


est très grande. 


La probabilité d’observer une énergie E peut être considérée comme continue. La 
densité de probabilité est: P(E)dE=d?(E)P(c)dE où P(c) est la probabilité d’une 
complexion d’énergie E. La densité P{E) présente un maximum très aigu, avec une 

dispersion très faible, à l’énergie la plus probable qui diffère très peu de E (voir marche au 
hasard). L’énergie la plus probable est obtenue en maximisant : 

1 Ô 

P(E)dE cc exp-[5 c (£ toMfe - £) + 5(£)], c.a.d, en résolvant : — [s c (E totale - E) + S( £)] =0 


soit : 


dS dS r 


8E dE r 


.S et 5 C : entropies de Z et Y considérés comme s’ils étaient 


En - E 


isolés. Cette dernière égalité exprime l’équilibre entre Z et Y que l’on peut encore exprimer 

1 dS 


par: p c où p z =-^~ et 

k dE 

5°) Valeurs moyennes : 

E =YEP(c)=—YEe- pE = <z> 

Y z Y zdp 


B = 

° k ÔE r 


En - E 


- = _ dLogZ(fi) 


dp 


L’écart quadratique moyen (écart type) entre l’énergie et sa moyenne : 


ct 2 (£)= E 2 -E 2 = — 


' 1 dZ^ 

<=> 


{zdpj 


v ' dp 


La relation liant la chaleur spécifique C v aux fluctuations de l’énergie cr (E) : 

1 dE 


1 1 

C v = — , or B = — => dp = -dT 

dT kT kT 2 


donc : 


C v =~ 


kT 2 dp 


Par conséquent : a 2 (£)= E 2 - E = kT 2 C v . 

• Entropie en distribution canonique: S c = -kTr(pLogp) avec: 

S c = - kTr\p(LogZ + fljf)] . Soit : 


P 


1 


-y 6H 


S c = kLogZ + k PU 


AP)_ 

où : U = E . 


Propriété : Parmi toutes les distributions, correspondant à la même énergie moyenne U 
(ou à une énergie moyenne inférieure), la distribution canonique est celle qui a la plus grande 
entropie. 

Remarques : 

i) Cette propriété généralise la propriété analogue pour la distribution 

micro canonique : Le thermostat fixe l’énergie moyenne U de Z (à travers (3). Compte tenu de cette 

contrainte, la distribution qui s’établit à l’équilibre thermique est celle d’entropie maximum. 

ii) Réciproquement, le principe de désordre maximum à l’équilibre, compte tenu des 
contraintes, permet de retrouver la distribution canonique. 

iii) L’entropie S (p) peut être considérée aussi comme fonction S ( U ) de l’énergie moyenne. 

q2 £ 

La concavité de l’entropie s(p) implique : — T < 0 . 

8U~ 
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6°) Ensemble canonique généralisé : 

a- Défïnition : Le système Z échange maintenant avec son complémentaire d? de 
l’énergie et du volume. Les conditions de contraintes s’expriment dans ce cas par: 

Z ^ -> Zu<^ 



K Jf c H toi 

E E c Etot — E + Ec 
V Vc V,o, -V + Vc 


b- La loi de probabilité : 

Suivant la même démarche que pour un ensemble canonique, on a : 


P(c) = 


w. 


\E totale E,V totale V J -S m {E,o,ale- E , V totale ~ V j 


W 


[e totale , b totale j 


oc e 


. Or : H ->■ E = E + pV , donc : 


E totale - E = E male - E , V totale ~ V = V mal ~ V => S, n [E totale ~ E , V totale 


V )= ( E totale -E, V iolal - V ) 


S,„ [Etotale - E , Y totale ~ V)= S„, (E^ , E ^ j ~ b (' ^ } = S,„ (E,^ , V^) ~ ~ 


avec : 


0 1 1 8S m 0 1 

j,n kT k ÔE et j "' Pm k 


dV 


r dS ' 


\ dV E 

'■ ' ’ totale t E, 


totale t ^ totale 


A la limite thennodynamique (que nous verrons dans la suite), il y a équivalence entre 
grandeurs micocanoniques et grandeurs canoniques. Donc : 

( E lotale — E , V totaI — V ) = 5 ( E totale — E , V total — V ) = 5 ( E totale , V totale ) — — — . 


Soit alors : P[c) = 


1 


, totale . V male) -/j[£(V) + pV~\ 


W (E totale , V ,otale ) 

Donc, la loi de distribution canonique généralisée est : 


p( c ) 


Z (T,p) 


-P[e(v) + P v] 


ou 

Z(T,p) = jdVY j e- fl[EiV) + pV] 


Or: Jdb^>(c) = l^ 

0 c 

définissant la fonction de partition canonique généralisée. 

b- Opérateur densité : p = 




Z{T,p) 


ou : 


eu 

Z = \dV Tr{e~ p *) 


IV- Distribution grand-canonique 


1°) Ensemble grand-canonique : 

Supposons que Z et ^ interagissent en échangeant de l’énergie et de la matière : c’est 
à dire que E et n (nombre de constituants) de Z ne sont pas fixes. Seules E et n sont 


constantes. Mêmes hypothèses sur ZuA 


[ E tota i e ~ E + Eq 


{ n ,o,ale =n+ n c 


E « Er 


n « n r 
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n étant variable, l’espace de Hilbert des états de Z est obtenu en faisant la somme directe 
©des espaces E pour un nombre n donné de constituants : E ©E ®...®E ©...= E. 

E est l’espace de Fock. Les opérateurs agissant sur cet espace sont : H , p et aussi 
l’opérateur N : opérateur nombre de particules, dont les états propres sont les états à n 
particules (n e N ) tels que : [iV, # ] = 0 . 


2°) Loi de probabilité : 

P(c ) où |c) est une complexion d’énergie E et de n particules est donnée par le 
principe d’équiprobabilité. Donc : P(c) ce W c ( E totale - E , n totale - n). 

L’hypothèse : E<s.E tolaIe et n «c n tolal conduit à : 


P(c) 


oc e 


an- pE 


o 1 dS C 

avec B — 

k 8E r 


et a--- 


\_dSç 
k drir- 


La loi de probabilité grand-canonique est par conséquent: 

est la grande fonction de partition. 


f(c) = 


Z 0 («./î) 


an- PE 


, ou : 


Z G (a,p)^ e 


an - PE 


3°) Opérateur densité : Sur la base ||c)} de vecteurs propres de H , la matrice densité est 

avec : Z c (a ,p)=Tr[e aN ~ pK ) ) 


diagonale : 


P = 


aN - PH 


Z c [a, p) 


<=> 




où H n est le hamiltonien partiel pour n particules et Z cn est la 


fonction de partition canonique correspondant à n particules. Z c {a , /?) est la transformée 
de Laplace de Z cn [p ) . 


4°) Valeurs moyennes : 


8LogZ c (a,p ) 

et 

dLogZ c (a,P ) 
n = 

dp 


8a 


• Le rapport : a/ P - p (~ E, est l’énergie de Lermi pour les métaux et semi conducteurs). 

• Entropie en distribution grand-canonique : — S (e , n ) = Log Z G [a , P) + PE - an . 

k 

Il suffit pour démontrer cette expression de considérer S (p) = -kTr[pLog p) avec: 
1 


Z G ( a , 


i ôs(e , n) 

» aN - pw . On en déduit: a=-~— - 


ôn 


i dS(E,n) 

et p= 

k dE 


Remarque : Ces relations montrent que l’équilibre thermique étudié par la distribution 

grand- canonique peut être décrit par (i)soit la grande fonction de partition Z G (a, fi) , soit l’entropie 

s(E,n) . 


V- Distribution d’équilibre en mécanique statistique classique 

1°) Ensemble microcanonique ou ensemble (N, V, E) 
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a- Définition : 

L’ensemble microcanonique est défini par une densité de probabilité uniforme dans le 
volume de l’espace des phases Q limité par les surfaces E et E + SE : 


Pm 


c si E<!h~(q, p) <E + SE 
0 ailleurs 


Pour le cas simple de système à N particules identiques enfennées dans un volume V, le 
triplet (N, V, E ) est appelé: paramètres fixes ou variables indépendantes de l’ensemble 
microcanonique. Ce sont celles-ci que l’on peut contrôler de l’extérieur. 


b- La constante c : J p m (q , p) dr est un nombre sans dimensions (nombre de 


microétats dans Q.). dT=Y\d i r t d" p, est un élément de volume de Q tel que : j dY = 

1=1 

nombre de points représentatifs du système à N particules (ou constituants) = nombre de 
microétats classiques à N particules = N\h 3N , où lé est le volume d’une cellule de Q. 
contenant un point et N ! est le nombre de pennutations de ces N particules. Par conséquent : 

[ p m dT = cft/r =1 <=> c = — r . 

J J N[h ™ 


Ainsi la constante c est définie comme étant le nombre de microétats par unité de volume 
dans l’espace des phases Q. . 


Remarque : dr=üd 3 r i d 3 


M= 


Kg.m 


=(js) =>[c] =(/.y) =(action ) 


c- La fonction de partition microcanonique Z m : 



Z m ( E,N,V)= le nombre de microétats classiques. 

La notation symbolique ( E , SE) indique que le domaine d’intégration est la ‘’peau” : 
E<tK(q,p) <E + 5E . 

d- Le potentiel thermodynamique microcanonique S m : 

Théorème : A chaque fonction de partition, on associe un potentiel thermodynamique 
(appelé ainsi parce qu’on peut en déduire toutes les autres grandeurs physiques). La fonction 
de partition microcanonique Z m a pour potentiel thermodynamique l’entropie microcanonique 
(en fait c’est la negentropie - S m qui est le potentiel thermodynamique) telle que : 

S m (E,N,V)= kLogZ n ,(E,NfV) 
où k est la constante de Boltzmann. 


Propriété : Les dérivées de S m par rapport aux trois axes principaux définissent : 


i)la température inverse J3 m : 


ii)la pression p m : 


? 1 _\ ÔS m {E,N,V) 

m kT m k 8E 


P m {E,N,V)p m {E,N,V) 


\ ÔS m {E,N,V) 
k dV 
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iii)le potentiel chimique p m : p m ( E,N , V)ju m ( E,N , V ) ^ E,N ’ V ) . 

k oN 

2°) Ensemble canonique ou ensemble ( N ,V,T) : 

L’énergie E de l’ensemble canonique fluctue autour d’une valeur moyenne E dont la 
valeur est déterminée par la température T du thermostat, qui sera aussi la température du 
système =>Dans l’ensemble canonique il y a un nouveau paramètre indépendant T (sans 
indice) qui remplace E => Le triplet de variables indépendantes est donc : (N ,V,T) . 

Remarque : On montre que pour des systèmes macroscopiques T m coïncide avec T. 

a- Définition : L’ ensemble canonique est défini par : 


Pc(q,p,t) = -L e ~ mq ’ p ’ l) 


Z C (N,V,T) = - 


1 


f 


-PH(q.p,t) ,3 


N]h 3N J ~ d 'i d r N d Pl d Pn 

b- Le potentiel thermodynamique canonique est l’énergie libre (canonique) 
F c définie par: 


F c (N,V,T) = - kTLogZ c ( N,V,T ) 


Ses dérivées par rapport aux trois axes principaux définissent : 


i) l’entropie canonique: 

ii) la pression canonique : 

iii) le potentiel chimique : 


S C (N,V,T) = 

Pc {N ,V,T) L - 
p c (N,V,T) = + 


3F C (N,V,T ) 
ST 

dF c (N,V,T ) 
ÔV 

ÔF C (N ,V ,T) 
dN 


Un changement simultané de N ,V, T donne une variation de l’énergie libre de Helmholtz F c : 
ÔF ÔF 3 F 

dF r = — -dT + — —dV + — —dN = - S r dT - p r dV + u r dN qui est une relation connue en 

c ST dV dN 

thermodynamique. 

c- Ensemble canonique généralisé ou ensemble (N, p, T) 

Il correspond à une distribution canonique généralisée ou distribution e ’T - p”. 
Supposons que le système Z échange avec son complémentaire C de l’énergie et du volume : 
C est donc son réservoir qui lui impose sa température T (thennostat) et sa pression p, 
constituant par conséquent les deux variables extérieures imposées à Z ■ 

L Définition : L’ ensemble canonique généralisé est défini par : 


1 ~/3H(q,p,t) 


p c (q,p,t) = ~-e 
Zc 


~ 00 1 


e rf 3 r x d 3 r N d 3 p { d 3 p N 


df(q,p,t) ~^>E = E+pV : enthalpie (valeur propre de tK(q,p,t)). 


ii- Le potentiel thermodynamique canonique généralisé est l’ enthalpie libre 
G(N,T,p ) définie par: 


G(N,T,p) = -kTLogZ(N,T,p) 
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Ses dérivées par rapport aux trois axes principaux définissent : 

~ . , dG(N,T,p) ~ , 

S c (N,T,p) = v ^ ’ O S c -kLogZ + pk(E + pVj 


• l’entropie : 

• le volume: 


V (N , T , p) = + 


dT 

dG(N,V,T) 


1 dLogZ 


le potentiel chimique : p(N ,T,p) - + 

iii- Valeur moyenne de l’enthalpie : 


dp p dp 

dG(N,V,T) _ 1 dLogZ 


dN 

E + pV= - 


p dN 

dLogZ [N, T, p) 
dp 


3°) Ensemble grand-canonique ou ensemble (p, V,T) 

a- Définition : L’ ensemble grand- canonique est défini par : 


Z G (p,V,T) = £■ 


1 




f - 0h(n ,(q , pf ,t) 3 - j 3 ' i3 * ,3 1 V 

} e d r } d r N d p x d p N =2_,< 


fpTo N ! h 

où Z C n est la fonction de partition canonique pour N particules. 


pfoN y 

C Zv. 


CN 


b- Le potentiel grand-canonique est : D. c (p,V,T) - - kTLogZ G (p,V,T) 
Ses dérivées par rapport aux axes principaux : 
l’entropie grand canonique : 


î- 


ii- la pression grand canonique: 

iii- le nombre de particules moyenne : 


s G = - 

Pc (N,V,T) = - 


dn(p,V,T) 


dT 

dn(p,V,T) 


N G (p,V,T)= - 


dV 

dn(ju,V,T) 


dp 


c- Dispersion du nombre de particules : 


N 1 

où n = — est la densité de particules et k t = — 
V V 


r dV ^ 


\ d P J T 


N 2 - N = NkTnK r 


est la compressibilité isotherme. 


4°) Equivalence des ensembles d’équilibre 
a- Limite thermodynamique : 

On étudie les systèmes tels que N est très grands. On pourrait penser que plus N est 
grand plus cette étude est difficile ! Paradoxalement, le cas limite N= oo redevient un cas 
relativement simple, et cela grâce aux méthodes statistiques. Ainsi sert-on en physique 
statistique de concept de la ‘’limite thermodynamique” : 

A la limite d’un système infini, les grandeurs intensives restent fixes. Autrement, V — >co on a : 
N E 

i) — et — restent fixes => N , £— »o o pour un ensemble microcanonique, 

N 

ii) — et T restent fixes => N — >cc pour un ensemble canonique, 

iii) p et T restent fixes pour un ensemble grand-canonique 

Représentation suggestives de la limite thermodynamique dans un ensemble microcanonique : 
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T T T 


0.125 

1 

8 

V (en cm 3 ) 

N , n22 E 

)=> — = 10" ; — = 4 fixes 
V V 

10 22 

8. 10 22 

64. 10 22 

N 

0.5 

4 

32 

E(en J) 



b- Extensivité des potentiels themiodynamiques : 

Définition : Une grandeur physique est dite extensive si elle est proportionnelle à la 
taille (le nombre de particules N ou le volume V) du système tout au moins dans la limite d’un 
très grand système. Autrement dit : A est extensive si 3 a constante telle que pour N — » oo on 
peut écrire : 

A — aN + [termes qui croissent moins que N ] — a N ou encore : 

— =a + [termes qui tendent vers 0 quand N — » go ] ou encore : 

Propriété : Les potentiels thermodynamiques sont extensifs. 

Exemples : - S ,E , F sont extensives, alors que - p , p , T sont intensifs. 


lim — -a . 

n — ^co yy 


Rq : Si A dépend de plus qu’une seule grandeur extensive, alors à la limite thermodynamique, 
a ne dépend que des rapports de ces grandeurs. 

Exemple: L’entropie microcanonique S m (A,E,£')peut s’écrire: S m (N,V,E) -Ns(n,s) où 
N E 

n = — ; £ ~~ : énergie par particule et s est l’entropie par particule. 

De même : F c [N , V , T) = Nf (n , T) où f est l’énergie libre par particule. 


c- Equivalence des ensembles : 

La densité pest très différente dans Q selon que l’on considère un ensemble 
microcanonique (m), canonique (c) ou grand canonique(G-c). Néanmoins, les trois ensembles 
sont équivalents, sous conditions mentionnées ci-dessous, 
i) Equivalence (m <=> ç) : 


Paramètres fixes 
(‘’de contrôle”): 


Grandeurs déduites : 


m 


N,V, E 


T m 

Mm 


> fonctions de N , V, E 


c 


N,V,T 

S c 


M c 


- fonctions de N , V, T 


Dans l’ensemble (c) , on peut ajuster la valeur de E c en faisant varier T. 
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Supposons que l’on choisit T telle que E c = E , où E est la valeur prescrite d’un ensemble 
microcanonique. Alors on aura : T m = T , S m = S , p m = 
ii) Equivalence fc <=> G-c) : 

c G-c 


Paramètres fixes(‘’de contrôle”) : 


Grandeurs déduites : 


N,V,T 

Y 


M,V,T 


S c 

Mc 


fonctions de N ,V, T 


AC 


fonctions de /j,V,T 


Dans l’ensemble (G) , on peut ajuster la valeur de N G en faisant varier p . 

Supposons que l’on choisit p telle que N G = N , où A est la valeur prescrite d’un ensemble 
canonique. Alors on aura automatiquement (dans la limite thennodynamique) : 
E c — E g , S C -S G , p c — //,... 

Une fois cette équivalence des ensembles établie, on peut calculer une grandeur physique 
dans l’ensemble de son choix. En général les calculs sont plus faciles dans les ensembles (c ) 
et (G) que dans l’ensemble (m). 


5°) Théorème de l’équipartition de l’énergie classique 

a- Enoncé Si l’énergie est un polynôme homogène du deuxième degré par rapport p 
variables, chacune correspond à un degré de liberté, à l’équilibre thermique l’énergie 

1 -fl ' 

moyenne est —k B T par degré de liberté. L’énergie moyenne totale est : E=p —k B T 

2 12 


tv Conséquences : 

E est répartie uniformément sur les degrés de liberté et elle est proportionnelle à la 

QË 

température T. Par conséquent la chaleur spécifique C - — est une constante. 


Remarque : Ce théorème n’est valable qu’à la limite classique de la mécanique quantique. 


Démonstration du théorème : Soit le polynôme £(x, , , x p ) diagonalisé par une 

transformation linéaire : x t = ^ a tj y ( où a tJ éléments d’une matrice A diagonalisant à la fois 

j 

l’énergie cinétique T et potentielle Y. Alors, le polynôme £(x, , , x p j devient: 


E = ^ b t y ‘ . La fonction de partition canonique : 


Zc = 


N\lr 


WM 


-PE 


dy x dy 


1 , P . C 2 ^ 

Z c se réduit à un produit d’intégrales : Z c = ^ ]~[ J e~ pb,y ‘ dy t . Or : J dxe 


Alors : LogZ c — Logfl + este . D’où : 


- = _ dLogZ c 


dp 


P 

2 P 


c.q.d. 


exosup.com 


page facebook 


34 


Chapitre 4 : 

REVERSI Bl Ll TE Ml CROSCOPI QUE ET I RREVERSI Bl Ll TE 

MACROSCOPIQUE 

I- Principe zéro de la thermodynamique 

1°) Enoncé : "Deux systèmes, en équilibre thermique avec un troisième sont aussi en 

équilibre thermique si on les met en contact”. 

2°) Equilibre thermique de deux systèmes 

• I, et E, deux systèmes mis en contact, c’est à dire, couplés à l’aide d’une interaction 
faible produisent un échange d’énergie. 

• Le système globale I, u E 2 est isolé. Le macroétat de I, u E 2 est stationnaire. Il en est de 
même pour le macroétat de toute partie de I; ul 2 : en particulier Ij et I 2 ; ils sont en 
équilibre thermique. 

• La condition d’équilibre de E, et E 2 est équivalente à la condition d’équilibre de S, u I, 
isolés, à savoir que son entropie S to t est maximum. 

En général : S tot < S, +S 2 (sous- additivité) 

Mais l’interaction étant faible => S , ( E, ) et S 2 (E 2 ) sont alors leur entropie, qu’ils auraient 
s’ils étaient isolés, et : S tot - 5, (£))+ S 2 ( E 2 ). 

La condition d’équilibre devient : S, ( E, ) + S 2 ( E 2 ) est maximum V(E,, E 2 ) avec la 

. i • ô'5, dS 2 

contrainte suivante: E, + E 2 =Cste => la solution : — - = — — . Soit p, = p 2 . 

dE x ôE 2 

Remarques: 

i) D’une part, l’interaction doit être suffisamment faible pour que E, et I, puissent être considérés comme 
isolés, d’autre part, il ne faut pas qu’elle soit nulle pour qu’il y ait échange d’énergie. 

ii) La condition d’équilibre : [\ = fi 2 détermine les énergies Ej et E 2 connaissant l’énergie totale. L’état 
macroscopique d’équilibre de chaque système est alors déterminé. 


3°)Genéralisation à deux systèmes échangeant de la matière 


La condition d’équilibre est : S i (E l ,N l ) + S 2 (E 2 , N 2 ) est maximum avec les deux 


contraintes : £, + E 2 -Cste et + N 2 -cte . => < 
températures et des potentiels chimiques. 

4°) Démonstration du principe Zéro : 

S, en équilibre avec Z 3 : /?, = 
Etat initial : < E, en équilibre avec I 3 : [E, - 
S, et Z, isolés entre eux. 




: Égalité des 
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S, et I 2 interagissent. Puisque initialement J3 { = J3 2 = J3 3 , 

Etat final : < ils s'atisfont déjà la condition d'équilibre; ils sont en équilibre thermique 

entre eux et leur interaction ne modifie pas cet état initial. 

5°) Echelle des températures : 

• Le paramètre [1 définit une échelle de température telle que deux systèmes en équilibre 
thermique soient à la même température 

• Toute fonction monotone de j3 définit aussi une échelle de température. L’échelle des 

_j f ôs y 1 

températures absolues correspond à la fonction : T = (k/3) = — 

\dE ) 

• Dans la pratique une échelle de température est obtenue par la mesure d’une grandeur 
macroscopique, dite grandeur thermométrique fonction de P et donc de T. 

Exemple; On mesure la pression P d’un volume V = cste d’hydrogène gazeux, et on définit 

PV 1 

l’échelle de température ordinaire: t = . R et a sont déminés par deux point fixes; 

R a 

Le point triple ; t = 0,01 °C et le point d’ébullition t = 100°C de l’eau à la préssion 

atmosphérique. L’hydrogène et l’enceinte qui le contient est un thermomètre. 

La température d’un système est obtenue en le mettant en contact avec un thermomètre. 

Celui-ci doit être assez petit pour ne pas perturber notablement le système. 

Remarque : La température ordinaire est une grandeur pour laquelle est simplement 

définie l’égalité ; tçz l )=t(ï . 2 ) , si s, et z 2 en équilibre thermique, mais le rapport des deux 
températures n’est pas défini. La température ordinaire n’est pas mesurables, elle est 
seulement repérable. Par contre le rapport des températures absolues est bien défini. La 
température absolue est grandeur mesurable. 

II- Premier principe 

Le premier principe exprime essentiellement la conservation de l’énergie au travers 

des concepts de travail et de chaleur. 

1°) Energie interne: U=E = Tr pH 

Conservation de l’énergie ; 

a- Système presqu’ isolé (système microcanonique) ; U = cte 




= -^-7>(p£ . Soit: ^- = -^-Tr^HpH- pdOf) = 0. 
b- Systèmes fermés : I et I' échangent de l’énergie mais pas de matière : 


= 0 => — = Tr 
dt dt 


d p 
dt 


A U j. + A U z . = 0 . 


2°) Définition du travail et de la chaleur 

a- Transformation adiabatique lente 

Soit un système X et TC (x , , x 2 ,...., jc. ,....) son hamiltonien où les paramètres x t varient 


lentement avec le temps. Alors : 


d‘K= V ÔJ{ -dx i = Y X.dx, . 
; ÔX; ' V 
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i) Définition: X i est appelé "force généralisée” associée à la variable x t . Le couple 
(Xj. , x t ) sont conjugués canoniques. X. est une observable de I . 


Exemples : 


x t 

* 

dtK 

V 

-p 

-PdV 

B 

-M 

-MdB 

y 

P 

pd(f> o 


ii) Théorème adiabatique : 

L’évolution de I est générée par une perturbation IV (t ) variant lentement dans le 


temps telle que son hamiltonien devient : 


dC (t)= tK 0 + W (t ) où dC 0 -cîe 


avec: { 


à t - 0 


w( f ) 


x 0 =x(o) 


| £ o) 


t ^ 0 

«■(») 

E(I) 

|E(0) 


Enoncé du théorème adiabatique : ‘’La probabilité de transition d’un état propre de J{ 0 
vers un autre état propre de .7/„ est négligeable pour une perturbation lentement variable 
dans le temps ’ ’. 

A t = 0 , p est diagonal sur la base j £j (0)^| de , car le système est alors en équilibre 
thermique : p(0) = (°)K (^ (°)| = 2>aAu(°) avec : p Â =p Xi ( 0). 

À A 

Le théorème adiabatique les p Â sont constantes. 

A t quelconque, on a : p(t) - Y\ E x (0 )Px { e a (Oh Z PxPxx (0 • 

A A 

La distribution de probabilité p À sur les états propres de l’énergie reste la même. La matrice 
densité constante sur la base de ces états propres. Elle est toujours diagonale : 

p^=( e M e ,)=t.P‘{ e »\ e *}{ e a e a 

A 

= Z /h W = =P Mfl ■ diagonale 

A A 

I reste en équilibre thermique. 


Remarque : La matrice densité p reste invariable; l’entropie S, c’est à dire, l’information manquante reste 
constante. 

iii) Définition du travail : 

Dans une transformation adiabatique lente, infinitésimal, le travail dW échangé par le 
système avec son milieu extérieur (ME) est : 


dW -Tr(pddT) . 


En effet, la variation de l’énergie interne dU au cours d’une telle transformation est réduite au 
travail, soit: 
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U = Tr (pX ) - Z ftAi = Z Pa < I X I ^ >= Z C'A => 

ûW = dV = ^p À dE À = ^p À <Ei\dX\Ei>-Tr(pdX) 

x i 

Explicitement: dW = 'Y J Tr(pX j dx i ) = ^X i dx i avec X, = Tr[pX.) 

i i 


Le travail reçu est le produit de la force moyenne X t par de le déplacement dx t . Ce travail est 
associé au déplacement des niveaux d’énergie. 

Force moyenne X associée au paramètre x en distribution canonique: 


X = 


dX 

ôx 


et p = 




X = Tr(pX) = -Tr(—e~ pH )=- — Tre . 
Z dx pZdx 


en utilisant : dTr(f(A )) = 7>( — dA ) où f est une fonction de l’opérateur A, on a : 

dA 


X =- 


1 dLogZ 
P dx 


: Equation d’état. 


Exemples: Variation de volume : 

— > 

Variation de B : 


dX _ — _ 1 dLogZ 

dV~ p dV 

— =-M 1 dLogZ 

dB p dB 


iv) Définition de la chaleur : Soit E isolé est initialement en équilibre thermique. A t on 
obtient un couplage faible avec un autre système par une interaction V (t ) qui est une 

perturbation fluctuant maintenant rapidement au cours du temps, mais restant petite. La 
propriété énoncée pour une variation lente n’est plus valable => Existence de transitions entre 
différents niveaux d’énergie. Les probabilités p À des états propres de X(t)-X 0 + V {t) 
sont modifiées. V (t ) produit une variation dp de l’opérateur densité. La variation de 
1 ‘énergie interne est : dU = ^ E À dp À = Tr ( X 0 d p ) 

/l 


Défi La variation de l’énergie interne, dans une transformation infinitésimal générée par une 
perturbation V (t ) fluctuant rapidement au cours du temps, tout en restant petite, est égale à 


la chaleur dQ échangée par le système avec son milieu extérieur telle que : 


dQ-Tr[ X 0 d p) . 


dQ correspond à un déplacement de la distribution de probabilité des états | E À > . La chaleur 
reçu est positive quand l’énergie interne augmente. Les états de plus grande énergie 
deviennent plus probables. 

Remarques : (i) L’évolution de p n’est plus donnée par : i h —=\ jt n , p] , puisque cette 

dt 


évolution est due à la perturbation v (t ) . 

(ii) Transformation quelconque : La variation de l’énergie interne accompagnant une telle 
transformation est : du = Tr [d ( /?#)] = Tr ( pdor ) + Tr [ptd p) = =dw + dQ. 


5°)Enoncés équivalents du 1 er principe 

a) ‘’Dans une transformation faisant passe d’une état d’équilibre à un autre, la variation de 
l’énergie interne est égale a la chaleur plus le travail reçu : A U = A W + AQ ’ 
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b) ‘ ’La sommes dW + dQ ne dépend que de l’état final et de état initial. Elle ne dépend pas 

du chemin suivi : Us - U, = À W + AQ ’ 

c) ‘’Dans une transformation adiabatique (dQ=0) la variation de l’énergie interne est égale 

au travail reçu : AU = AIE . 


6°)Transformation à pression constante 

dQ = dU + PdV = dH , avec H =U + PV est l’enthalpie du système. Chaleurs spécifiques à 

-dW 


volume constante : C v 


f dQ ] 


f St/A 

[dT ) 

V 

Uf J 


et à pression constante : 


(dQ^ 

- 

{ dH \ 

l dT ) 

r 

Ut J 


7°) Application à un gaz parfait classique 

3 5 3 5 

U = -RT ; R = Nk ; PV = RT ; H = —RT alors : C v =-R ; C P =-R 
2 2 2 2 


=> C P - C v -R : Relation de Mayer. 
Transformation adiabatique lente : 


jdU =C v dT = -PdV 
j PdV + VdP = RdT 


<=>< 


dT = ^[PdV+VdP] 

CV 

PdV + —[ pdv + VdP ) = 0 


« 


dP C P 
— + — 
P Cy 



£jl 

Cy 


y « PV y = Cte 


III- Deuxième principe 


Le second principe associe l’irréversibilité des transformations macroscopiques à une 
perte d’information c’est à dire à une augmentation d’entropie. 

l°)Transformation quasi-stationnaire : Soit une transformation : A — » B 

état équilibré état équilibré 

Les états intermédiaires entre A et B ne sont pas des états d’équilibre. Sinon, il n’y aurait pas 
d’évolution de A vers B , puisque un état d’équilibre (intermédiaire) est un macroétat 
stationnaire qui se maintient indéfiniment. Or, le formalisme de la physique statistique 
présenté ici ne décrit que des états d’équilibre thermique. On considérera donc seulement des 
transformations dites quasi stationnaires. 

Déf. : Les transformations quasi stationnaires sont des transformations où les états 
intermédiaires sont suffisamment proches des états d’ équilibre. 

Propriété : 

* Pour qu’une transformation soit quasi stationnaire, il faut que l’évolution du système 
soit plus lente que le temps de relaxation des états intermédiaires vers les états d’équilibre 
thermique voisins. 

* Il en est ainsi si la résultante des forces extérieures est très petite et si les gradients 

de température dans le système et son milieu extérieur sont très faibles. 

Exemple : La compression d’un gaz est très lente si la force extérieure exercée par le piston 
excède très peu la pression cinétique du gaz sur le piston. Pour tous déplacement du piston, 
même petit, le gaz évolue de manière presque instantanée vers l’équilibre thermique. 
Pratiquement tout au long du déplacement du piston le gaz reste en équilibre thermique bien 
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qu’il évolue. De même un très faible gradient VT implique des températures presque partout 
égales. Pour tout transport de chaleur produit par VT , le système s’adapte presque 
instantanément à un nouvel état d’équilibre thermique. 

— > _ 

Rqs : (i) Les forces de frottement : F = -kv qui détruisent l’équilibre thermique s’annulent 
— > —> — » 

avec v et la perturbation de l’équilibre produite par F disparaît avec v . 

(ii) A la limite où la vitesse d’évolution du système est nulle; le sens de cette évolution peut 
être inversé : il suffit de modifier très peu les forces pour changer le signe de la résultante ou 
d’une petite variation de T pour changer le signe de VT . Une transformation quasi- 
stationnaire est réversible. 

— ^ 

Inversement, une transformation réversible à tout moment => une résultante de F et 

— ^ 

V T presque nuis. Une transformation réversible est toujours quasi-stationnaire. 

Une transformation adiabatique (exemple: la détente adiabatique d’une gaz ou encore la 
désaimantation adiabatique d’un seul paramagnétique) est quasistatique. 

Deux systèmes initialement isolés et en équilibre thermique, mais avec des températures 

— ^ 

différentes, leur mise en contact thennique => V T important avec transfert de chaleur de 
l’un vers l’autre : ils évoluent très rapidement. La transformation qu’ils subissent n’est pas 
quasi-stationnaire ; elle n’est pas réversibles. 


2°) Variation d’entropie statistique (dans une transformation réversible). 

L’entropie statistique est définie par : 5 = -kTr(pLogp) . 

Sa variation est : dS - -kTr\J^pLog p)' d p^~ -kTr\j^Log p + \)d p^ = -k,Tr(dp Logp) 

pour une variation dp de l’opération p . 

Chaque état intermédiaire d’une transformation quasi-stationnaire est un état d’équilibre 
thermique, alors : p = ^ en supposant que seule l’énergie est échangée avec le 

milieu extérieur (ME). D’où: dS = -kTr'ÿd pi^-LogZ - /Uf)] = kTr[dpLogZ + fid ptt]. 


Donc : 


dS — flkTr^Jfd /?) = k j3dQ . 


(S : entropie statistique) 


Rqs : ( i) En thermodynamique l’entropie et la température absolue sont définies 
conjointement par la relation : dQ = TdS , pour une transformation quasi-stationnaire. 


L’entropie statistique est donc identifiable à l’entropie thermodynamique. La température 
absolue est alors définie par : T - l/k/3 (k = Constante arbitraire) . 

Dans une transformation quasi-stationnaire, la variation d’entropie, c’est à dire, l’information 
manquante sur le système, est contrôlée; elle est proportionnelle à la chaleur reçue. 

(ii) En général dQ n’est pas une différentielle totale exacte, c’est à dire, n’est 

pas la différentielle d’une fonction de Z (fonction de partition). Mais ds l’est. Le facteur t/r 
est un facteur intégrant pour la différentielle dQ . 

Propriété : dans une transformation adiabatique lente, la variation d’entropie est nulle. 


3°) Réversibilité microscopique 

Les phénomènes à l’échelle microscopique sont invariants par renversement du sens 
du temps . Ceci veut dire que : 
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A la suite des états: 

par lesquels passe le système aux instants: 
correspond une suite inversée : 


KO; KO; 


h-i ^ h ^ h+i 

K _ 0; K _ *.--i); 


aux instants: -t i+l < -t i < - t i t 

et que cette suite inversée est aussi solution des équations d’évolution du système. Cette 
propriété d’invariance résulte de la nature des forces qui régissent l’évolution des systèmes 
microscopiques. 

a) Point de vue quantique : 

Par suite de la nature complexe de l’équation de Schrôdinger, la suite innervée 
est définie par : y/ ' (t) = y/* (-t) 


Si i//(/) est solution de l’équilibre de Shrodinger : 


ih d ^- = mt)y/(t) eC 
dt 


ys'(t) satisfait à : 


ifi 


dt 




Condition pour une invariance par renversement du sens du temps : 

= ]. 

Remarque : La symétrie par renversement du temps est bien vérifiée à l’échelle 
microscopique, quelle que soit la nature de l’interaction entre particules. La seule 
exception connue est l’interaction dite superfaible, d’un type analogue à l’interaction 
faible responsable de la radioactivité / 3 . Cette interaction superfaible , non invariante par 
renversement du sens du temps , n’intervient que dans des phénomènes très rares tels que 
la désintégration : k° l -> . 

b) Point de vue classique : 

L’invariance par renversement du sens du temps à l’échelle microscopique 


apparaissant en M.Q, existe aussi en M.C, puis que celle-ci n’est qu’une approximation de la 

M.Q. L’intérêt d’étudier cette invariance en M.C est de déterminer quelles sont les forces à 

— ^ 

l’échelle microscopique qui possèdent cette propriété. Si q(t ) = {q x {t), ,<?„(/)} est 


solution de l’équation de Newton: 

— > — » 

inversée q\t) = q{-t) satisfait à : 


m - 


d 2 q, 

dt 2 


= Fi(q, q , t) 


d q , 

avec q = . La suite 

dt 


m- 


, d 2 q] 
dt 2 


= F t (q\ ~ q',~t ) 


d q' 
d(-t ) 


= -q 


Il y a invariance par renversement du temps si : 


F^q, q , t) =F i (q\ 



• Il en est ainsi pour les forces F qui ne dépendent ni du temps ni de la vitesse. 

• Les forces de frottement proportionnelles à la vitesse brisent cette invariance. C’est 
une première manifestation de l’irréversibilité microscopique. 


Remarque : La force de Lorentz F = qv a b exercée sur une particule de charge q et de 
vitesse v dans un champ magnétique b semble aussi briser cette invariance par 
renversement du temps. En fait, il n’en est rien, car le reversement du temps doit être 
appliqué nom seulement à cette particule, mais aussi au courant donnant naissance à b. Ce 
qui change le singe de T’ et de b , mais pas celui de f . 


exosup.com 


page facebook 


41 


4°)Irréversibilité macroscopique 

• La symétrie du reversement du sens du temps est brisée à l’échelle macroscopique. 

En effet, la quasi-totalité des systèmes macroscopiques ont une évolution qui n’est pas 
réversible. 

Exemples: la détente de joule ; la transformation du travail en chaleur ; la conduction de la 
chaleur... Cette irréversibilité apparaît, en particulière, dans l’évolution des être 
vivants => d’où la distinction entre le passé et le futur ! 

• Cet apparent paradoxe (réversibilité microscopique et irréversibilité macroscopique) 
tient à l’impossibilité de préparer un état pur quelconque d’un système macroscopique. On a 
toujours un mélange statistique incohérent d’états purs. 

En effet, si on savait préparer un état pur, on pourrait tout aussi bien en préparer un autre, 
ayant des vitesses opposées et donc l’évolution serait exactement l’inverse. 

Exemple : Détente de Joule : Toutes les molécules, sont à l’état initial dans le compartiment 
(1) , de volume V, , et en équilibre thermique. La case (2) , de volume V 2 , est vide. 



A t =0 , on enlève la paroi entre (1) et (2) . Les molécules se repartissent ensuite dans les 

deux cases. Au bout d’un certain temps, le gaz remplit tout le volume. Si il était possible à un 
instant t > 0 de changer le sens des vitesses de toutes les molécules, le mouvement inverse 

serait observé et toutes les molécules retournent dans (l) .Ceci est concevable si le nombre de 


molécules était petit, mais devient de plus en plus irréaliste quand ce nombre augmente. Au 
contraire le gaz tend vers un état d’équilibre dans lequel il occupe tout le volume offert. 


t <0 

t = 0 t= T » 0 

pTTl 

i 


o ° 

O 

O 

_ O 

* O 


si"- 


Détente inverse : (on inverse le signe des vitesses des molécules) 

t = -T <s: 0 t =0 



Etude de la détente de Joule dans l’espace des phases : 

L’équilibre thermique est représenté par une densité de probabilité uniforme dans un 
hypervolume Lj , sur l’hypersurface : E = cste . Après enlèvement de la paroi, Lj se déforme 
et se ramifie en projetant ses ramification dans tout l’hypervolume L 2 disponible de l’espace 
des phases. 
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Dans cette déformation le volume T, reste constant conformément au théorème de 

Poincaré (invariants canoniques). Cependant, l’hypervolume final ressemble grossièrement à 
une densité uniforme dans I\ (si on ne tient pas compte des détails) correspondant à l’état 
d’équilibre des molécules dans let 2. L’état final est donc pratiquement identique à cet état 
d’équilibre. 


Détente inverse : Considérons l’hypervolume r, obtenu à partir de r, par inversion du 
signe de toutes les vitesses. Tout point de r’ est un microétat conduisant à l’évolution inverse. 
T’ comme F’ est contenu dans f, . La probabilité d’observer une évolution inverse est donc : 


r',' _ vf _ 1 
r 2 - v 2 n ~ 2 n 


où V l est le volume de la case (l) et V 2 : volume de la case (1 + 2) . Cet 


événement est donc hautement improbable pour un système macroscopique (N > 10 23 ) . 


Conclusion : La détente de Joule conduit à un état final voisin d’un nouvel état d’équilibre 
et pratiquement indifférenciable de celui-ci. L’évolution inverse n ’ est pas impossible . elle est 
seulement très improbable : la détente de Joule est irréversible. 

Remarques : (i) Paradoxe de Zermelo : Un théorème de la mécamique analytique (théorème du retour 
de Poincaré) énonce que l’évolution de tout système mécanique est quasipériodique : le système revient 
aussi près qu’on veut de son état initial au bout d’un temps assez long. En fait ce théoerème ne conrtedit 
pas ce qui percéde (paradoxe de Zermelo) car le temps necessaire pour observer le retour est immense 
(plus grand que l’âge de l’univers) pour un système macroscopique. 

(ii) Le chaos moléculaire : Le remplacement de l’etat final complexe par l’état d’équilibre est 
équivalent à l’hypothèse de l’existence d’une interaction aléatoire qui déplace au hasard les points de 

l’hypreuvolume Fj dans l’hyprovolume T, . Cette interaction aléatoire hypothétique est incorporée 
dans le potentiel d’interaction entre molécules, qui décrit les collusions. C’est l’hypothese du chaos 
moléculaire . 


5°) Variation d’entropie d’un système isolé 
a- Point de vue classique 

La détenete de Joule constitue un exemple de système isolé en dehors de l’équilibre. 

Très précisément l’entropie S d’un tel système est constante. Mais pratiquement, on doit 
considérer qu’elle augmente ! !. 

Dans cet exemple : S 1 = kLogF v Or le théorème de Poincaré => l’état final 

a même entropie : S[ = kLogT\ = kLogF 1 . 
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Il y a conservation de l’information. Cependant l’état final n’etait pas différenciable en 
pratique de l’état d’équilibre, on doit considérer que l’entropie finale est plustôt : 

S 2 = kLogT 2 et alors S 2 >5^ 

En négligeant les détails qui différencient l’état final et l’état d’équilibre, on a perdu de 

l’information et S a augmenté. Les petites interactions avec l’extérieur en ont fait perdre aussi 
L’augmention apparente de l’entropie résulte de l’imprécision et de la grossièreté de nos 
expériences, du point de vue microscopique. Elle ne traduit pas une modification du microétat 
où se trouve le système. Autrement dit la densité final p[ est remplacée par p 2 de l’état 
d’équilibre. p 2 varie plus doucement dans l’espace des phases que la densité réelle p\ . Ce 
remplacement est une ‘’érosion " de la densité de probabilité. La densité "érodée " p 2 est 
appelée densité semi-finie ou densité grossière. 



Cette érosion peut être pratiquée en découpant l’espace des phases Q en petites 
cellules, dont les dimensions sont grandes à l’échelle microscopique, mais petites à l’échelle 
macroscopique. Dans chaque cellule la densité réelle est remplacée par une densité uniforme 
ayant la même intégrale dans le volume de la cellule. 

La valeur moyenne : A = J pdl de A(q, p) prend la même valeur avec la densité 

semi-finie qu’avec la densité réelle, pourvu que la grandeur A(q, p ) soit fonction douce de q et 
p dans Q . 


b- Point de vue quantique : 

En mécanique quantique aussi, l’entropie 5 =-kTr(pLogp ) est constante . 
L’évolution de l’opérateur densité donnant p(t) à partie de p(t = 0) est exprimée (en 


représentation d’interaction dite de Heisenberg) par : 
L’invariance de S par transformation unitaire de p 


p(t) - e p{ 0)e h 

=>S(f) = S(0). 


L’augmentation apparente de 5 peut être expliquée par le modèle suivant basé sur une 
hypothèse dite stochastique définie par: 

Le hamiltonien H du système est séparé en deux parties : ttf = :tf„ + V 

Jf 0 => évolution lente du système 
V => évolution rapide et complexe 

V Z) l’interaction avec l’extérieur et entre molécules. Elle est supposée être équivalente a une 
petite perturbation aléatoire. Elle est remplacée par un ensemble statistique de perturbations : 
{v; , p, } ; où Pi sont les probabilités associées à ces perturbations. 


L’évolution du système n’est plus prévisible. C’est un processus stochastique. 

• / \ ^ i / \ , . — [tf. + — 

L’opérateur densité à t est : p[t) = ^ p.p. (rj où p t (r) = e h p( 0) e h 


■K + w)]' 


c’est 


l’évolution de p correspondant à l’hamiltonien 


K 0 +V t et à l’entropie S (f>, ) = S (f>(0)) . 
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La concavité de S => > £(/?(())) . 

L’augmentation de S résulte de l’introduction d’un hasard dans l’évolution du système. 

Conclusion : Du point de vue macroscopique, ce système isolé évolue de telle sorte 
que son entropie augmente. Il tend vers un état d’équilibre thermique pour lequel 
l’entropie est maxima . 

6°) Enoncés du second principe : 

a)” Dans une transformation adiabatique d’un système macroscopique l’entropie ne peut 
diminuer. Elle reste constante si la transformation est réversible ; elle augmente si elle est 
irréversible ’ ’. 


b )“’ Dans une transformation menant d’une état d’équilibre 1 à une autre état d’équilibre 2, 
dans lequel un système macroscopique reçoit des quantités de chaleur A <2, d’une suite de 
thermostats, de températures respectives T , la variation d’entropie d’un système satisfait à 


l’inégalité de Clausius 


AS> Z 


A Qi 


La transformation de chaque thennostat est supposée quasi-satationaire. Les T. sont les 
températures des thermostats et non celles prises par le système au cours de la transformation 
(elles sont différentes en général ) 

L’égalité n’a lieu que si la transformation du système est réversible. Les 7^ sont alors celles 
que prend successivement le système. 

Une transformation réversible de l’état d’équilibre à f à l’état d’équilibre à T 2 implique un 

échange de chaleur avec une infinité de thermostats de températures T. tels que : f<T i <T 2 . 

2 2 


La variation d’entropie est alors : 5 2 - 5) = J = J C (r) où C(T ) peut être 


dT 


C v ou C P selon la transformation. 


c) Enoncé de Carathéodory : ‘ ’ Dans tout voisinage d’un état d’équilibre A quelconque 
d’un système macroscopique, il existent des états d’équilibre thermique adiabatiquement 
inaccessibles” . 



Les états inaccessibles sont ceux ayant une entropie plus petite que celle de A. 

Exemple : Pour un gaz parfait classique dans le diagramme P-V, la détente adiabatique 
réversible à partir de A est représentée par la courbe : PV r - Cste . Cette courbe est la frontière 
qui sépare la région accessible de la région inaccessible. 
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d) Enoncé de Clausius : ”11 est impossible de construire un système X, décrivant un cycle de 
transformations et n’ayant d’autre effet qu’un transfert de chaleur d’une source froide vers 
une source chaude.” 

Dans un cycle le système X revient à son état initial. 



Ti>T 2 


e) Enoncé de Kelvin Il est impossible de construire un système X qui fournirent du 
travail en ne recevant de la chaleur que d’une seule source. ” 

Cette forme du second principe implique l’impossibilité du ‘'mouvement perpétuel”. 



Remarque : Les énoncés b), c), d) et e) peuvent être déduits de a). De même a) peut être déduit de l’un 
quelconque des autres. 


7°) Troisième principe (Principe de Nernst) 

La probabilité p À d’occupation du niveau d’énergie E x est : 



e -fi(E x - Eq) 


„-P( E X - Eo) 


E 0 énergie du niveau fondamental. Quand T -> 0 , alors p x — » 0 pour A* 0 ( probabilité des 
états excités). Seule p 0 est non nulle ( p 0 : probabilité d’occupation de l’état fondamental). 
L’entropie au zéro absolu est: S 0 =kLogg où g= W est la dégénérescence de l’état 
fondamental. En général cette dégénérescence est très petite devant le nombre N de 
constituants du système. S 0 est très petite devant l’entropie S aux températures élevées, qui 
est proportionnelle à N. Pratiquement on peut poser : S 0 = 0 . 


Enoncé : ”A T = 0 K, l’entropie d’un système quelconque est nulle et minima, 

quelques soient les paramètres extérieurs [v , B,...j et quelque soit l’état physique ou 
chimique du système. ’ ’ 

Remarques : (i) On peut imaginer des systèmes particuliers qui ne satisferaient pas à ce principe. 

En pratique, aucun n’a été trouvé dans la nature. Par exemple, un ensemble paramagnétique de spins, 
Mans une induction magnétique nulle a une entropie : S=kLog 2 , indépendante de la température, et 
qui ne tend pas vers zéro quand T -» OK . En général un tel ensemble devient ferromagnétique ou 
antiferromagnétique, à basse température. Les spins s’alignent par suite de leur interaction mutuelle, et 
l’entropie s’annule. 

(ii) Le troisième principe est un phénomène quantique. 11 est dû à l’existence d’un état fondamental peu 
dégénéré. Par contre en mécanique classique, l’entropie n’est pas bornée inférieurement. 

Conséquences du 3 eme principe 

(i) Impossibilité d’atteindre le zéro absolu : 

Le zéro ne peut être approché que par des transformations adiabatiques (en absence de 
thermostats de température inférieures à celles qu’on sait réaliser!!). Conformément au 
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second principe(énoncé de Carathéodory) l’état à T -0K est adiabatiquement inaccessible, 
car c’est l’état d’entropie minima. 

(ii) A T-OK, toutes les chaleurs spécifiques s’annulent : 

Dans une transformation donnée, la variation d’entropie entre deux états d’équilibre 

thermique, de température 7j et T, respectivement, est : A S = [c — = S (T 2 )-S (T t ) où C est 

i T 

la chaleur spécifique de cette transformation. Quand T — > , l’intégrale donnant A 5, doit 

C 1 

rester fini, puisque S(0) = 0. Ceci n’est possible que si — croit moins vite que — ; c’est à 
dire C-»0. 

(iii) Aux très basses températures, les chaleurs latentes de chansement d’états 
décroissent plus vite que T , quand T — >0 K : 

La chaleur latente d’un changement d’états physique à T est la chaleur reçue au cours de 
cette transformation : L(T)-T AS . 

AS est la différence d’entropie entre l’état initial et l’état final. AS — r >0 — > 0 , donc L(T) est 
un infiniment petit d’ordre supérieur à un en T . 

Exemple : La chaleur de fusion de l’hélium He 4 a un point triple pour T= 1.74 K . 

A T >1.1 4 K : l’hélium solide fond en donnant He, liquide et la chaleur de fusion est 
proportionnelle à T . 

A T<\.14K : le solide fond en donnant He„ liquide superfluide. La chaleur de fusion décroît 
plus vite que T . 



8°) Potentiels thermodynamiques: 

a- Transformation adiabatique : dp = 0 (théorème adiabatique) => p uniforme 

=> ensemble microcanonique. 

Le potentiel thermodynamique est - S (appelé parfois négentropie). 

D’après le second principe, ce potentiel ne peut augmenter : - AS <0 . 

L’équilibre stable correspond à un minimum de ce potentiel, c’est à dire un maximum 
d’entropie. 


b- Transformation isotherme, isochore : 

Le potentiel thermodynamique est l’énergie libre 


F=-kTLogZ(x,j3) = U-TS , 


où x désigne tout paramètre extérieur (tel que le volume) et Z(x,(3) est la fonction de 
partition canonique. 

Evolution : dF < dW - SdT (l’égalité n’a heu que pour une transformation réversible). 
Dans ce type de transformation, on a dW -dT = 0=> F ne peut augmenter : dF < 0 . L’équilibre 
stable correspond au minimum de F . 
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c- Transformation isotherme, isobare : 

Le potentiel thermodynamique est l’enthalpie libre G(T, P) défini généralement par. 

G=-kTLogZ(x,j3)=H-TS . 

où H est l’enthalpie, x désigne tout paramètre extérieur (tel que la pression) et Z(x,/3) est la 

fonction de partition canonique généralisée. 

Evolution : dG< Vdp-SdT (l’égalité n’a lieu que pour une transformation réversible). 
Dans ce type de transformation, on a dp = dT = 0=>Gne peut augmenter : dG< 0. L’équilibre 
stable correspond au minimum de G . 


9°) Propriétés de l’énergie libre F 

a- F (x, T) est une fonction caractéristique : Ceci signifie que toutes les grandeurs 
thermodynamiques (valeurs moyennes ) peuvent être calculées à partir de F (x, T ) ; 

r q F \ 


Force généralisée : X = 


dx 


Jt 


U=- 


r F ' i JdF^ 

-F -T 


W 1 


dT 


équation d’état. - En 

Chaleurs spécifiques : 



f dF^ 

: S=- 


U T J 

C x - -T 

' d 2 F 

3'T' 2 


Energie interne : 


ô (/r) 

b- F [x, T) est une quantité extensive. 

c- Dans une transformation isothenne le travail fourni par un système au milieu extérieur 
est ; - dW < - dF . La diminution de l’énergie libre est le travail maximum que peut fournir 
le système; ceci justifie son appellation d’énergie libre. 

d- La probabilité d’une complexion est : P(C)= e p[t ~ £) . 


10°) Potentiels chimiques 

En distribution canonique N est un paramètre extérieur constant. Cependant, ce 
paramètre peut être varié comme le volume V ou le champ magnétique B . 

Cette variation peut être obtenue, par exemple, en introduisant ou en extrayant des particules 
au travers d’une paroi semi-perméable, ou bien encore par une réaction chimique qui crée ou 
détruit des particules. Cette situation est différente de celle donnant lieu à la distribution 
grand-canonique dans laquelle N fluctue autour d’une valeur moyenne N suite à des 
échanges aléatoires avec un thermostat. La variation de N envisagée ici est au contraire 
contrôlée et non erratique. Dans une telle transformation où N varie, la variation de l’énergie 

r dF^ 


interne est : d U = dW +dQ + pdN ou p 


\dN jy 


est le potentiel chimique, c’est à dire , la 


force généralisée associée au paramètre extérieur N . 

La quantité pdN est le travail reçu par le système pour faire varier le nombre de ces 
particules de dN . L’énergie qu’il faut fournir pour introduire une particule supplémentaire 

dans le système comporte d’une part le travail ^ e t la chaleur dQ=T . 

dN dN 

La pression P d’un gaz exprime la tendance du volume à augmenter. De manière analogue, le 
potentiel chimique exprime la tendance des particules à sortie du système, ou à y être 
absorbées dans une réaction chimique. 

Exemples : i) Gaz parfait classique isolé : le potentiel chimique est l’énergie moyenne 

- - ôS 

(U = N s) moins la chaleur par particule : p - s - T — . 

dN 
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ii) Le Gaz d’électron à T -0K : Le potentiel chimique est égal à l’énergie de Fermi E F (0) . 

Potentiel chimique et énergie de Fermi : En distribution canonique, dans une 
transformation sans variation des paramètres extérieurs autres que le nombre de particules N, 

la variation de l’énergie interne est: dU = dQ + fidN . Donc : dS- - ^dN 

( dS 

supposant la transformation réversible). D’où : fi = -T\ — 

l dN 


( en 


Sx, U 


En distribution grand-canonique le paramètre a est défini par : 


1 ds(E, N ) 

a = 

k 8N 


En admettant l’équivalence entre ces deux distributions pour N- IV , on a : a = fdfi . 

Le potentiel chimique est alors identifié à l’énergie de Fermi introduite dans la statistique de 
Fermi-Dirac et de Bose-Einstein (voir chapitre suivant). 


11°) Grand-Potentiel 

Soit un système en équilibre avec un thermostat avec qui il échange aussi des 
particules. On doit utiliser la distribution grand-canonique. On définit alors le grand potentiel 


analogue à l’énergie libre F par : Q(x, F, //) = -kT LogZ c (x, T , ju) = U -TS - juN 


où x est un paramètre extérieur. L’équivalence entre distributions canonique et grand- 
canonique implique: 


Q 


(x, T, fi) = F (x, T,N)-juN 


G(x, T, /u) est un potentiel thermodynamique, donc : ùQ< Xdx - SdT - Ndju (l’égalité 

n’a lieu que pour une transformation réversible). Dans une transformation isotherme (dT = 0) , 
isochore(ùE = 0) d’un système en communication avec réservoir de particules (dfi = 0) , le 
grand-potentiel O ne peut augmenter. L’équilibre correspond au minimum de O. Ceci 
généralise le cas des transformations isothermes, isochores et à nombre constant de particules, 
pour lesquelles l’énergie libre F est le potentiel thermodynamique. 

Q(x, T, /u) est aussi une fonction caractéristique: X = 



5 = 


l Sx , 

T,» 

U T; 


iV = - 




dfi 


et U = Cl -T 


ST, 


ÔQ 

dT 


F 


r dQ} 


dfi 


ST 


12°) Expressions intégrales des fonctions thermodynamiques 

P grandeur intensive, i.e. indépendante du volume. Or : P = - 


D’où les expressions intégrales : 
G- fiN ; et H = TS + fiN . 


-PV 


Cl=-PV ; F - -PV + fiN ; U = TS - PV + fiN ; 


13°) Inégalités thermodynamiques 

Soit un système tel que : F = este et T = este . Pour toute variation à partir de 
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l’équilibre stable, on a dG = dU - TdS + PclV >0 (l) . Développons au second ordre en 5 et 


, j ... dU dU 1 

V : dU — dS dV -i — 

dS ÔV 2 


d U dS 2 + ^rdV 2 + 2— M 


dS 2 


dV- 


dSdV 

52 t 


.Or: ™=T. 

dV dS 


dU dU dU 

Remplaçons dans (1) dU par son expression, on a : — -dS 2 + 2 dSdV + — - dV 2 > 0 . 

dS 2 dSdV dV 

Or, la théorie des fonnes quadratiques implique pour qu’une forme comportant deux variables 
dS et dV soit toujours positive, il faut que ses coefficients vérifient les deux conditions : 



d 2 U 

d 2 U 



d 2 U 

d 2 U 


dP 

dP 

i) 

dV 2 

d 2 U 

dSdV 

d 2 U 

> 0 et ii) 

Oj 03 

^ cT 

V 

O 

O 

dV 2 

d 2 U 

dSdV 

d 2 U 

0 

O 

A 

dV 

dT 

dS 

dT 


dVdS 

dS 2 



dVdS 

dS 2 


dV 

dS 


< 0 . Donc : 


i) 


d 2 U dT 


dP dT dT dP 


dP dT 


= — = — > 0 , ii) < 0 => < 0 


dS 2 dS C 


dV dS 


dV dS 


dT dP 
dV~dS 


f dT^ 


K dV; 


< 0 (relations de Maxwell). Donc obligatoirement : 


dV dS 
f dP ^ 


car 


\dV; T 


< 0. 


Résumé : (i) =>C V >0 , ceci résulte en fait de la concavité de l’entropie, (ii): 


r dp^ 


\dv j T 


< 0 : la 


croissance du volume à T= Cste s’accompagne toujours d’une décroissance de la pression, 
(iii) Les états d’équilibre qui ne satisfont pas ces deux inégalités sont instables. 


14°) Relations thermodynamiques 

Soit un système dont l’équation d’état est caractérisée seulement par deux variables 
d’état : (P , V) ou ( V , T) ou (5 , V ) ou (5 , P) . 

a- Pour chaque type de transformation correspondent différents coefficients de réponse : 


Coefficients de dilatation : 


f dV^ 


a 


V 


Coefficient de variation de la pression : (3 = - 

ï 

1 f dV^ 

Coefficient de compressibilité : x = ~~ 


\8T j p 
1 f dP ' 


a 


isoentropique 


r dV^ 


V 


\dT Js 


dT 


Jy 


Chaleurs spécifiques : C v = 


f dQ ' 

\dT J y 


= T 


\dP j T 
r dS^ 


\dT J y 


et C P = 


r dQ} _ T ( dS^ 


\dT j P 


\dT j P 


b- Les 1 er et 2 eme principes imposent des contraintes sur différentes transformations 
(isobare, isochore, isotherme, isentropique,...)qui se traduisent par des relations entre ces 
coefficients de réponse. Formellement ces relations sont obtenues en imposant que les 
différentielles dU , dS , . . .soient des différentielles totales exactes ( d.t.e). 

Exemple : Les conditions de d.t.e pour dU , d F , dG et dH donnent les relations de 
r dS^ ( dP\ ( dT') f dP) (dT') ( dV') , ( dV') ( dS^ 


Maxwell : 


\dV j T 


\dT J y 




\dS J y 


\dP J s 


\ dS j P 


et 


\dT j P 


\dP j T 


Une relation de C p - C v est obtenue sachant que: dQ=C p dT + T 


f dS^ 


\dP j T 


dP =C,,dT + IdV 


avec: dP = 


r qp\ 
\ST Jy 


dT + 


( Qp \ 

\3VJt 


er 


dV . Alors : C P - C v = PVT afi ou encore : C P - C v =VT — . 

X 
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Chapitre 5 : 

DI STRI BUTI ON DES VI TESSES DE MAXWELL 


-I- Introduction : Modèle du gaz parfait. 


• Considérons à T ordinaire, un gaz rare dilué. Ses N atomes se comportent comme des 
objets sans aucune structure et leur potentiel d’excitation est très élevé. Ils restent 
gelés dans l’état fondamental. Les collisions entre atomes sont des chocs élastiques. 

• Ces atomes sont représentés par des points matériels de masse m (sphères dures 
=> chocs élastiques). 

• La statistique classique est une bonne approximation puisque le gaz est très dilué et la 
température T est élevée. 

• Le libre parcours moyen (lpm) d’un atome est grand devant ses dimensions : 
lpm ~10“ 4 cm à T = 300 K et P - 1 atm . Les collisions sont ainsi rares. On peut 
négliger le potentiel d’interaction interatomique devant l’énergie totale du système. 

• Cependant, le nombre de collisions (~10 9 / sec /molécule dans l'air ) est suffisant pour 

permettre au gaz d’atteindre l’équilibre rapidement. 

• Les chocs des atomes sur une paroi permettent des échanges d’énergie avec le milieu 
extérieur. 


II- Distribution de probabilité dans l’espace des phases 


1°) Densité en phase 

Le gaz est supposé en équilibre thermique avec un thermostat. Il y a échange d’énergie 

N / 2 , _ 

=^> distribution canonique à l’approximation classique. L’énergie : £ = ^ /2m + V ( r ) 


Habituellement V (r ) représente l’action des parois : V (r j 


fo, 

loo. 


1 


i = 1 ' 

à l'intérieur 
ailleurs 

N 


Dans l’espace des phases, la densité de probabilité est : pdY=—e ~ pE d " Pi / Nih ' N 


i = 1 


2°) Calcul de la fonction de partition Z 


Z = 


\ - ( j dp,) 3 - [Inm/PŸ ; 


V= Volume du gaz 

car : 


J e " r dx = yjn/a . Donc : 


r 7 r' 

Z = — -z 
NI 


avec : z=-^j ; X- hj \j(2xmkT) ; où z est la fonction de partition pour un 


À 


seul atome et À est une longueur de l’ordre de la longueur d’onde associée aux atomes du gaz. 
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III- Distribution des vitesses de Maxwell. 


Considérons le cas d’un seul atome (N = 1) en équilibre avec le reste du gaz, qui joue le rôle 
du thermostat pour l’atome considéré. Pour cet atome, on a : 


pdT=-e~ pe 

z 


d 3 r d 3 p 
h 3 


avec 


— 2 



1°) Densité de probabilité des impulsions et des vitesses 

A N = \, Y intégration sur r donne: f(p)d 3 p oc e 

Par un changement de variables on obtient : / (v) d 3 v 


~ pp2/2m d 3 p. 


C n \ 3 / 2 

pm 

2 n , 


- P p 1 jim 


d 3 v 


2°) Densité de probabilité de la position 

-Mr) 

L’intégration sur p donne: f[r)d 3 r = — — —d 3 r. 

je ' r ' d 3 r 

Elle est constante dans la bouteille contenant le gaz et nulle en dehors. Cependant l’expression 
obtenue est aussi valable pour un potentiel des forces extérieures quelconques, par exemple la 
pesanteur. On obtient ainsi la formule barométrique : f[rj oc e~ Pmx si v[r^ = mgx est le 

potentiel de pesanteur. Cette dernière formule donne la variation avec l’altitude x de la 
densité de l’atmosphère supposée en équilibre thermique. 

3°) Autres distributions dérivées 

a- Distribution d’une composante de la vitesse : /(v.) dv, = (ftrn/27rj e^ p,n '^ 1 dv_ 
Cette expression est obtenue en intégrant sur v x et v v . 
b- Distribution du module v de la vitesse ; 
c- / (v) <A= ' - e- pm ' pl Anv 1 dv : loi de Rayleigh, 

f (v) dv — r ^ 0 . - > 0 et passe par un maximum pour v = v + valeur la plus probable. 



Applications : v 2 = 3v„ 2 ; v = (8 /zrm/?) 1/2 ; v + = (2 /ntfif 2 
s c = p 1 jim = 3/2/? : l’énergie cinétique moyenne ne dépend pas de la nature du gaz 
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00 1 _ 1 _ 1 00 

Formulaire : I n =j x n e~ a * 2 dx= — T (- — ); jx n e~ f,x dx=n\ ii <n+l) =T{n + \) n~ (n+v> -, n*0 

r sont des fonctions de Riemanntq: r(x + l) = xr(x) ; x>0 de valeurs r(3/2) =4nj2 et 
r(5/2)=3V^/4 . 

La vitesse moyenne est du même ordre que la vitesse du son : v son = (y / Pm f~ où 

y = C p jC v ~ 1 . Cette vitesse est d’autant plus grande que la masse m est plus petite. 

Dans une atmosphère planétaire, une proportion assez importante de molécules légères ( H 2 ) 
peuvent avoir la vitesse de libération ( 11 km/ s sur terre). Ceci explique la disparition de ces 
molécules. C’est ainsi que toute atmosphère a disparu sur la lune. 


IV- Notion de théorie cinétique des gaz. 


En considérant le mouvement des molécules, on peut étudier un certain nombre de 
problèmes liés aux gaz. En particulier dans le cas proche de l’équilibre ; on suppose les 
vitesses réparties selon la distribution de Maxwell : C’est la théorie cinétique des gaz. 

1°) Effusion : L’effusion est l’échappement des molécules à travers un trou de la paroi 
suffisamment petit pour ne pas perturber la distribution des molécules, même près du trou, 
(celui-ci doit avoir un diamètre <s: au libre parcours moyen). 


Calcul approximatif du flux é à travers un trou 

Les molécules de vitesse v traversant un 
trou de surface dS pendant dt sont 
contenues dans un cylindre de volume : 


de la surface : 

v dS dt . Si n est le nombre de molécules 
par unité de volume (densité) ; 1/3 d’entre 

Vdt 

1 — 0 ds 


elles ont à peu près leurs vitesses normales à dS dont la moitié, soit 1/6 se dirigent vers le 
N n - 

trou. Donc è = , avec N = —vdSdt, s’écrit alors : 

dS dt 6 

Ce flux dépend de la masse m des molécules. Ce fait est à la base d’un procédé de 
séparation isotopique (enrichissement de l’uranium en isotope 235). 



2°) Pression cinétique 

• Une molécule frappant la paroi et renvoyée par celle-ci, exerce une force sur cette paroi. 
La force moyenne par unité de surface est la pression exercée par le gaz sur la paroi. 

• Le choc d’une molécule sur la paroi conserve son énergie cinétique moyenne et donc son 
impulsion. Celle-ci est mv avant choc, elle devient - mv après choc, soit une variation de 
A p = 2 mv 

T ■ . , , „ F „ „ Ap , , r , 2 mv 

• La pression est égalé: P = — . Or F =- £ - et alors P = . 

ds t ds t 

— 1—2 

• Pour N molécules, la pression est N P soit: N P- 2mv</> - -nmv <=> 

Ce résultat est indépendant de la distribution des vitesses et de la réflexion sur les parois. 
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V- Théorie cinétique des phénomènes de transport. 


Le but est la description statistique simple de quelques états stationnaires de transport, 
tels que l’écoulement d’un gaz, les conductivités thermique et électrique d’un gaz. 

Un état stationnaire est caractérisé par des grandeurs macroscopiques constantes dans le 
temps. Cependant, il n’y a pas équilibre thermique car il y a transfert de matière, de chaleur 
ou de courant au milieu extérieur qui n’est donc pas en équilibre avec le système considéré. 
Nous considérerons seulement le cas d’un gaz dilué tel que les approximations suivantes 
soient vérifiées : 

i) La durée d’une collision est négligeable devant l’intervalle de temps entre deux 
collisions. 

ii) Les collisions à plus de deux molécules sont négligeables. 

iii) L’ approximation classique est valable (A «: d) . 


1°) Temps de collision et libre parcourt moyen 

• P k (?) = probabilité d’observer k collisions dans un intervalle de temps t. C’est une 


distribution de Poisson : 



Car : i) les instants de collisions sont indépendants entre eux ; 

ii) la probabilité d’observer une collision dans dt infiniment petit est Adt (infiniment petite 
de 1 er ordre) ; 

iii) la probabilité d’observer plus d’une collision est un infiniment petit d’ordre supérieur. 

P 0 (?) =<? ~ À ‘ = probabilité que la molécule ne subisse aucune collision. 


• Le nombre moyen de collisions par unité de temps est A . 

Dans un état stationnaire, A ne varie pas, mais dépend généralement de la vitesse de la 
molécule. 

• Soit Q(t)dt= la probabilité que l’intervalle entre deux collisions d’une molécules soit 
compris entre ? et dt : Q (?) dt = P 0 (?)/} {t)-e ~ h Adt . Soit alors : Q(t)= Ae~ Xt . 


• Temps de collision r= intervalle moyen entre deux collisions : r=[?<2(?)d? = — ; 

0 ^ 

est aussi appelé temps de relaxation. 

Alors: Q(t)dt =e~ ,/r d(t/ r) . Soit / = v? chemin parcouru par une molécule de vitesse v 
entre deux collisions : Q(l)dl =e ~ l/l ° d (///„ ) où : 1 0 =vt est le libre parcours moyen. 


2°) Section efficace de collision 

• La probabilité p de collision d’une molécule sur un parcours d est proportionnelle à la 
concentration n du gaz et à la distance d parcourue : p = a n d . a est la section efficace 
de collision. C’est encore la probabilité de collision quand l’épaisseur de matière traversée 
comporte en moyenne une molécule par unité de surface (nd =1 ). 

Remarques : (i) cr a la dimension d’une surface (c’est la surface plane que présente une molécule 
cible) ; (ii) Les molécules "cibles” sont supposées immobiles. 

Modèle des sphères dures : Une sphère de par unité de sur f aC e : a = n (a, + a 2 ) 2 . Soit 
rayon a, rencontre une sphère de rayon a 2 
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dans le cas où 
cr = nd 2 . 


d 

a, = ci-, =— , on a : 
1 ~ 2 



• Section efficace et temps de collision : 

Calculons le nombre moyen A de collisions d’une molécule considérée comme projectile, 
animé d’une vitesse v, . Considérons d’abord les molécules cibles, de vitesse v 2 . Dans le 
système au repos des molécules cibles, la vitesse du projectile est : v, + v 2 . 


L’épaisseur du gaz traversée par le projectile et par seconde est 
La densité des molécules cibles étant : n f (v, ) d 2 v 2 , donc : 


v, a 


v, - v. 


le nombre de collisions d’une molécule projectile par seconde est : n f ( v 2 ) d 2 
Intégrant sur v 2 et moyennant sur v, , on a alors le nombre moyen A de collusions par 


seconde: A - n j f (v l ^jd 3 v l J/(' ; 2 )d 3 v 2 1 

de Maxwell, on obtient : A = n (fim/lxf J llv, - v. 


. Utilisant la distribution de vitesses 


^ m(v ’ + l -)/ 2 ad 2 v { d 2 V 2 . 


Soit : A =V2cr nv = r 1 et l 0 -v t-—= 

V2 


■L - 


1 


cr n 71^2 nd 

Exemple : Pour l’azote à T = 300 K et P=\atm, on a : d - 2.10 ^ cm ; 


(modèle de sphères dures). 


n - 2.4.1 0 19 cm 3 ; v = 5.10 4 * * * * cm/s ; cr =1.2.10 15 cm 1 ; A =2.10 9 s 


l n - 3.10 cm . 


On vérifie sur cet exemple que la durée d’une collusion [-d/ v ) «: r : (temps de collision) 
ou intervalle moyen entre deux collisions. 


3°) Viscosité 

a - Définition : Le phénomène de viscosité est caractérisé , pour un écoulement permanent 
d’un gaz dans un tube, par : 

• La vitesse d’écoulement est plus grande au centre du tube . 


z* 




-> 


-f 


Ux 


4 — 

ir 




• Cette vitesse est nulle au voisinage immédiat des parois. 

Pour simplifier, considérons un écoulement parallèle au plan Oxz avec une vitesse 

d’écoulement u ( z j dépendant de l’altitude z . La couche de gaz située en dessous de 

l’altitude z exerce une force de ralentissement sur la couche au dessus. 


Soit n cette force par unité de surface (pression)normale à Oz. Au premier ordre, cette force 
est proportionnelle au gradient de la vitesse ( = frottements) : 


du. 


TC 


OZ 


. tj est le facteur 
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de viscosité du gaz. Le mouvement d’une couche située entre z et z + dz : La résultante des 
forces exercées par le gaz en dessous et au dessus est : n + dn - n , soit : 

d 2 u 

dre = - // — f- dz . La condition de stationnarité : dre = 0 (pas d’accélération) 


d\ 

' dz 2 


dz 


= 0 : gradient des vitesses uniforme. 


b- Calcul approximatif de 77 : 

re résulte du transfert des molécules au travers du plan horizontal de côte z par 
agitation thermique. La quantité de mouvement des molécules venant d’en bas, qui est en 
première approximation : mu x (z - /„), est inférieure à la quantité de mouvement de celles 
venant d’en haut : mu x (z + / 0 ) (en l cie approximation). 

Supposons que le gaz est presque en état d’équilibre thermique, ceci suppose que la 
vitesse d’écoulement est très petite devant la vitesse moyenne des molécules à la température 
considérée (toujours réalisée au températures ordinaires). Le libre parcours moyen doit être 
aussi très petit devant l’enceinte contenant le gaz de façon à pouvoir considérer l’écoulement 
comme un phénomène collectif et non comme une diffusion incohérente de molécules. 

Alors le nombre de molécules venant du bas ou du haut est par unité de temps : nv/ 6 .Ce qui 
donne pour ;r : re - nm v\u x (z - /„) - u x (z+ Z 0 ) ]/ô 

Soit donc en comparant avec la force re de la définition : 


n 


nm v 


h h 


ou : 77 —nkT t . 


ou encore : 




Remarques : i) 77 ne dépend pas de la pression du gaz et de sa concentration n. Quand 
n diminue, il y a moins de molécules traversant le plan d’altitude z ; mais celles-ci proviennent d’une 
distance l 0 plus grande et ont donc un excès de vitesse plus grand. Ces deux effets se compensent, 
ii) 77 augmente avec T, contrairement au cas des liquides. 


4°) Conductivité thermique 

a- Considérons le régime permanent de transport de chaleur dans un gaz. Supposons 
que le gaz soit localement en équilibre thermique ; mais que la température croit avec 
ÔT 

l’altitude z : — >0 . La chaleur s’écoule du haut vers le bas. Au 1 er ordre, le flux a de 

dz 


chaleur est proportionnel au gradient de température : 


ôT 

Cl — — K 

dz 


. k est le coefficient de 


conductivité thermique du gaz. La condition de stationnarité d’une couche comprise entre z 


et z + dz est: clq = 0. Soit encor : 


ô 2 T 

~ôë 


= 0 : le gradient de température est uniforme. 


b- Calcul approximatif de k : 

Le nombre moyen de molécules venant du bas et traversant un plan de côte z par unité de 

temps est : — nv . Ces molécules apportent une énergie moyenne : e(z- /„) . De même, le 

nombre moyen de molécules venant d’en haut et traversant un plan de côte z par unité de 

temps est : — nv . Elles apportent une énergie moyenne : c:(z+ /„ ) . Le flux est donc : 

6 
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q 

a : 


1 -r- 

— n 
6 


v[s(z- /„) - s(z+ /„)] = - j n vl 0 • En comparant avec ce qui précède, on 


1 -/ C 
k — — n v L — 

3 N 


C Ô£ 

. — - — est la chaleur spécifique du gaz par molécule. 

1 c 1 î - c 1 

Remarques : (i) k = - nv = = - v =- est indépendant de la pression et de la 

3 N V2 ni y 3 N usjl 


concentration, (ii) K = f(sfr j : K croit comme \ff lorsque T augmente, (iii) 


k 3 k 

il 2 m 


5 °) Diffusion 

a- Considérons par exemple, la diffusion dans un gaz de molécules dites marquées ; 
c’est à dire, ayant un noyau radioactif. On suppose négligeable la différence entre la masse de 
ces molécules et celle des molécules ordinaires du gaz. Supposant en plus que la 
concentration n des molécules marquées dépend de l’altitude z seule. 

Au 1 er ordre, le flux J de ces molécules est proportionnel au gradient de la concentration : 



où D est le coefficient de diffusion. 


Considérons non pas un régime permanent, mais l’évolution vers l’équilibre thermique, n 
varie donc avec le temps et évolue vers une constante uniforme. La conservation du nombre 
de molécules pour une couche de gaz comprise entre z et z + dz implique l’équation: 


d 

— (ndz) - — dJ 
dr ’ 


dite : équation de diffusion. - dJ : 


élément de flux provoqué par la 


diffusion et —(ndz) : élément de flux provoqué par l’évolution dans le temps. On en déduit 


que la concentration obéit à l’équation différentielle : 


d n p d 2,? 

ôt dz 2 


b- Calcul approximatif de D : 

• Le nombre moyen de molécules venant du bas et traversant un plan de côte z par unité de 
temps est : vn(z-l 0 )/ 6. 

• Le nombre moyen de molécules venant d’en haut et traversant un plan de côte z par unité 
de temps est : vn(z+l 0 )/ 6. 

• Le flux est donc : v \n(z- l 0 ) - n(z+ l 0 ) ] /ô . 


On en déduit alors : 



Rqs : (i) D est inversement proportionnel à la pression ou à la concentration n. 


(ii)A pression constante D varie comme T car v= —kl 


varie comme 


Vr et 


p * 1 
P = nkT = este => n = — et /„ = —= — 

kT sJ2 n 


' k ^ 


V2 


Pu 


K m J 

l n varie donc comme T. 


(iii) — = — = — ( p : masse spécifique du gaz). 

rj mn p 


6°) Conductivité électrique 
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Considérons le déplacement des ions dans un gaz sous l’effet d’un champ électrique E . Au 1 er 
ordre, la densité de courant j est proportionnel à E : 

électrique. Le champ E accélère les ions entre deux collisions. La vitesse s’exprime en 

- e~E 

fonction du temps 6 écoulé depuis la dernière collision : v = — 6 + v 0 où 0 a la même loi de 

m 

probabilité que le temps entre deux collisions (vu précédemment). 

Supposons que chaque collision restaure en moyenne les vitesses d’équilibre : vo = 0 , alors : 


j = oE 


où cr est la conductivité 


2 

— eE ~ ne 

v = — r => j = nev = Et= ctE 

m m 


=> 


2 

ne 

cr = t 

m 


: Formule de Drude 


VI- Equation de transport de Boltzmann. 

La fonction de distribution / [r, v, tj peut subir des modifications dans le temps sous 

l’effet de plusieurs facteurs : effet d’un champ (de force), de diffusion, de collision,... 
Supposons qu’à un instant t la particule de vitesse v est en position r ; elle est représentée 

dans l’espace des phases par le point (r, vj entouré de l’élément de volume d 3 rd 3 v . 

A , t' = t + dt elle est en (r', v'J entouré de d 3 r'd 3 v' tel que : 

—, - - dr , - - , — , - dv , - F , 

r = r+ dr = r H dt- r +vdt => v = v + dv = v H dt =v H — dt 

dt dt m 



Or, le mouvement dans le temps d’un système est une transformation canonique infinitésimale 
de la configuration de ce système dans l’espace des phases (voir cours de Mec. Analytique : 
P3). Donc le volume d 3 rd 3 v de cet espace est conservé sous cette transformation. 

Alors : d 3 rd 3 v = d 3 r'd 3 v' . On a ainsi deux cas: 


a-En absence de collisions, le nombre de particules au voisinage de la particule considérée 


dans le temps est conservé : fyr',v',tj -f(r,v,tj 


d 3 rd 3 v = df d 3 rd 3 v = 0 


df . df dr , df dv , 
<=> — dt + Ar — dt + Ar — dt 
dt dr dt dv dt 


d 3 rd 3 v = df d 3 rd 3 v = 0 . 


<=> 

df * — • — 

— + r V r f + v V,./ 



dt 

dt 


Equation de Boltzmann sans collisions . 
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b-En présence de collisions la fonction de distribution réduite f[r,v,t ) est modifiée aussi 
dans le temps car ces collisions font apparaître ou disparaître des particules dans d 3 rd 3 v 


— ^ 0 . Donc : 
dt 


df * — • — 

-±- + r V r f + v Vvf 
dt 


v dt j coll 


: Equation de Boltzmann générale. 


Rq : (i) L’équation de Boltzmann est un cas particulier de l’équation de Liouville obtenue 
pour la densité totale : p(r,v,t\ ; car l’équation de Boltzmann est obtenue pour la densité 

réduite : f[r,v,t \ . 


(ii) Le terme de collisions 


peut être obtenu en utilisant des hypothèses 


s coll 


d’approximation. Dans l’hypothèse du temps de relaxation r : 


d L 

8t Jcoll 


f-f 


ou 


f° est la densité réduite de l’état d’équilibre et / est celle de l’état hors d’équilibre. 
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Chapitre 6 : 

SYSTEMES DE RARTI CULES I DENTI QUES ET I NDEPENDANTES 

I- Introduction 

• Il existe des systèmes formés de particules identiques faiblement couplées tels que les gaz 
parfaits, les électrons de conduction des métaux, l’hélium liquide, les photons, les phonons 
des vibrations des solides,.... Ces particules sont représentés par des points matériels de 
masse m (sphères dures => chocs élastiques). 

• L’identité des particules (fermions ou bosons) conduit à tenir compte du postulat de 
symétrisation de la mécanique quantique . 

Fermons : électrons, protons, neutrons, v , mésons // , He i , quarks,. . . 

Bosons : photons, mésons /r et k , noyaux paires, He 4 , H 2 , gluons, . . . 

• Les conséquences de l’indiscernabilité sont importantes lorsque il y a recouvrement des 
fonctions d’onde des particules. 

=> A la limite classique : pas de recouvrement, les particules sont discernables. 

• Dans certaines conditions l’indiscernabilité peut être levée. 

Exemple : Dans le cas des ions de Cu 2+ dans un sel paramagnétique, la localisation de 
ces ions sur les sites aboutit à un très faible recouvrement des fonctions d’onde d’espace. Les 
ions sont alors discernés par leur localisation : la mécanique quantique se manifeste seulement 
par quantification du spin et du moment magnétique. Ce cas est décrit par la statistique de 
Maxwell-Boltzmann. 

Les particules indiscernables seront décrites par les statistiques quantiques de Fermi- 
Dirac et Bose-Einstein. 


II- Particules identiques discernables et indépendantes : Distribution de 

Maxwell-Boltzmann (classique) 


1°) Introduction 

Un système I à N particules, son hamiltonien est : H =jf (l, N^j =^/i, (/) +^V iJ ( i,j ) 

i i j 


où li = h, V i (identiques) , VL est supposé négligeable et i = ensemble des coordonnées 

d’espace et de spin de la particule i (il peut contenir d’autres nombres quantiques pouvant 
spécifier la particule i tels que les moments orbitaux internes,...). Les états propres | y/) de 


K sont construits à partir des états propres 


Va' 


(7)) de 


h : — h : h 


Va 


P) 


Va 


(-))■ 


Une complexion de I est donc : y/ (I, 2, , = y a, (ï )} ® Va , (2 'jj ® .... ® \// ?s [n . 


L’état quantique de chaque particule étant parfaitement connu. La valeur propre 
correspondante est : E = ^ e x . Plusieurs particules peuvent être dans le même état. Soit 
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N À le nombre de particules dans l’état \y/ A ) ( nombre d’occupation de l’état y/- ) ) d’énergie 
Ci , alors : N = £n â et E = Ci • 

Définition : t/ne configuration est par définition l’ensemble { Af } des valeurs 
possibles de N. . Pour une configuration donnée on connaît seulement le nombre N À de 
particules dans y/. ) , mais on ne connaît pas lesquelles parmi N qui sont dans cet état. 

Propriété Pour une configuration [N À } il y a W complexions possibles telles que : 





Démonstration qualitative : N\ est le nombre total de pennutations des N particules entre 
elles, qu’il faut diviser par N À ! permutations de particules qui sont dans le même état | y/f) 
puisque ces N. ! pennutations ne donnent pas une complexion de I différente. La division 
est appliquée pour tous les états \y/ x ) possibles, d’où l’expression de ff({AC}) ci-dessus. 


2°) Distribution de Maxwell-Boltzmann (M-B ) 

a- Cas d’une seule particule 

La particule est supposée en équilibre avec un thermostat, on utilisera donc la distribution 
canonique. La probabilité d’occupation de y/fj d’énergie e x est donnée par la loi de 

distribution de M-B: = f (s x ) = -e ~ PSl avec : z = ^ e ~ Psi . 



• La superposition incohérente des états \y/ x ) a pour opérateur densité : p = ^ f | y/, ) (i// A 

A 

• Valeurs moyennes : s = 'Y J f À £ À ; £ 2 =^f À £ À 2 . 

A A 

• p = l e -^ avec h = J^£ À \y/ À )(ip À \ . 

Z A 

• La connaissance de z permet de donner : s = î cr 2 (f) = - ^ et s=kLogz + kj3s . 

Remarque : Z porte sur tous les états i// Â ) ; elle peut être remplacer par E à condition de tenir compte 

A ^ 

de la dégénérescence g (yy ) des niveaux £ x . Par exemple : 


1 e x 


b- Système de N particules 
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a ) Les N particules sont indépendantes, donc la probabilité d’une complexion P(C ) 
est égale au produit des probabilités d’occupation de chaque particule : 


P(C)=-e 

Z 


Z 


.-g-*'» 

Z 


<=> 


P(C)=fe 

Z 


-PE 


avec : 




P ) L'opérateur densité de I : p = p x ® p 2 ® .... ® p N ; soit : 


1 


1 


1 


— B h, 1 —B lu 1 —B h, 

p = — e p 1 — e F * — e p 3 


1 


,-PK_ „ - 


1 -AÏ," 


=P= Z 


donnant ainsi les grandeurs 


Z Z Z Z 

N 

suivantes: La fonction de partitionZ = z N ; l’entropie S = '^ j s i = N s et l’énergie de 

î =i 

S : E - N s telle que cr 2 (£) -Na 2 (s ) . 


y ) Expression de P(C ) en fonction des nombres d’occupation N À . Posons 


alors: P(C) = «xf 2 •....• x N f •...<=> 


fW=7ll 


z=zn 


PM 1 


1 


*Valeur moyenne du nombre de d’occupation : Ni = ^ P(C)N À =— ^A^x 2 

c Z c 


N, N, N, 


1 5 

X , 

Z dx 


- y, xf • x^ 1 x N f «. . .Soit: 


1 c 


Ni = X; 


dLogZ 


Sx, 


iV 


. Soit encore : Ni - — e : M-ZÎ. 


*Fluctuations de N, : 


(V.) = 


ô N , 


Sx, 


soit: a 2 (N À ) = N f x ( 1- f\ ). 


Conclusions :i) La population Ni d’un état y/ - ) décroît exponentiellement avec l’énergie. 
ii)Les fluctuations relatives de cette population deviennent négligeables pur N grand. 


II- Particules identiques et indiscernables 
- fermions et bosons - 


1°) Introduction 

a) Principe de Pauli : La fonction d’onde de N particules doit être symétrique ou 
antisymétrique, par rapport à toute permutation de ces particules entre elles. 

• Elle est symétrique pour des particules de spin entier : Les bosons. 

• Elle est antisymétrique pour des particules de spin demi-entier : Les fermions. 

b) Conséquences du principe de Pauli : 

i) Deux fermions ne peuvent être dans le même état => les valeurs possibles de N À sont : 

fO, 1 pour fermions 
À [O, 1, 2,... pour bosons 

ii) Dans les deux cas, il y a une complexion {Af }) pour une configuration \ N À } donnée. 
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En symétrisant les NI fonctions d’onde y/ (1,2, , de S, on obtient les fonctions 

d’onde UAfJ) de N particules (bosons ou fennions) : 

iR})=h«o ^ o^dpz^ ^ ’ ^)H pn^ ] 2 z p K 1 2, , n ) 

p 1 y / à J p 

|{^}> = k^) ,ÏŸ))=(l/iV!)^^ P H ï ’2, ,N) 


P ' P 

P désigne les permutations des N particules et s p = ± 1 désigne la parité de la permutation P . 


2°) Statistiques de Fermi-Dirac et de Bose-Einstein 
a- Probabilité d’une complexion. 

Nous utilisons ici la distribution grand canonique (pour couvrir le cas général où il y a 
échange d’énergie et de matière avec le milieu extérieur). 

• Les complexions |C)= |{Afi}) forment la base de l’espace de Fock de I ; leur 

1 


probabilité est : P(C ) = — e 
Posons = e a ~ PSi , P(C ) s’écrit aussi : 


ou 


N = et K = y.'V.c, . 


p ( c )=^n 


-‘G A 


• La fonction de partition est : 

- P e *) N . 


z g = no + r > x x)' 


avec : r) = 


\ 1 pour fermions 

I - 1 pour bosons 


En effet : Z c = ^ e 


N, N 2 AT, 

Or: 


-z n 

N/ 


f \ 

N 


soit : 


va y 


f \ 

Z<‘ 

r \ 

Z< 2 


f V 

z^ 

=n 

r \ 

z^ 

K N > y 

K N 2 J 


l JV. J 

a 

y 




Finalement : 


fl + x k pour fermions (N- =0, 1) 

1(1- x x ) ' pour bosons (N À =0, 1,2,3,...) 

Z G — n(l P X x) aV OC P F -D , B-E —il • 


b-Valeur moyenne des nombres (V occupation N , - Fonctions de distribution de 

Fermi-Dirac et de Bose-Einstein : 

On remarque que P(C) en fonction de N À a la même forme que P(C) en distribution 

dLog Z 


de Maxwell-Boltzmann => N À aura la même forme : N a - x k 

n ^= x a-^—Yj L °s { 1 + 7 7 x a)" =P x a-^—Yj L °s { 1 + P x a) 

Sx À a dx k À 

P _ X A _ 1 


ôx. 


. Soit : 


p x x - - - , 

1 + P X A 1 + P X A X A + P 


Soit encore : Nx=f(e x ) = 


ou : x x = 


A 


1 + P X a 1 - P f A 

Les distributions de Fermi-Dirac et de Bose-Einstein : f(s)~ 


p( e - fi) 
e y ’ + p 


■> Pf-l 


= ± 1 
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Remarques : i) f (s) joue le rôle de la distribution des vitesses pour un gaz parfait quantique. 
p = ~p = £ F ( 0 ) + °{t 2 ) (pour métaux et semi-conducteurs). s F : énergie de Fermi. 

ii) Toutes les grandeurs thermodynamiques de I s’expriment en fonction de f k , puisque 

LügZ G = îjY J Log{ 1 +ï]X x )=-r 1 Y J Log(\ -rjf k ) . 

2 1 

c - Statistique de Fermi-Dirac : p = + 1, donc : f k = , } . = f(s k ) 

e + 1 

Propriétés de f{s k ) : 

a ) Soit P ; (l) = probabilité d’avoir une particule dans | y/ k ) alors : 

fx = Nx = £^ = 1.^(1) + 0.^(0) = (1). 

i 

Ainsi f À est aussi la probabilité d’occupation de l’état \y/ k ). 

P) f [s) présente un point d’inflexion F à l’énergie de Fermi : /u =s . 

Le domaine de / [s ) décroît rapidement à une largeur de 2 //? = 2 kl . 

Définition : On définit la température de Fermi T F par : p - k T /: . 



A T = 0 K , tous les états | y/f) tels que : s k < p sont occupés (=> P k (l) = 1 ) et tous les états 
| y/ k ) avec s k > p sont inoccupés ( => P k (l) = 0 ). 

Quand T/', le gaz de Fenni reçoit de l’énergie et les états | y/ k ) d’énergie s k < p se 
dépeuplent au bénéfice de ceux d’énergie e k > p . C’est ce que traduit l’aplatissement 
du genou de f{s). 

y ) Symétrie particule-trou d’un gaz de Fermi : L’absence d’une particule dans | y/ k ) 
est interprété par la présence d’un trou dans cet état (et inversement). Or f{s) est symétrique 

par rapport au point F [p, V.) . En effet, si on prend p = 0 (origine) => f(s) = 

et f{-e) = \-f(e). Autrement dit : un trou d’énergie s , par rapport à l’origine p = 0 , a 
une probabilité \ - f (e) qui est égale à f(~e) probabilité d’une particule d’énergie 
-s . De ce point de vue : ‘ ’un trou d’énergie s, est équivalent à une particule d’énergie 
relative - s” . 


exosup.com 


page facebook 


64 


d- Statistique de Bose-Einstein : 77 = -1 => f B _ E (£■) = 


1 


, P (s - m) 


e ' ' ’ ' - 1 

Propriétés de / (s) : 

a) Le potentiel chimique des bosons est négatif p< 0 pour £ = 0 (énergie minimum) 
alors que pour la distribution de Fermi-Dirac : p peut être > 0 ou < 0 . 

En effet, f(s) doit être >0 ^ - 1>0 Vf .=><?“ ^ > 1 <=> e < 1 <=>//< 0 . 

p)s>n. 

En effet, f(s) donne la répartition moyenne des particules dans les différents états ; elle 
n’ a de sens ici que pour : s > p . 

r) Le phénomène de condensation de Bose”. 

A T-OK , toutes les particules sont dans l’état fondamental d’énergie £■ = 0 . Dans le 
voisinage de T = 0K , une partie importante du gaz de bosons peut encore se trouver dans 
l’état fondamental. 



e- Etude des deux statistiques 
a ) Fluctuations : N- ^N À => N = Z 


et 


, . . d Ni , , . , . a 2 (N A \ \ + f À ' pour B - E 

Sx, Ni -1+ fl pour F -D 


Les fluctuations relatives ne s’annulent pas V le nombre total de bosons ; ceci 
résulte des échanges de particules avec le milieu extérieur. Cependant pour un gaz de 
fennions, les fluctuations a 2 (N À ) sont nulles pour les états certainement vides (f À - 0 ) et 
certainement occupés ( /, = 1 ). Pour les bosons ces fluctuations sont nulles pour les états 
certainement vides (f À = 0 ) uniquement. 

P ) Energie et entropie- Grandeurs thermodynamiques 

• £ = Z /V ^-C. = ZAC 1 


• Q. - -kT Log Z G , or LogZ G =- t]^Log(l ~r/f À ) => Q = kTpJ^^S ~ 1 fi) 

• L’entropie : 


s = - k T;[fx Lo zfi + 7 ( 1 - nfx ) L °s { 1 -nfx)] 


En effet: - D - pN] = f, - m 7 L og (l -;//,)- ///,] ; or: 

fi = no. => e o P s i = PE + L °S (1 - r/fi) - Logf À <p> 


piSi-fl) 

r u ' + rj 


£1 = M + kTLog (1 - rj fi)~ kTL °Sfi 
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5 = + - -V/i) - kTLogf À ) - kTijLog (l - *7 jf A ) - ///*] 

* A 

= èZ[^ T (A _ 7 )^og(l - » 7 /i) - kT f À Logf À ] = ~ kJ^[(f À - 77 )Log(\ - rjf x ) - f À Logf À ] 

A A A 

= ~ k Ts[f* L °Zf* + 7 A l ~ Vfx ) Lo s( 1 - 7/J] 

X 

On retrouve alors la relation de shannon pour l’entropie d’un gaz de Fermi avec f À 
probabilité d’occupation de l’état | i// x ) par une particule et (l - j) ) par un trou. 


Remarque : Comme en statistique de Maxwell-Boltzmann : E -> Z g (s x ) ; où g(e À ) est la 


À c. 


dégénérescence des niveaux s À . 

Y Approximation des grands volumes 

Dans la pratique les dimensions des systèmes étudiés sont macroscopiques. Les 
niveaux s x sont très serrés ; les sommes discrètes peuvent être approximées par des 

intégrales : 


Z^ Z#(Ci)-> \g(s)<D(s)ds 


Exemple : 


cD{s) = 


dW 

ds 


N= J g(s)Œ>(s)f(s)ds et E= J g(s)Œ>(s)f(s)sds . 


est la densité des niveaux d’énergie (nombre de niveaux dans un 


intervalle unité d’énergie); W est le nombre d’états d’une particule dans l’espace des 
phases O et s 0 est l’énergie de l’état fondamental. 


Application : Gaz de particules libres (fermions ou bosons). 

i) Densité d’état : Pour une particule, on a : 

W = Jùr= = Z \dp x dp y dp z |ù 3 ^ f 


h' 


h 


-J dp x dp y dp z = Z J 4 n p 2 dp 


1 

n / 2 ^ 2 \ 

Or, pour une particule libre : s = - — = — p x + p + p, . Donc : 

2m 2m ' ' 


4 /rp 3 V 4V n(2msY 2 V(2 ms)/ 2 


3(2 ntiy 2>(2nK) 


6 


®(*) = 


dW V 


ds 4 n“ 


(2m/ h 2 Y s 1/2 


Remarque : On démontre donc généralement que ® (t) ne dépend ni de la forme du puits, ni 
des conditions aux limites sur ses bords ; et ceci dans la limite des grands volumes. 

— N 

ü) Le nombre N (et concentration /î = — ) des particules ; on pose s 0 = 0 : 

V (2 tt0 3/2 


N= \ g(s)(D(s)ds= N = Jg(ff) 
g zm | j 


4 k 1 


4n 

cV 2 


.1/2 


^2 mf 

h 1 J J 


V n 2 


-ds 


l_ 

e v ’ + r/ 


ds 


+ TJ 


iii) Densité d’énergie : Le même calcul donne 


E 

u — — = 


8 


V 4 n 1 


2 m 


J fi(£-ju) 

0 <? 1 + 77 


ds 
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Si on change la variable s — » ps- x , alors : n= — 

4 K~ 


2 mkT Y 

{ 


1 


dx et 


+ 77 


4n~ 


2mkT 

Tr 


\V2 


a] 


dx . 


+ q 


Formules donnant le potentiel chimique p : n= n 0 (m)F^(a) et u - —n 0 (m)kT F y (a) 


ou : zi, 


{m)=g 


{ 2 


A 3 et A =h(3mkT) 2 = longueur d’onde associée à une particule 


v -> y 


de masse m et d’énergie cinétique — kT . 


F n («) 




r(n + 1); 


dx 


sont appelées moments d’ordre n de la distribution de 


Fermi ou encore intégrale de Fermi. 


iv) On peut aussi calculer , par le même procédé, toutes les autres propriétés du gaz 
de particules élémentaires de Fermi et de Bose ; en particulier le grand potentiel : 


O = r/kT^Log^ - tj f x ) - q kT^ g (s) Œ) (s) Log (l - qf(s))ds . Soit : 


fl = - rjkTg 


V 


4 n 2 


fêfï 

s^Logll+qe ^ = g 

V 

' 2m Y 

U J { 

\ ) 3 

4 rr 

[tr J 


J B(s-u) 

0 <? + d 


de . 


Or : Q = - PV ; donc : 


PV = — E 
3 


: relation qui est aussi satisfaite dans le cas de la 


distribution de Maxwell-Boltzmann. En effet, la valeur de Boltzmann de PV est 
— 3 

PV = NkT <=> E = — NkT . On retrouve ainsi la relation de Clapeyron. 
f- Limite classique 

Si «<s: 0 , les distributions statistiques de Fermi-Dirac et de Bose-Einstein 

V s k > 0 qui seront identiques à la 
N 


f\ deviennent : f À = f(s x ) -e ^ = e a Ph 


distribution de Maxwell-Boltzmann : 
N 


N*=fM = ~ e 

z 


-p^i 


avec la condition de 


i(a) . 


normalisation : — - e a 

z 

Conclusion : A la condition a« 0 om N^z (ou encore a (a) <s: l), les trois 
statistiques F-D, B-E et M-B sont identiques. Ceci veut dire que la limite classique est 
obtenue pour les trois statistiques quand «<sc 0 ou N <tcz (ou encore a (a) <tc l). 

On vérifie bien que ««: 0 correspond bien à la limite classique A<z d ou 


77 n Q — 


3mkT 

h 2 


(limite de validité de la mécanique classique). En effet : a<r 0 
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F„ ( a ) — e a => 77 «c ti 0 



avec n ü (m} = g — /î. 3 



et 


A. =h{3mkT^ 2 ) => ( « <tc 


Ainsi la condition a «: 0 est satisfaite quand la concentration n est faible et la température T 
est grande (les effets quantiques sont alors négligeables et on utilise la distribution de M-B). 


Définition : Un système de particules identiques et indépendantes est appelé gaz dégénéré si 
l’approximation classique n ’ est pas valable, c’est à dire : a> 0 ou 77 » n Q . Au contraire à 

la limite classique, il est appelé gaz non dégénéré, ou gaz parfait classique ( a «c 0) . 

Ne pas confondre la dégénérescence d’un gaz avec la dégénérescence des niveaux 
d’énergie U 

— N — (3 £ 

Remarques : (1) En statistique de M-B à la limite classique, on a : N . = — e À « 1 => Les complexions les 


plus probables C + sont celles telles que : N À = 0 ou 1 . En négligeant les probabilités des autres 
configurations, la fonction de partition du système est approximativement : 


fait, si les particules sont indiscernables, il faut diviser par N ! car les N ! complexions ne sont pas distinctes, 


f 



Fig. : La distribution 1 pour particules libres : 
Fermions (1); Bosons (2) et Maxwell-Boltzman (3) 


Z = I W ({ /V. }),f' • -x^ Â »... avec W ({ AQ}) = CSLr * c * ~ 0 ou 1. Soit donc : Z = N ! E e . En 

c + C, 


N 



N\ C + 


c. 


La grande fonction de partition en statistique de M-B à la limite classique est pour des particules 



N N N ! N N\ 


N\ 


On vérifie qu’elle est identique à la grande fonction de partition des statistique de F-D et B-E, pour a < C 0 . 
iii) Phénomène de condensation de Bose : 




• fl (a) <SC 1 pourr grande et n faible (limite classique). 

• Quand T\ (et tend vers 0) et n forte, a (a)/' et donc a («)\ et alors pour l’état 


fondamental (s 0 = 0 ) => N 0 = 



S. 
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Ainsi quand a(a)\ et devient <1=> N 0 <0 ce qui est impossible. Donc a (a) admet 1 comme valeur limite 
(critique) et alors N 0 devient oo . Mais 1V 0 =A=> \[ m A^a) = \et A(a) ^ 1 => A(a) = 1 + g avec 


g T ^ g > 0 . Ainsi : 


N 0 ->iV 


r-> o 


A 


. 1*0 ' 


0 + g) 


J £ i 


- 1 r->o 


e 


~P £ i 


->0 


Les états excités sont pratiquement vides ; la plus grande partie des particules se trouve ainsi condensée dans 
l’état fondamental : "condensation de Bose”. 


IV- Application aux électrons de conduction dans les métaux 


1°) Introduction 

Quelque 75 éléments chimiques sont des métaux (dans les conditions normales). 
Ils possèdent des propriétés communes : 

• Une faible résistivité : p < 10 “ 4 D.cm ( pazT ) 




Une grande conductivité thermique : k 

Vérifient la loi de Weidmann-Franz : k = LT a ( approximativement) où 


- 

3 l e) 


est le nombre de Lorentz et a - p 


• Un grand pouvoir réflecteur de la lumière 

• Une grande ductilité (ou malléabilités). 

Ces propriétés sont dues au type de liaisons des atomes du cristal : 

Liaison métallique : Les électrons de valence perdus par les atomes ionisés fonnent un gaz. 
Ces électrons ne sont pas localisés dans le cristal : ils ne sont pas liés seulement à quelques 
ions mais à l’ensemble du cristal. Inversement ces électrons sont responsables de la cohésion 
du cristal. 

Exemple : Na a un seul électron de valence : Na — » U 2 2s 2 2 p 6 3s 1 . Dans un cristal de 
Na les N électrons de valence de la couche 3s sont liés à l’ensemble du cristal ; alors que les 
10 électrons ( ls 2 2s 2 2p 6 ) internes d’un atome lui restent attachés ( électrons liés à l’atome). 


2°) Modèle du gaz d’électrons libres (Théorie de Sommerfeld) 

Certaines propriétés des métaux s’expliquent bien par le modèle dit d’électrons libres 
basé sur les hypothèses suivantes : 

i) Les N électrons de valence interagissent faiblement entre eux ; ils se meuvent 
indépendamment les uns des autres dans le potentiel moyen créé par les ions et les autres 
électrons. 

ii) Les variations spatiales de ce potentiel sont négligeables. Cependant la résultante 
des forces dues aux ions est grande (à la surface) et est dirigée à l’intérieur du cristal : les 
électrons de valence sont libres dans un puits de potentiel limité par la surface. 



iii) Les électrons de valence suivent la statistique de Fermi-Dirac. Ils sont responsables 
de la conductivité électrique : Ils sont appelés aussi " électrons de conduction” ; ils forment 
un gaz quantique parfait : gaz de Fermi. 
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a- Niveaux d’énergie des électrons libres 

• Le puits de potentiel est supposé cubique de côté L. 

• La fonction d’onde d’un électron libre (indépendante de t) est : 


2 r 

— ]= e ,k ' r d’énergie: s [k] = — — 

yjÿ b \ ) 2m 


(dépend du vecteur d’onde k ). 


• Les conditions aux limites du puits => les valeurs possibles de k et de s . Par 
exemple, dans le cas des conditions de périodicité on a: 

2 n 


ik-L i 

e = 1 =>i 


= — n-i ; i=x,y,z et n { e N 


( r* \2 n 2 2 2 2 2 

£i = (2k) n ; n = n r + n, + n . 


2m L 


b- Densité des niveaux : ©(e) = 2m/h 2 ) ~ f 1/2 . 


\ 3/2 


3°) Etude du gaz d’électrons libres 

Les électrons sont des fermions => N À =f (^ ) = 


1 




e ' v " +1 

a- A T = OK : V » les états I y/ n ) sont vides pour e n > ju 0 et occupés pour s n < ju 0 




f*0 

i)Le potentiel chimique /u 0 est détenniné par: N=2^Œ>[s) de — nombre total 


d’électrons avec: g = 2 = 25 + lest la dégénérescence de spin 


N 

( ju 0 x qq eV) n e = — est la concentration des électrons. 

A) 

ii) L’énergie du gaz : E 0 = 2 J ®(f) s ds soit : E 0 


Mo = 


2m 


(3* 2 n e ) 


V 


(2m> 


5k 


n 2 


2 ^ 
Mo 2 


- E 0 3 3 

£o = — = — Ll n - — k . 

N 5^° 5 F 


o v« y 

énergie moyenne par électron. 


T — A 

F ~ k 


est la température de Fermi qui est de l’ordre de 10 4 K . 


Cette énergie (T F ) est purement cinétique et elle est importante, contrairement au gaz parfait 
classique pour qui elle est nulle. Ceci résulte du principe d’exclusion de Pauli interdisant plus 
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N 

d’un électron sur chaque niveau et donc oblige les électrons à se répartir sur les — niveaux 

les plus bas à T = 0 K : la dégénérescence du gaz est totale. 

1 dLogZ a 1 ^ dLogZ G ôx- 
p SV ~ P x ex. ÔV 


iii) La pression du gaz : 


P = — 


avec x À = e a 




8LogZ c de. 
ôx x ÔV 


= -z* 


de 


2 e 


ôs ^ o 

, — Or, 2ZLL- _ puisque s, = 2n~ 

1 ÔV ÔV 3 V À mÛ 


X 1 et 


Na 


V =L\ Soit : 


P ———u 
3 


ou : u — — . 
V 


Conclusion : Cette pression s’exerce sur le réseau cristallin. Elle est considérable P = 10 5 atm . 

2 

L’équation d’état du gaz parfait quantique s’écrit PV = — £ 0 : indépendante de la température. 

b-: Le gaz est dégénéré : 0 < T <sc T, : Comme T F est très grande ; on est en situation 
normale. Le gaz est toujours fortement dégénéré. En effet, //varie peu avec T. Par 

conséquent: ju = — <=> a = P pi - — =— »1 : gaz fortement dégénéré. 7, est aussi 

P kT T 

appelé température de dégénérescence. 

00 f ( x) 

Pour 1 , F (oe) = J r _ a ^ dx est approchée par son développement au voisinage de 1’ 


de la façon suivante : u =x - a=>du = dx ; u\“ a quand x|” alors : 

. , f (u + a) ° r f (u + a) f(u+ a) 

F (a) = | — -du= | — '-du + ’-du . 

V/J e u +1 J e u + 1 J 

-a c ^ 1 -a c ^ 1 n 


u , 1 

e +1 


r j mt a /• 

Or, changeant u^-u: I — — -du = I 

6 l n 


\f(u+a) a çf(a-u ) 1 

■ 1 1 du ,. Sachant que : =1- 


— u .1 

e + 1 


—u u , i 

e +1 e + 1 


? fla-u) %f(u+a ) p , . } fla-u) }f(u+a ) 

=> Fia) = ’-du + ’- du= \ fia- u )du - ’-du + ’-du 

’ { e~" +1 { e u +1 { y ’ { e" +1 { e“ + \ 

a a 

Changement de variables : t -a- w=> dt= - du^> f / («- u ) du - J / (t ) dt et alors : 

0 0 

. , r °rf(u+a)-f(a-u) 

F (a) = I f{x)dx + I — — -du car a» 1 . Développement de Taylor de f 

o o e + l 

f(u+a) = f(a)+ uf'(a)+^-f”(a) + ^—f "'(«)+ 0(u 4 ) 

2 6 

f (a - u ) = /(«)- uf'(a) + y/"(a) - y/"'(«) + o(w 4 ) 


Donc : 


F(a) L 

f(x)dx + 2 f'(a) 

00 1 

f — - — du + -/"'(«) 
J 0 e “ + 1 3 V ’ 

oo 3 

f — ///< 

{«“ +1 


e" + 1 


du sont des intégrales de Lermi dont les valeurs sont tabulées et données par : 
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Z* 1 


A " ' A 


; x >0 


3 


r(^) vérifient : r(jc + l) = xr(x) avec : r = et et £( x ) sont ^ es 

2 4 

fonctions de Riemann telles que: <f(2) = — ; ^(4) = —; = 2.612; 4”(/^) = 1-341 . 


Ce qui donne finalement : 


’(«)=/ f(x)dx + ?-f'(a)+ ^-/"'(«) + 


360 


On en déduit : * La concentration des électrons : 


*La variation de l’énergie de Fermi avec T : ju = Mo 


1 

f 2m) 

3 

2 2 

K 2 

f kT ) 

2 

n. = — - 


M 2 

1 + 


+ ... 

In 2 

U 2 J 

8 

l M ) 



1-^1 

12 


r kT ^ 


K Mo j 


+ .. 


Cette variation est faible ; elle est de l’ordre de 5.10 3 pour T g [0, 300/6] et si ju Q - 5 eV . 



c- Pour des températures suffisamment élevées : T^>T F 
Alors m <£ o = 0 . Dans ce cas la distribution de F-D peut être remplacée par celle de M-B. Le 
gaz est non dégénéré (limite classique). 

(3 Ce) 


d- Chaleur spécifique des électrons libres : Seuls les électrons aux voisinages du 
niveau de Fermi, dans la bande , peuvent acquérir aisément de l’énergie (~ 1 ef) et peuvent 

alors effectuer des transitions des états occupés < ju) vers les états inoccupés (s À > ju) . Les 

autres électrons sont figés à cause du principe d’exclusion de Pauli (les transitions pour ces 
électrons sont improbables pour des processus mettant en jeu des petites énergies comme 

l’énergie thermique kT-^eV , à température normale). Donc le nombre effectif des 

3 T 

électrons ayant transité est N eff « N tot , soit N eff - 2 © ( /v ) kT =>N efr - — N — . 

2 Tp 



exosup.com 


page facebook 


72 


Q = 


f ds A 


, ôr 

V Jv 


E= 


V ( 2m 


2k- \ h 


_d_ 

' srr 

3 

\ 


f 3 

^, ot kT 


\ 


3 3 T 

<=> C v = —N,„, k — .Un Calcul exact donne: 


kTF V/2 ) F y( a h E o 


1 

Ul 

r j \ 

2 


3 

( TT~ T A 
71 1 

1 + 

12 

•J- F , 

+ ... 

=> 

C v - — N k 
v 2 

3 T 

\ J F J 


Conclusions : i) C v est environs 100 fois plus faible que pour un gaz parfait classique 



ii) C v varie linéairement avec la température T. 


Aux températures ordinaires, elle est masquée par le terme atomique. Par contre aux 
températures très basses (=0/f) , elle est plus importante que le terme atomique T 3 ) . On la 

mesure donc dans cette gamme de températures. 
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Chapitre 7 : 

ELEMENTS DE THEORI E DE BANDES D'ENERGI E DES 
ELECTRONS DANS LES SOLI DES CRI STALLI NS 

I-Introduction 

La théorie des bandes a été introduite pour répondre à la nécessité de décrire les niveaux 
d’énergie des électrons dans un réseau cristallin et aussi d’expliquer les différences entre les 
solides métalliques, isolants ou semiconducteurs. La théorie élémentaire de Sommerfeld étant 
incapable d’expliquer ces différences. 

La différence particulière entre ces différents matériaux est donnée par la résistivité 
électrique p très différente, allant de : 1(L 6 Q cm à 10 22 Q cm . 

Nous traiterons le cas simple d’un réseau à une dimension 

II- Description qualitative des bandes d’énergie 

Facteur d’échelle A (A = ca tel que c e '.H et a est le paramètre du réseau) : 

Toutes les longueurs caractéristiques d’un réseau sont étudiées en fonction de A . Initialement 
A » 1 , les atomes sont alors isolés, sans interaction mutuelle. Ensuite A diminue progressivement 

jusque à A = 1 . 

1°) Le cas : A » 1 : Les N atomes sont isolés et indépendants. 

a- Rappel : Orbitale = état propre de l’énergie d’un électron dans le potentiel central 
créé par le noyau et la distribution moyenne de (Z - l) autres électrons. Chaque orbitale 

\<t) n _ , m ) est fixée par n,l,m : nombres quantiques. m ) -» s n , z 71 avec n et /. 
Chaque s n ; est 2(2/ + l) = g fois dégénéré, car le potentiel ne dépend pas du spin (facteur 
2) et il est invariant par rotation (facteur 21 + 1). 

Les g états correspondant à une même énergie s n , forme une couche électronique : 

K, L, M, ..., ou en notation spectroscopique : \s ,2s2p, 3s... 

b- Ensemble de N atomes isolés et indépendants 
Une orbitale = | </> n ; m ^ où q = 1, , N désigne l’atome correspondant à cette orbitale. 

Les N orbitales : correspondent à la même y, ; => Les niveaux s nl sont 

maintenant gN fois dégénérés. 

2°) Le cas : A —> 1 : N atomes interagissent. 

Aux voisinages d’un atome, le potentiel global est : 

• La somme du potentiel central de cet atome et du potentiel dû aux atomes voisins 
(potentiel d’interaction). 

• Il n’est plus central. 

• Il n’est plus invariant par rotation. 

• Il a la périodicité du cristal. 

• Les électrons peuvent passer d’un atome à un autre par effet tunnel => Ils ne sont plus 
liés à un seul atome, mais à l’ensemble du réseau cristallin. 
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• L’interaction entre deux atomes voisins peut être considérée comme un couplage entre 
leurs orbitales. Les tennes de couplage dans le hamiltonien jf sont des termes non 
diagonaux: (</>„,,, m , q \ x \ avec q±q'. 

Ce couplage modifie les et lève leur dégénérescence d’ordre N (pour n,l,m 

donnés et pour une composante S_ du spin de l’électron donné). 

• Chaque \</>„ iL ) d’un atome isolé donne N niveaux distincts sur lesquels peuvent être 

placés 2N électrons (à cause de la dégénérescence de spin). Ils fonnent une sous- bande. 

• Les 21 + 1 sous-bandes d’une même couche correspondant aux mêmes n et l, 
constituent une bande d’énergie dite bande permise (fig.). 



Fig. : Structure de bandes pour un cristal (à paramètre a variable) de Sodium. 

o 

La ligne verticale en pointillé correspond à la valeur mesurée a =3.7 A 

• La largeur d’une telle bande dépend de l’interaction entre les atomes. Elle est très petite 
pour les niveaux d’énergie plus basse (interaction faible devant le potentiel créé par le 
noyau). Elle devient importante pour les niveaux d’énergie plus élevés ; c’est à dire pour 
les couches externes. 

• Les bandes correspondant à deux couches peuvent se chevaucher 

Exemple : s, p externes des métaux alcalino-terreux (voir figure pour le sodium Na). 

• Les bande permises sont séparées par des bandes interdites BI (intervalles vides : il n’y a 
pas de solution de l’équation de Schrôdinger dans les BI). Leur largeur est de l’ordre de 
quelques électron- Volts pour les couches externes. 

• La disposition des bandes permises externes et leur remplissage par les électrons sont 
responsables des propriétés électroniques des solides cristallins qui se traduisent par une 
classification en métaux, isolants et semiconducteurs. 

3°) Classification : métal - isolant - semiconducteur 
a- Métal : Un métal est caractérisé par : 

• Une bande permise incomplète à T = OK , dite bande de conduction (BC). 

• Le potentiel chimique p est situé vers le milieu de la BC. 
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Exemples : i- Na : Dans un cristal il y a N électrons dans la BC à T - OK . Elle est à moitié 
pleine (elle peut accepter 2 N électrons). 

ii- Mg : La couche 3s 2 => 2N électrons dans la BC à T = OK qui est donc 

complète. Mais 3s et 3 p se recouvrent et fonnent une seule bande fortement incomplète. 


b- Isolant : Un isolant est caractérisé par : 

• La dernière bande remplie est complète à 0 K . C’est la bande de valence (BV). La bande 
de conduction (BC) est vide. 

• Le potentiel chimique p est situé dans la bande interdite(BI), entre la BV et la BC. 

• L’intervalle A s entre la BV et la BC est grand ( As = 5 eV pour le diamant). Ce qui 
explique que la BC reste pratiquement vide à température ordinaire. 

• Les électrons de la BV ne peuvent être excités par un champ électrique. 

Seuls les rares électrons de la BC peuvent conduire l’électricité 

=> La résistivité électrique p est très grande (10 14 à 1 0 22 Qcm ). 

Cependant p\ quand T Z : de plus en plus d’électrons de BV peuvent être excités vers la 
bande de conduction . 



c- Semi-conducteur : Un semiconducteur est caractérisé par : 

• Isolant à T - OK avec As entre la BV et la BC est petite ( As <1 eV). 

Si Ge SnTe 
As(eV) 1.1 0.7 0.2 

• La partie inférieure de la BC est faiblement remplie ( n - 10 12 à 10 17 électrons / cm 3 ) => 
possibilité de conduction électrique par ces électrons. 

• Les électrons de la BV peuvent aussi être excités vers les états vides de la BC. Ceci est 
équivalent à l’excitation d’un trou. 

• Cette excitation de la bande de valence produit aussi un courant électrique sous l’action 
d’un champ électrique. La conductivité d’un semiconducteur résulte donc, au total, de la 
présence d’électrons dans la BC et de trous dans la BV : 

P„,é, < P sc Z 10- 2 à 10 9 Qcm) <p isolant . 
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• p peut être diminuée de moitié, contrôlable par purification poussée du cristal puis par 
ajout d’impuretés choisies (ces impuretés peuvent augmenter les électrons dans la BC ou 


les trous dans la BV). 


T Cf K 

C m 

I occupas J 


BAW3L' DK COWCVCTICW 


BAMDK DK VAiffACF 


III- Chaîne linéaire d’atomes 

1°) Modèle simple 

• Soit une chaîne très longue d’atomes espacés d’une distance a et ayant une seule 


orbitale \é) par atome d’énergie s 0 sans recouvrement : 


e 

'q \Yq' / — °qq' * 


On diagonalise le hamiltonien W dans le sous espace engendré par {)$,)} • Le coulage 


est limité aux voisins : Jf , = 

11 \ r q 


H 


H = 


?/ \ i 

s 0 -A 0 • 

-A s 0 -AO 

0 -A • • 

• 0 


K 


1+ 1 


-A pour q'=q± 1 et 


il °o 


£ n . 


0 


• • • 0 
• • • • £ 0 -A 

0 • • 0 - A £ 0 

L’équation aux valeurs propres : H \y/) = e\ y/) avec \y/) - 

q 

A(C y _, + C q + 1 ) = (e 0 -e)C q avec: q= 1 , 2 ,.... 


donne : 


qui est un système d’équations linéaires couplées en supposant la chaîne infinie 
• les solutions convenables : 


iqk a 


C 9 =* 


avec: ~—<k<—, correspondant à la 

a a 


relation de dispersion £(k ) pour l’énergie £ : £{k) = £ 0 - 2 A cos ka 


On obtient une bande d’énergie dont la largeur est : 4 A , proportionnelle au couplage. 



• Les fonctions d’onde dans l’espace vérifient : y/ k (x + pci) = e ,pka y/ k (x ) avec pe Z . 

La probabilité de présence est donc périodique. Ceci exprime qu’un électron est hé à 
l’ensemble de la chaîne et non à un seul atome particulier. 

• On en déduit que les fonctions d’onde sont des fonctions d’ondes de Bloch : 


Va {x) = e ikx u k 


(x) 
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où u k (x ) est périodique. u k (x ) vérifie aussi : H u k (x )= s u k (x ) : Ceci est équivalent à 
une transformation de jauge. 

• Conditions aux limites (périodiques) : Pour une chaîne linéaire de N atomes, on a : 
y/ k U N - l)u ) = (0 ) => Les valeurs possibles de k sont : 


k„ = 2n- 


ou 


L = (N — 1)a et n — — 


N -1 


N -1 


Il y a autant de niveaux dans la bande d’énergie pennise qu’il y a d’atomes dans la 
chaîne. On retrouve ce qui précède: chaque orbitale atomique donne N niveaux dans la 
bande permise correspondant à cette orbitale. 


Cas limites : 


k = 0 et k - ± — 
a 


Les états sont des ondes stationnaires : 

v , ( x )=Z(- 1 ) ? ^( x )- 


v 0 (*)=IXM et 

<7 ' a <7 

Ceci résulte de la réflexion par les atomes de toute onde progressive, qui satisfait alors à la 
condition de Bragg : Pour ces valeurs de k les ondes réfléchies par deux atomes sont en 
phase. Les ondes incidentes et réfléchies interfèrent pour donner une onde stationnaire. 

• Vitesse de groupe : y/(x , t) = J u k (x)e ,(A ' ~ ‘ :U ' l ' n,] D(k)dk : paquet d’ondes localisées 

ou superposition incohérente d’ondes de Bloch. La position du maximum du paquet est 

donnée par : — [kx- s(k)t jh + a k ) = 0 , où a k est la phase de D(k)u k (x) . Cette position se 
ûk 


déplace avec la vitesse du groupe V telle que: 


1 ôe(k) 2 Aa . , 

V, = = smfai 

h dk ti 


en négligeant la modulation due à celle de D ( k)u k (x) ( a k est indépendant de x). 

frk 

Cette V est * — : la vitesse de groupe d’un paquet d’ondes planes. 
m 


Dans les cas limites : k = 0 et ± — , L ç = 0: pas d’onde progressive. 

Un électron se déplace donc dans le cristal comme s’il était libre 
(approximativement), mais avec une vitesse différente. 


• Masse effective : F est la force extérieure ; son action donne une variation de 

ôs dk 

l’énergie: de - FV dt-ds = — dk^>F = h — . Par ailleurs : 

* dk dt 


dV g _ 1 d 2 s (k) dk 
dt ti dk 2 dt 

effective m e : 


L’équation de Newton : 




1 d 2 £(k) 

h 2 Ôk 2 


2Aa 2 

= — - — cos ka 
h 


définit alors une masse 


(i) m e est inversement proportionnelle au couplage A entre atomes. Par conséquent, 
les électrons fortement liés des couches internes ont une valeur de m e plus grande que 
pour les électrons des couches externes. 

(ii) m e > 0 dans la partie inférieure d’une bande permise ; alors qu’elle est <0 dans la 
partie supérieure de cette bande permise (lié à la concavité de s (k)). 

(iii) m e des électrons de conduction dans un métal est plus grande que m des électrons 
libres . 
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(iv) m e des électrons de conduction est ~ m des électrons libres dans un 
semiconducteur. 

• Densité des niveaux : Aux voisinages du minimum de la courbe de dispersion, 
l’énergie peut être approximée par une parabole : 


électrons libres, avec une masse m e et d’énergie {s{k) - s a ). La densité des niveaux est 


Conclusion : Aux voisinages d’un bord d’une bande permise, les électrons se 
comportent approximativement comme des électrons libres d’impulsion fik , et de 
masse effective m e ^ m , dans un puits de potentiel uniforme caractérisé par une 
origine s a . Ce sont des quasi-particules . Ce modèle sera utilisé pour les électrons de 
la BV et les trous d’un semiconducteur. 

2°) Généralisation : Les propriétés d’une chaîne linéaire se généralisent à un réseau cristallin 
à trois dimensions telles que : 

• Au nombre d’onde k^>k(k x , k y , k z ) : vecteur d’onde. 

• Chaque orbitale atomique donne naissance à une bande d’énergie permise 
comportant un nombre de niveaux égal au nombre de mailles élémentaires du cristal ( ce 
n’est pas toujours le nombre total d’atomes, car une maille peut contenir plus d’un 
atome). La relation de dispersion est une fonction multiple^,, ycj ; l’indice n numérote les 
bandes permises. 



alors : 






p = -ftk 



En plus l’énergie s n [kj dépend aussi de la direction du vecteur d’onde k . Ceci 

résulte de l’anisotropie du cristal. 

Dans l’espace du réseau réciproque (espace des vecteurs k ) les surfaces 
ne sont qu'approximativement des sphères. 
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Figure - Représentation de* court*?* d'énergie constante pour un 
réseau carré k deux dimensions. On voit que pour k petit mi obtient un 
carda (une «.phare à trois dimensions'!. 

• Les bornes d’une bande permise dépendent aussi de la de la direction de k . Il peut 
même arriver que deux bandes se recouvrent, deux niveaux de deux bandes successives 
peuvent avoir la même énergie mais pour deux vecteurs d’onde k différents. 

• V g et m e dépendent aussi de la de la direction de k . 

y g =|v*f (*) 

- dV ** - 1 17 

1 _ I d-£ (k) dt (w, ) iv y 

( l %) xy h 2 8k x dk y 


• L’énergie aux voisinages d’un extremum est donnée par : 


K k v . 

1 ôk,ôk„ 


y 


Elle dépend de la direction du vecteur d’ondes k . 



III- Fig. : Relation de dispersion pour deux directions du vecteur d’onde 

IV- Les semiconducteurs 

1°) Modèle simple 

A T = OK , le cristal pur comporte une BV complètement pleine et une BC 
complètement vide. La largeur de la BI est as ~ 1 eV . 
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• BC et BV sont isotropes. On suppose aussi que BV admet un maximum et BC un 
minimum. 

• A T * 0K , la conduction électrique résulte de la présence d’électrons dans BC et de 
trous dans BV. 

• Domaine extrinsèque ( à basse température) : la conductivité est dominée par la 
présence d’impuretés. 

• Domaine intrinsèque ( à haute température) : les impuretés jouent un rôle négligeable. 



• Dans le domaine extrinsèque, les électrons de conduction dans BC sont fournis par des 
impuretés(donneurs) qui possèdent un excèdent d’électrons par rapport aux atomes normaux 
du matériau : on a alors un semiconducteur de type n . Les trous de la BV proviennent de la 
capture d’électrons de cette bande par les impuretés (accepteurs) présentant un déficit 
d’électrons: le semiconducteur est alors de type p. 

• Dans le domaine intrinsèque, les électrons de conduction et les trous proviennent de 
l’excitation d’électrons de la BV vers la BC. 

• Dans les deux cas, les électrons de conduction peuplent les niveaux les plus bas de la 
BC, au voisinage du minimum (k = 0 ). De même, les trous sont au voisinage du sommet de la 
bande de valence BV. 

• L’énergie s (k j est alors donnée aux voisinages de ces extrema par : 


As + 


(*)- 


n 2 k 2 

2m e 

h 2 k 2 
2 m „ 


-> BC 


BV 


2°) Propriétés des trous 

L’absence d’un électron dans la BV est équivalent à la présence d’un trou dans la BV. 
Ainsi, au trou on associe des propriétés qui sont celles d’une particule : impulsion, énergie, 
vitesse, masse et charge électrique. 

■ La transition d’un électron de la BV d’énergie s e i [ k e j au niveau s e (à; (, ) de la BC est 

équivalente à la création d’un électron dans BC et d’un trou dans BV ; c’est à dire à la 
création d’une paire ” électron- trou”. 

■ Inversement, la transition d’un électron de la BC vers un niveau inoccupé de la BV est 
équivalente à l’annihilation d’une paire électron-trou” . 
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a- Vecteur d’onde d’un trou k P : Création de ” électron- trou” => absorption de 
l’impulsion : 


hk - h(ke ' - ke ) 


. /r , 7 \ ^ 

k p — k e 

= hyke + k P \ 



b- Energie d’un trou 

L’énergie absorbée au cours d’une création d’une paire 1 ’ électron-trou” sera (par 

, ,,, . . . E e {ke') - S e {k e ) 

conservation de 1 energie) : W = < _ _ <=> s Ak P ) = - s e (k e ) 


£ e (ke ') + £ p (kp) 

Avec l’origine au sommet de la BV => s p (k p )> 0 =>La surface s (k p ) est symétrique par 
rapport à l’origine de s e (k e ) : 



c- Vitesse : 


d- Masse : 



<=> 


V = V 

P e 


m 


1 S 2 s n 


1 d 2 s„ 


ri dk; 


<=> 


m p = - m e 


h 2 dk 2 

Dans la partie supérieure de la BV, m e < 0 donc m p > 0 . En général, m e dans BC et m p dans 
BV ne sont pas égales. Pour un semiconducteur, elle sont toutes deux plus petites que la 

m e p 

masse m d’un électron dans le vide : — 0.1 . En effet, une faible As de BI => une forte 


m 


largeur de BC et de BV => fort couplage => m e faible. 


e- Charge électrique : - q^s + v p AB) <=>- fi-j- = q p (£ + V e As) . Or. 


dt 


= -\e\(s + Ve A donc : 


dt 


<1p = + 
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Dans un champ électromagnétique un trou se comporte comme une particule de charge 
positive. En particulier, £ déplace le trou en sens inverse de l’électron, mais provoque un 
courant électrique de même sens ( j = o £ ; axe 2 ). 

Conclusion : Les trous se comportent comme des électrons de charge et d’énergie 
opposées. Ce sont les antiparticules des électrons. Cependant, il s’agit d’électron et de 
trou liés dans un cristal. Tout se passe comme si ils n’interagissent pas avec les 
noyaux atomiques. Ce sont des quasiparticules. 

Enfin, les trous comme les électrons suivent la statistique de Fermi-Ditrac : 


fp( £ p)= l - fe{ £ e) 


1_ 

/f +£ p ) 


+ 1 


3°) Domaine intrinsèque 

a- Concentration des porteurs n (électrons de conduction) et /;( trous) 


00 

n= 2 J cD e (s)f e (s)ds 

A E 

CO 

P = 2 l CD p{ £ )f P { £ ) d£ 
0 

1 f o A 3 / 2 

/ x 1 2777 / . \ 1/2 

J (£ " A£) 

1 f o A 3 / 2 

/ \ 1 AÎYl i/o 

77 = n 0 F f (« - A) 

p = p {) F yi ( - a) ; a= fp et A - j3As 

n 0 = n(m e ) 

Po =n(m p ) 


où: 77(777)= 2 


f 2n mkT^' 2 


et m = m e , m . Les concentrations des porteurs s’écrivent en 


fonction de l’intégrale de Fermi : F l/2 (cc) = f 


PI 2 


e x ~ a + 1 


- dx . 


OÙ: 2DF = 2 j©0)/0) 


de 



b- Position du niveau de Fermi : Considérons le cas où les électrons de BC et les 
trous de BV constituent deux gaz non dégénérés alors: n <sc ? \ et p <sc p Q , soit : //:» kT 
et Ae -//» kT donnant n - n 0 e ~ p ^ E ~ ^ et p - p 0 e~ flM . 


L’énergie de Fermi est déterminée par : n = p . D’où : 


As 

3 kT T m 

P = — + 

A LOg 

2 

4 m e 
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Conclusions : 

• Si les masses effectives m e et m p sont égales => p se trouve au milieu de la BI, pour 








toute température T. 

Sinon (cad si m e ^ m p ) => p varie linéairement avec T. 

Quand T /*, p se rapproche de l’une des deux bandes. Un des deux gaz devient 
dégénéré à partir cl’une certaine température (As - p< k ou p< kT). 

Les concentrations des deux types sont égales : 


n — p = n Q 


(>n„ m p ) 


1/2 




augmente donc avec T. 


c- Conductivité intrinsèque : a . 

• cr est proportionnelle à la concentration des porteurs et à leur mobilité 


u = 


V 


S 


telle que : 


CT = 77 ll„ + 



• La mobilité varie peu ; elle diminue légèrement avec T. Elle est limitée par les 
interactions avec les défauts cristallins et les vibrations des noyaux atomiques. 

• Exemple à T = 3 00 K , on a les valeurs de u (unités : cm 2 /v-S ) : 




Ge 4500 3500 
Si 1300 500 


cr varie donc comme la concentration : 


cr oc e 


-flAe/2 


( As énergie du gap : As ~ leV ) 

Pour T < 300 K , on constate une déviation importante à la loi cr ~ e~ pul2 ( voir fig.). 
C’est le domaine extrinsèque : les impuretés ont un rôle qui devient important 



4°) Domaine extrinsèque 

Dopage : Un cristal de Ge ou de Si peut être purifié à un très haut degré (concentration 

limite d’impuretés : 10 12 cm' 3 ). On peut ensuite introduire une impureté de type donné en 
quantité contrôlable, par exemple, de l’arsenic (As) ou du bore (B) dans un cristal de Ge. 
Alors il y a substitution d’un atome de Ge par un atome As ou B présentant une différence de 
charge du noyau et aussi du nombre d’électrons de valence avec le Ge. 

• As est pentavalent ; il peut être considéré comme l’ion As + tétravalent comme Ge. A 

cet AC est lié un 5 eme électron de valence. Ce dernier sera facilement excité vers la 
BC du cristal. L ’ As est un donneur. 
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• B est trivalent, il peut être considéré comme un ion B tétravalent comme Ge auquel 
on enlevé un 4 eme électron de valence. Ce qui est équivalent à envisager B comme B 
auquel est lié un trou. Le 4 eme électron manquant pourra être facilement fourni par la 
BV. Le bore est un accepteur . 

a- Niveaux donneurs : 

i) Pour les électrons de conduction, un réseau de Ge avec un As + dans un site donné est 
équivalent à un puits de potentiel uniforme (ayant les dimensions du réseau), plus une charge 
+ \e\ localisée sur le site. Cette charge positive est donc équivalente à un noyau d’hydrogène 

auquel est lié le 5 eme électron de valence pour former un atome hydrogénoïde neutre Le 
niveau fondamental de cet atome hydrogénoïde constitue un niveau supplémentaire qui 
s’ajoute aux niveaux d’énergie du puits de potentiel constitué par le cristal. C’est un niveau 
donneur. Il est très localisé dans l’espace (contrairement aux niveaux du puits de potentiel) : 
un électron sur un tel niveau est lié à un ion As + .Sa fonction d onde d espace est limitée au 
voisinage du site. Son extension est en général beaucoup plus petite que la distance entre 
impuretés : il n’y a pas de recouvrements entre niveaux donneurs voisins. 


ii) Calcul de l’énergie d’un niveau donneur : 


Niveaux de l’atome d’hydrogène : 


, , 1 me \2>.6eV 

2/L~ (4 nsjdf = n 

L atome hyperboloïde à une fonction d’onde d’espace d’extension dans le cristal qui 
recouvre environ un millier d’atomes. Il faut donc remplacer m par m e et multiplier £ o 
par la constante diélectrique £ du cristal. 

L’énergie d’ionisation de cet atome est : 


-13.6 

m s 


en eV. 


Pour Ge: 


m. 


-= 0.2 ; e =15.6 


A e j - 0.013 eV . 


m 





Le rayon de Bohr a caractérise l’extension dans l’espace de ce niveau : 

4 7is n 1i 2 m m 

a = —e — - — 

me m e m e 




L : longueur du cristal. 
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b- Niveaux accepteurs : 

Même raisonnement pour (f? 

résultats : 

• L’énergie d’ionisation: 


+ 1 trou) qui est un antiatome hydrogénoïde conduit aux 


2 9 

m £ 



soit : A e a - 0.01 eV 


Remarque : A s d ou Af o sont de l’ordre de l’énergie thermique à température ordinaire ( 1/40 eV ), les 
électrons ou les trous sur de tels niveaux sont donc très facilement excités vers la BC ou la BV. Ceci 
explique la conductivité extrinsèque des semiconducteurs à basse température. 


c- Statistique des électrons sur les niveaux donneurs : 

Etant donné que les niveaux donneurs sont très localisés dans l’espace, il ne peut y avoir plus 
d’un électron sur chacun de ces niveaux (sinon, les électrons, même de spin opposés, seraient 
très proches et la répulsion Coulombienne devient plus forte que l’énergie s d d’ionisation). 

Il n’y a donc qu’un seul électron avec deux orientations de spin possibles. 

La statistique est par conséquent différente de celle de Fermi-Dirac habituelle ( quand on peut 
mettre deux électrons de spin opposés sur un même niveau). 

f \ 

• La grande fonction de partition d’un niveau donneur d’énergie s d : z G - ^ Y\ x x l 

PUV a 

avec 2 = 1,2 (il y a deux niveaux de même énergie avec spin T ou i pour un électron). 
Par conséquent : {N À } = {0,0} ou {0 , 1} ou {l , 0} et z G = x d x d + x° d x d + x d x d 

-" a - pSi ; soit: 


z G = 1 + 2x d 


Pour N,, atomes donneurs : 


Z i' 

G ~ Z G 


( 1+2 *,)\ 


Le nombre moyen de niveaux donneurs occupés est : 


dLogZ c 

ôa 


= N„ 


dLogz c dx d 
ôXj d a 


soit : 


AL 


+ î 


Probabilité d’occupation d’un niveau donneur : 


fi = f ( £ d ) = — r 


1 

,p(td - MÏ 


+ 1 


d- Statistique des trous sur les niveaux accepteurs 

• Le nombre moyen de niveaux accepteurs occupés par un trou est : 

„ = ^ 

" + 1 

où N a est le nombre d’atomes accepteurs et s a est 1’ énergie d’un niveau accepteur. 

• Probabilité d’occupation d’un niveau accepteur par un électron : 
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/ =1 — — 
Ja N„ 


1 


2 e 


P(£i - e) 


+ 1 


• La condition d’équilibre entre les concentrations n et p porteurs (électrons de 
conduction et trous de la bande de valence) est : 


n + N a -n a = p + N d - n d 


soit : 


2j <D e {s)f e {s)ds + N a fa = 2]<D p {e)f p { £ )de + N d (l- f d ) 


Cette équation doit être résolue (en général numériquement) pour détenniner l’énergie de 
Fermi . On en déduit ensuite les concentrations n et p des porteurs. 


5°) Semiconducteurs de type n 

Prenons le cas simple de semiconducteur dopé par une impureté de type n : 

• N a - 0 et N d * 0 . 

• As est assez large et T assez basse pour que p « 0 (dans BV). Les électrons de 
BC sont porteurs majoritaires et les trous porteurs minoritaires. 

• La condition d’équilibre devient (condition de neutralité) : n - N d - n d => 

les niveaux donneurs jouent le rôle de la BC dans la région extrinsèque. A^est 
très petit devant le nombre de niveaux de la BC=> Le nombre d’électrons excités 
ne peut dépasser N d =^> il y adonc saturation. 












AT = 0 K 



Quand T /* => n/' et atteint la saturation pour n = N d , puis reste constant dans 
la région de saturation : C’est la région intéressante dans l’utilisation, car la 
conductivité y dépend peu de la température. 

A température plus élevée commence la région intrinsèque où l’excitation 
d’électrons de BV domine. 

La condition d’équilibre est résolue, en prenant l’énergie à s d - 0 : 

N 

n = n 0 C /2 {a -S) et n = N d - n d = 


avec : a- p p et 5 = (3 As d . 

La relation correspondant aux points d’intersection de : 

y = ^ F in( a -8) 

< d 

y = {le “ + l) 
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69V K 



a- Région de saturation : 

Tl 

• Elle correspond à a <tc 0 , soit : 1 et n , - 0 . 

N d 

• On a aussi : a <s: S c ad p<s. As d . /u est très au dessous des niveaux donneurs. Ceci 
exprime encore qu’il y a ionisation totale des atomes donneurs. 

• Le gaz d’électrons de la BC est bien non dégénéré. 






t A 


M 




/*■’ — 


b- Région d’ionisation : 

N n 

Elle correspond à a » 0 , soit : n <tc IV , , or : n - N. - n, = => «c 1 

d d 2e a + \ N d 

Le gaz d’électrons de conduction si (a -£)«:- 1 , c.a.d. pour kT <sc As d - /u . 

Le niveau de Lenni se trouve entre la BC et les niveaux donneurs, et la température est 
très basse. 



i 


• La variation de la concentration n des porteurs, comme dans la région intrinsèque, obéit à 
une loi d’action de masse(admise) : n -yj N d a(T) e~^ A£d/2 où a (T ) est une fonction 
douce de T. 

6°) Jonction : p - n 

Une moitié d’un cristal semiconducteur est dopé avec une impureté de type n (As par 
exemple), l’autre moitié avec une impureté de type p (B par exemple). 

A l’équilibre thermique, il se produit un saut de potentiel à l’interface entre les deux régions. 
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A l’interface, les concentrations n et p présentent un gradient élevé. Des électrons passent de 
la région n vers la région p, où ils se recombinent avec des trous. Inversement, des trous 
passent de la région p vers la région n, où ils se recombinent avec des électrons. 

Aux voisinages de l’interface, la charge positive des ions As + n’est plus compensée par celle 
des électrons de conduction 

=> il apparaît de ce côté de l’interface une charge d’espace positive. 

De même, le charge négative des ions B n’est plus compensée par celle des trous 
=> il apparaît de ce côté de l’interface une charge d’espace négative. 

Ces charges d’espace produisent un champ électrique £ dirigé vers la région p , donc un 
potentiel électrique plus élevé dans la région n. 



exosup.com 


page facebook 


89 


Chapitre 8 : 

GAZ DE BOSONS EN NOMBRE NON CONSERVE 
STATI STI QUE DE PLANCK 

-I- Introduction 

• Les statistiques canonique et grand-canonique reposent sur le principe de la conservation 
du nombre de particules. Ce nombre est conservé pour l’ensemble formé par le système 
étudié et le thermostat C avec lequel il interagit. 

• Certains systèmes possèdent des modes de vibration qui sont ceux d’oscillateurs 
harmoniques. Les quanta d’énergie de ces modes peuvent être considérés comme des 
particules indépendantes. Leur nombre n’est pas conservé, car les échanges d’énergie 
créent ou annihilent des quanta. 

Exemple : Les photons, ou quanta du champ électromagnétique, sont des particules réelles. 
Les phonons, ou quanta des oscillations des atomes dans un cristal, sont considérés comme 
des quasi-particules. Ces particules fonnent un gaz de bosons en nombre non conservé. 

Remarque : La conservation du nombre de particules massives (m ^ 0) résulte de la conservation de la 

masse. Il n’en est rien en Relativité : à cause de l’équivalence masse-énergie, des particules peuvent être 
créées en fournissant de l’énergie ou annihilées en libérant de l’énergie. Ces particules doivent être des 
bosons, pour conserver le moment cinétique ( les fermions doivent être créés par paires). 


-II- Statistique de Planck 


• La seule quantité conservée est l’énergie => Seuls les échanges d’énergie avec un 
thermostat sont à considérer => Il faut utiliser la distribution canonique. 


• De manière équivalente, on peut utiliser une statistique grand-canonique avec une énergie 
de Fermi nulle ( a = 0 ). 

Justification : Prenons pour origine l’énergie du vide (gaz de zéro particules), l’énergie de 
Fermi p du gaz doit être telle que : p < 0 . Pour interpréter la non- conservation des 
particules on peut considérer le thermostat C comme un réservoir de particules en nombre 
infini. Ceci implique que ju c > 0 (la condition de dégénérescence du gaz). Or condition 


d’équilibre 


M = M c = 0 


1°) Fonction de distribution de Planck : 

• Fes opérateurs H : hamiltonien du gaz de bosons (en nombre non conservé) et N sont 

tels que: N =^N À et H = JXo | • 

2 2 

J -P'5L N i e i\v'x){v , x\ 

• F ’ opérateur densité : p = — e l 

• Fa statistique de Planck est obtenue en reprenant les formules de la statistique de Bose- 
Einstein avec la valeur a - 0 . 

• Fa fonction de partition : Z = ^ - p £ 
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• La distribution de Planck : 


AC = f { £ z ) = 


1 


, Pe 


e - 1 

• Valeurs moyennes : £ = '^, £ xf{ £ x) et A r = ^f{ £ ^) 

A A 

2°) Limite classique : s À 2 > kT 

La limite où les particules sont parfaitement localisées (longueur d’onde associée 
négligeable) et où elles suivent la statistique de Maxwell-Boltzmann ne peut être atteinte 
au voisinage du niveau s = 0 , car la concentration y est très grande (voir figure). Autrement 
dit, la condition a «: 0 pour que la limite classique soit atteinte n’est pas satisfaite, puisque 
ici : a = 0 . 

Toutefois cette limite est atteinte pour des énergies élevées telles que s À » kT . Cette limite 
partielle concerne une partie du gaz d’autant plus grande que T est basse. 



-III- Assemblée d’oscillateurs harmoniques 


1°) Propriétés : 

Considérons un système possédant M modes de vibrations non couplés. 


H = 


ou 


^ = + si). 


â=i 2 

a > ; l : pulsation du mode À et P. et Q À sont les coordonnées conjuguées réduites : 
impulsion et position. 


L’énergie totale est : 


E = £ (AC + 1/2) tuo. 


Les quanta Tuo x en nombre N À pour le mode A constituent un gaz de bosons (en nombre 
non conservé) indépendants. 


La fonction de partition est : 

1 




, -p*«>xl 2 


1 - e 


- pt“°i 


Na = 


, Pt“°i 


- 1 


et 


U=E = X 


, Pf«°A 


-1 


+ 1/2 


fuo, 


L’énergie fondamentale : E 0 = 1/2 est obtenue à T = 0K (énergie du point zéro). A T 

Al 

très élevée E est proportionnelle à T : E - MkT . 


exosup.com 


page facebook 


91 



2°) Approximation classique des grands nombres quantiques : 

Si un mode comporte un grand nombre A, de quanta, son énergie A, tuo i est très 
supérieure à l’espacement lia> i entre niveaux. Ce mode est donc celui d’un oscillateur 
classique dont l’énergie est kl (deux degrés de liberté). 

Ceci se produit pour : kT » hco i car N : - kT / tu» j est grand. 

Résumé : L’approximation classique est maintenant valable dans les deux cas suivants : 

i) tw) l » kT -> Les quanta se comportent comme des particules localisées (la 
longueur d’onde associée est négligeable). 

ii) t cOj « kT Les oscillateurs se comportent comme des oscillateurs 
classiques. 

Remarque : Cette limite est différente de la précédente. 


IV- Corps noirs 


1°) Photons : 

• Une onde électromagnétique, d’amplitude : A^r, t ) = A 0 je d 3 k 






se propage avec la vitesse c = co/k (vitesse de phase) . 

Le champ électromagnétique de cette onde est quantifié, on le décompose en une 
assemblée d’oscillateurs harmoniques. Le photon est le quantum d’énergie : s = ho = bco 

et de quantité de mouvement : p =hk ou p = hu/c . 

Le spectre des modes d’un champ électromagnétique dans une enceinte cubique, de côté 


L, est donné par: 



avec : n 2 = n x 2 + n 2 + n z 2 et n i e N . 


Ce spectre est obtenu en imposant la condition aux limites : ” Amplitude nulle sur les 


parois de V enceinte” . Le spectre en fréquences est : 



Chacun de ces modes est double, car il y a deux états de polarisation indépendants : 
polarisation rectiligne suivant deux directions transverses orthogonales entre elles, ou de 
manière équivalente : polarisation circulaire droite ou gauche. En effet : 

£(r,t) = q e ,iyk " où q est le vecteur de polarisation. La transversalité de £ 


<=> V • £ — 0 ^>£*/7 = 0^>3 seulement 2 directions de polarisation orthogonales à k . 


2°) Equilibre du champ électromagnétique dans une enceinte : 

• Le champ électromagnétique crée un courant dans les parois de l’enceinte. Il s’y produit 
une dissipation de l’énergie. Il y a donc absorption du champ. Inversement, l’agitation 
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thermique des électrons des parois y fait régner un courant électrique désordonné. Ce 
courant émet une onde électromagnétique dans l’enceinte. 

• Les photons du champ de l’enceinte constituent un gaz de bosons en équilibre thermique 
avec les parois qui jouent le rôle du thermostat. 


3°) Densité de modes : 

Pour chaque état de polarisation le nombre W de modes de fréquences inférieures à u 

2 L 

est égal au volume d’un quadrant de la sphère de rayon : n 0 = — u. Soit donc : 

c 


f , 


W - — 


— k n. 


V- 


4 

3 


= -/z-L 3 ^- 


La densité des modes : 


, V dw A . vf 

<d(o] = 4nL — 

w du c 3 


<=> 


<D(o) = 4 nv~ 


V 


4°) Densité spectrale d’énergie : u (v) 

Définition : La densité spectrale d’énergie u(v) est la densité d’énergie par unité de 


{v) = -(D(v)f(v)hv. 


volume et par unité de bande de fréquences : 

Le facteur 2 est dû à l’existence de deux polarisations transverses. Soit : 

(Loi de Planck). 


‘00 = 


8tt/2 


3 Bhu 

c e p 


- 1 



• Cette loi donne la répartition de l’énergie suivant la fréquence, c’est à dire, dans les 
différents modes. Cette répartition ne dépend que de la température (spectre thermique). 

kT 

• La position du maximum u m varie linéairement avec T : u m =2.822 — 

h 

(Loi du déplacement de Wien). 

• Pour : h v » h u m (limite classique) on a la formule de Wien : u (v) e ~ phu (M-B) . 

c 

• L’approximation classique des grands nombres quantiques : hv<£ h v m , donne la formule 

de Rayleigh- Jones : u(v) ~^^-kT u 2 . 

c 

Rq : La constante de Planck n’apparaît pas dans cette dernière formule purement classique. 
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5°) Puissance rayonnée par un trou: (mesure de u ( v)dv ) 

La mesure de la puissance rayonnée par une ouverture dS d’un corps noir et captée 
selon un angle solide dÇl(= d 2 co) pennet la mesure de u(v)du (vérification de la loi de 
Planck). 

• Les rayonnements qui peuvent sortir sont ceuxcontenus dans le volume : c dt cos 9 dS 

• La puissance émise est : 


d(p- u[v)du x c cos 9 dSx dQ./4n 

densité d'énergie volume / 1 angle solide relatif 

par unité de volume 



(Pu* 


• La radiance R est la puissance totale rayonnée par un trou de surface unité ( dS = l) est 
R = c|||m(v)i/(;cos0 dÇl/4n <=> R - —^ u ( v)dv = — \ ^ J - — 


/7 3 c 2 J 0 e * - 1 


L’intégrale : J 


r x dx n 


- 

R = a T 4 

le* - 1 15 



: loi de Stefan-boltzmann (1879). 


cr = 


x 4 k 4 
60 h 3 c 2 


= 5.67 10 12 Watt/ cm 1 {°K)' constante de Stefan-Boltzmann. 


La radiance très faible à basse température, augmente très vite avec T ; elle passe de 
0.046 w/cm 2 pour T = 300/6 à 460 W / cm 1 pour T = 3000K. Sa mesure permet de 
déterminer la température de l’enceinte. 

Remarques : i) Une enceinte percée d’un trou, de petite dimension est appelé ‘'corps noir”. Une onde 
électromagnétique entrant par le trou a une probabilité négligeable à en sortir. Elle sera absorbée lors des 
multiples réflexions sur les parois. Le trou est totalement absorbant. C’est le cas aussi d’un objet noir mat pour la 
lumière visible. 



ii) L’Univers considéré comme un corps noir est rempli par un rayonnement ‘'fossile” correspondant à T-3K 
(ondes centimétriques ; découvert en 1965) et qui nous renseigne sur l’époque où la matière était ionisée : 
T >3000 K (0.25eV). 11 y avait alors équilibre thermique entre rayonnement et matière ( à travers les 
interactions rayonnements-électrons; ces derniers sont libres) => Le rayonnement ne peut se propager : opacité 
de l’Univers à tout rayonnement électromagnétique. Cette interaction empêche par ailleurs la contraction 
(pression du rayonnement). A T = 3000 AT cessation de l’ionisation : les électrons se regroupent avec les protons 
pour former des atomes de l’hydrogène => galaxies, étoiles, . . . Les rayonnements sont restés après la formation 
des galaxies et des étoiles avec des fréquences réduites à cause de l’effet Doppler. 

Refroidissement de l’Univers : T est inversement proportionnelle à la dimension de l’Univers, car VT 3 =cste en 
admettant une dilatation adiabatique de l’Univers. 

On peut aussi connaître l’âge de l’Univers (en connaissant la vitesse de son expansion). 
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6°) Pression de radiation: 

Un champ électromagnétique a une densité de quantité de mouvement: (l/c 2 ) £ AH . 

Ceci est équivalent à dire que cette densité peut être attribuée aux photons, chacun d’eux 
ayant la quantité de mouvement : hv/c . 

L’absorption, l’émission ou la réflexion d’une onde électromagnétique par une paroi 
se traduit par une variation de la quantité de mouvement, donc par une force exercée sur la 
paroi ; c’est la pression de radiation. 

Comme pour un gaz d’électrons de conduction, la pression d’un gaz de photons est : 

p=- IX 


dV 


. Pour un gaz de particules ultrarelativistes : s À = h o, = hA 


— - hA — 

2 L 2V 1/3 


ds , 


dV 


1 S. 1 V -1 E 

~ — => P= / N, S , = — <=> 

3 V 3V z f 3V 


P=“ 

3 


Le calcul de la pression ne dépend pas de la distribution statistique des photons dans 
les différents modes possibles. On peut aussi calculer P en traitant tous les photons comme 
des particules localisées, à la limite classique : hvts> kT . Le calcul est semblable à celui d’un 
gaz parfait classique. La seule différence résulte de la relation existant entre l’énergie s et 

f p - s/ c pour un photon 

l’impulsion p : < 

I p = y/ 2ms pour une particule non relativiste 

Le résultat est le même pour l’équation d’état: P = u/3 . 


V- Les phonons 


1°) Introduction qualitative 

Les atomes d’un réseau cristallin peuvent effectuer des petits déplacements autour des 
sites qu’ils occupent. Chaque site constitue une position d’équilibre stable pour un atome. 
L’énergie potentielle U de l’atome est minima en ce point. Au voisinage, U peut être 
approximée par une forme quadratique du déplacement r(x, y, z): U =U 0 + ax 2 + by : + cz 

avec un choix convenable des axes. Ceci correspond à un oscillateur harmonique à 3- 
dimensions. Les N atomes du réseau constituent un ensemble de 3 N oscillateurs harmoniques 
à une dimension. Cependant, ces oscillateurs ne sont pas indépendants. Ils sont couplés entre 
eux. Cette situation est analogue à celle des électrons d’un solide étudié précédemment. Ces 
3N oscillateurs couplés sont caractérisés par (3 N - 5) modes de vibration, en retranchant les 
trois translations et les deux rotations de l’ensemble du cristal, soit approximativement 3 N 
modes pour N grand. Ils sont analogues aux N niveaux d’une bande permise des électrons 
d’un solide. 

Le cas le plus simple est celui d’un réseau cubique avec un atome par maille. Les 3N 
modes se décomposent en N modes longitudinaux, pour lesquels les atomes oscillent dans la 
direction de propagation des ondes de vibration et 2N modes transversaux pour lesquels ces 
oscillations sont transverses. Ces ondes de vibration sont des ondes acoustiques se propageant 
avec des vitesses différentes c, et c t suivant la direction de vibration. 

Dans le cas général d’un réseau comportant n atomes par maille et p mailles, on 
distingue 3p modes acoustiques pour lesquels les n atomes d’une maille oscillent en phase et 
3(/7 - l) p modes optiques, pour lesquels ils oscillent en opposition de phase. Ces deux types 

de modes ont des relations de dispersion co(k^ différentes : 
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/-\ | c k — » modes acoustiques 

V J k*o ± ck modes optiques 

Les 3N modes de vibration peuvent être considérés comme les modes de 3 N 
oscillateurs indépendants entre eux. Les quanta de ces 3 N modes sont appelés : phonons. Ils 
constituent des quasi-particules analogues aux électrons de conduction et aux trous dans les 
semiconducteurs. 

Contrairement aux photons, quanta du champ électromagnétique, les phonons ne peuvent se 
propager dans le vide ; ce sont des particules virtuelles. 

A l’équilibre thermique, les phonons suivent la statistique de Planck (bosons en 
nombre non conservé). 

• L’énergie d’un phonon du mode A est : tuo } . 

• Le nombre moyen de phonons de ce mode est : Na = ptwj ^ . 

_ 3 N 

• L’énergie totale des 3N oscillateurs : U = E = + 1/2 )tico À 

A=l V 

• Le seul problème à résoudre est de déterminer le spectre de fréquences : {co À } . 


2 °) 




Théorie d’Einstein 

Le spectre est supposé fonné uniquement d’une seule fréquence : u E = colin (fréquence 
caractéristique d’Einstein ~10 12 à 10 13 Hz) en considérant les oscillateurs indépendants, 
discernables et identiques. 


L’énergie moyenne est : 


_ f 

E = 3NTuo 

v 



- 1 


+ 1/2 . 


• La chaleur spécifique : 


C v = 3 Nk(j3hco) 2 


P tuo 


(e 


phco 


-l ) 2 










_0 e /2T 


L’énergie libre : 


F = -3NkT LogZ avec Z = 


S e /t 


1 - e 

Approximation classique des grands nombres quantiques: Tuo <tc kT 


et 9 e = lw E /k . 


C v - 3 Nk 


1 - — 

12 


r 0, ' 2 


+... 


v 1 y 


C v = 3 Nk 


(Dulong et Petit) 


Limite classique des phonons localisés : tuo » kT 


(a \ 


C v - 3Nk 


-6 e /T 


v 1 y 


Ce modèle simple permet d’expliquer pourquoi C v — r >o -> 0 ; le nombre de phonons tend 

vers zéro, de manière analogue au nombre de photons émis par un corps noir. 

Cependant, l’accord de la théorie d’Einstein avec l’expérience n’est que qualitatif. 


3°) Théorie de Debye 

• Aux très basses températures, seuls les modes de vibrations de fréquences basses (c’est à 
dire de grandes longueurs d’onde) sont excités. Quand cette longueur d’onde est beaucoup 
plus grande que les distances interatomiques, on peut considérer le cristal comme 
homogène et élastique. Dans ce milieu des ondes de déformation élastique peuvent se 
propager. Elles correspondent aux modes acoustiques. On néglige ici les modes optiques. 
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Les ondes se propageant dans une direction Ox avec la vitesse c= co/k ont une amplitude : 
, <(*•* - ® f ) J solution de l’équation d’onde : 


A(x, t) = Aq 


d 2 A 


8 2 A 


dr 


- c 


dx l 


= 0 


• Dans un volume, ayant les dimensions d’un cube de côté L, il s’établit des ondes 
stationnaires dont les modes sont donnés par les conditions aux limites habituelles : 


A(0, t ) = A(L, t ) pour chaque direction. C’est à dire : k x = n x ; k y = n y et k, = n, . 


TC 


K 


• La densité des modes longitudinaux : 

• La densité des modes transverses : 


/ \ V co 

(®) = ~ 
2.n c, 


i \ - V co 

H = 2 ITT — 

ITC C, 


en tenant compte des deux directions transverses possibles ; c’est à dire des deux directions 
de polarisation des ondes transverses (dégénérescence deux). 


• La densité totale des modes des deux types : 
en définissant la vitesse moyenne par : 




-3 

c 


1 2 

"T + ~ 
c, c, 


t 


• Le modèle de Debye est correcte pour les fréquences basses. Aux fréquences élevées 
©(&>) est différente de celle donnée par ce modèle . 

Pour tenir compte du nombre fini de modes (3 N), la densité ! 'D{oj) est arrêtée à la 


fréquence maximale telle que : 3N = J Œ>(a>)dco = 


V o), 


2n 2 


OU CO, 


= c( 6k 2 n) ! . 


La longueur d’onde associée A m = 2Kcjco m est de l’ordre des distances interatomiques. Elle 
est du même ordre que la fréquence de coupure co c . 

En effet, à cause de la structure atomique du milieu, la longueur d’onde ne peut devenir 
inférieure à la distance interatomique. Le spectre des modes présente une coupure à co c . Plus 
précisément cette limite est atteinte lorsque deux atomes voisins oscillent en opposition de 
phase. co c est l’analogue du sommet de la bande permise pour les électrons d’un solide. 



<u(27T. 10” * ') 


Fig. : Densité des fnodes ( D ( CO ) [ — p ( (O ) sur la fig] . 

C 0 D ( fréquence de Debye ) — 00 m 

Elle satisfait à la condition de Bragg : ‘’Le déphasage de l’onde entre deux atomes voisins 
est égal à une demi-longueur d’onde”. 
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Dans le modèle de Debye, l’énergie moyenne est : e -3NkT 


0 , 


\ ^ t \ 3 \ t 3 dt 

je, = 


0 D - îuo M /k est la température de Debye; elle est calculée à partir de la mesure des vitesses 
du son dans le solide : 0 D ( Pb ) = 90 K ; 0 D (AI) = 400 K et 0 D = 200 K pour le Diamant. 

3 x 

• A haute température j » q d , on a : x « 1 et A (x) — - J rdt = 1 . 

X n 

La contribution des phonons à la chaleur spécifique est alors : C v - 3 Nk . On retrouve 

la loi de Dulog-Petit. 


3 00 T 

A basse température T <s:0 D , x » 1 et A (x) — T |— -dt = 


K 

5r* 


La contribution des phonons à la chaleur spécifique devient alors : 


12 4 , 

f T 3 

3 

C v - — rc A Nk 



v 5 

\&d J 



qui est en pratique souvent vérifiée par les solides. Une seule mesure de C v 
n’importe quelle température (dans ce domaine). 


9 n et C,, à 



• Aux voisinages de T - 0K , C v varie linéairement avec la température. Ceci est dû 

aux électrons de conduction dont la contribution devient prépondérante. 

• Les valeurs expérimentales sont en bon accord avec le modèle de Debye. 

Aux températures élevées, on trouve un développement de la forme : 



1 

(d,A 

2 

3 Nk 

1 


+ ... 


20 

[t J 



formule semblable à celle obtenue d’après le modèle d’Einstein. 
• L’accord entre les deux modèles est obtenu à condition que: 



«-) i , 

1/2 

ii 

0E 


IL 
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